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FOREWORD 


This investigation was performed for the NASA Manned Spacecraft 
Center Structures and Mechanics Division. Dr. Donald M.Curry 
was the technical monitor, and Dr. Kenton D. Whitehead was the 
project manager. The study was conducted by Dr. Whitehead and 
a project team consisting of Dr. K. T. Shih - Thermodynamics, 
Mr. G. L. Getline - Dynamics, Mr. R. S. Wilson - Stress, and 
Messrs R. H. Trelease and S. T. Hitchcock - Weights/Cost 
Analysis. All work was done at the San Diego Operation of the 
Convair Aerospace Division of General Dynamics with the excep- 
tion of consultation provided by Mr. J. D. Anderson of the Fort 
Worth Operation on the acoustic fatigue computer program. Re- 
sults of the study are published in two volumes; the Final Report 
(Vol I) and User's Manual (Vol II). 
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SUMMARY 


A study was performed to assimilate and develop computational techniques for the de- 
sign optimization of thermal protection systems for the space shuttle vehicle. The 
resulting computer program was theh used to perform initial optimization and sensi- 
tivity studies on a typical thermal protection system TPS to demonstrate its application 
to the space shuttle TPS design. The program was developed in Fortran IV for Convair 
Aerospace's CDC 6400, but it was subsequently converted to the Fortran V language to 
be used on the MSC Univac 1108. Documentation for the study is reported in two-vol- 
umes -the Final Report and the User's Manual. The latter contains input instructions 
and a sample problem to illustrate use of the program. 

The major effort of the investigation consisted of the development of the computational 
techniques and programming of the subsequent methodology. The program itself was 
effected in modular fashion to allow continuing improvement and update of the perfor- 
mance prediction techniques. The program logic involves subroutines which handle the 
following basic functions: (1) a driver which calls for input, output, and communica- 

tion between program and user and between the subroutines themselves, (2) a thermo- 
dynamic analysis which includes prediction of both the aerodynamic heating rates and 
the resulting heat transfer and temperature response of the TPS, (3) a thermal stress 
analysis which predicts the internal stresses and creep rates of the TPS by a discrete 
element analysis which structurally models the TPS subject to both external forces 
due to aerodynamic pressure and thermal stresses caused by heating, (4) an acoustic 
fatigue analysis which predicts both the noise excitation due to a number of external 
sources and the fatigue life of the panel, and (5) a weights/cost analysis which deter- 
mines the weight and manufacturing cost of the system by identifying and evaluating 
these parameters for each of the TPS's components parts. In addition, a system total 
cost is predicted based on system weight and historical cost data of similar systems. 
Each of the major components of the program described above is complemented by 
other subroutines which provide specialized calculations for the analyses. 

Two basic types of input are provided, both of which are based on trajectory data. In 
the first, vehicle attitude (altitude, velocity, and angles of attack and sideslip) is input 
and external heat and pressure loads are calculated. In the second, heating rates and 
pressure loads are provided to the program as a function of time. Standard program 
output includes heating rates, temperature, and stresses for the discrete elements of 
the TPS analyzed as well as dynamic stresses and the number of stress reversals for 
the panel and its weight and cost. A panel redesign technique is included to increase 
the panel thickness to transfer mechanical loads and to increase insulation thickness 
to protect the underlying load-bearing structure. In a subsequent investigation these 
redesign iterations are being refined. 
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Optimization and sensitivity studies are performed by the user by varying panel size, 
material properties, and configuration (six different metallic panel cross-section 
geometries are provided) in a series of computer runs. The program sizes panel 
and insulation thicknesses. An optimum design is then identified as the one giving 
either minimum weight or cost as a function of the parameters being varied for the 
investigation. Sensitivity studies are performed by noting the change in system weight 
or cost due to the variation in some independent variable such as trajectory or heating 
prediction method for an optimum panel configuration. 

As the final task of this study, recommendations are made for computer program 
improvements which include new thermal protection systems (both active and passive) 
and improved computational and iterative techniques. 


xiv 



SECTION 1 


INTRODUCTION 


The objective of this contract was to develop a computer program for use in optimizing 
the thermal -structural design of the thermal protection system (TPS) for manned space- 
craft in terms of cost, weight, resue, and performance, both thermal and structural. 
The study was divided into five subtasks including documentation. 

1.1 TASK 1 - DEVELOPMENT OF THE COMPUTATIONAL PROCEDURE 

This task was based on an existing capability for performing state-of-the art prediction 
of the aero thermodynamic environment and the simultaneous structural and temperature 
response of a simply-supported metallic TPS panel. The input parameters to the pro- 
gram included either (1) a given trajectory for various body locations, or (2) specific 
pressure and heating conditions. 

The computer program is composed of a number of subroutines that are called in turn 
by the driver program. These subroutines perform the following functions: 

a. Input/output communication with the program user. 

b. Prediction of the aerothermodynamic environment with free stream conditions 
based on the 1963 Patrick Reference Atmosphere. 

c. Determination of the TPS structural temperature response including resizing of the 
insulation to satisfy system temperature constraints. 

d. Performance of a thermal stress analysis of the TPS panel cross section and panel 
resizing to preclude panel failure due to ultimate tension and compression, yielding, 
crippling, and elastic stability. 

e. Determination of TPS weight and cost per unit area, the latter including manufac- 
turing, engineering, inspection, and refurbishment costs. 

The thermal analysis is performed first. For the initial configuration, an explicit 
finite difference statement of the energy equation predicts heat transfer and temperature 
response of the panel, the insulation, and the underlying structure. If the temperature 
constraint of a part of the TPS is exceeded, insulation is resized to satisfy the con- 
straint. Next, a section of the panel is isolated and analyzed as a simply-supported 
beam by a discrete element thermal-stress prediction procedure for the steady loads 
which include thermal stresses and aerodynamic pressure; the system is resized should 
structural failure occur. A sonic fatigue analysis is performed for the panel which has 
been thus far sized for thermodynamic and aerodynamic loads. The panel is assumed 
to resonate at its fundamental frequency due to any combination of noise sources such as: 
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a turbulent boundary layer, the booster engine, and the jet (flyback) engines. The analy- 
sis of the jet engines can also include the effect of the jet engine exhaust scrubbing the 
panel of interest. The resulting applied dynamic stresses are then compared to the 
allowable stresses of the S-N data using Convair Aerospace's cumulative fatigue dam- 
age theory, and the number of stress reversals endured by the structure at the applied 
level is printed out. This can be related by the program user to the fatigue life of the 
panel by evaluating the data point in terms of the allowable S-N curve. 

The resulting TPS design is next analyzed to determine system weight and cost per unit 
area. The weights procedure adds up weights of component parts of the TPS such as 
the cover panel , insulation , substructure , and fasteners . The cost routine , in its cur- 
rent state of development, estimates TPS unit costs by two methods: (1) system manu- 
facturing costs are predicted by identifying manufacturing costs for each of the compo- 
nent parts of the TPS, and (2) theoretical first unit costs and both recurring and non- 
recurring operational costs are identified as functions of material weight and complexity. 
Results printed out in the computer program include heating rates, structural tempera- 
ture, stresses, and material property values, combined dynamic stresses and stress 
reversals, and weight and costs per unit area of the TPS and its component parts. Pro- 
gram discussion is presented in Section 2. 

1.2 TASK 2 - DESIGN OPTIMIZATION STUDIES 

The computer program developed in Task 1 was used to perform a design study for a 
typical metallic TPS configuration to demonstrate the validity of the procedure in terms 
of cost, weight, reuse, performance, and mission requirements. The panel demon- 
strated was located on the bottom centerline of the booster. The configuration chosen 
was an open semi-smooth corrugation (a corrugation wavelength of six inches and a 
depth of 0.5 inch) of Rene' 41. The temperature constraint on the system was a 300°F 
temperature limitation on the underlying cryogenic tank. The results of the analysis 
along with descriptions of the thermodynamic and thermal stress models and other per- 
tinent input data and results are shown in Section 3 . 

1.3 TASK 3 - SENSITIVITY STUDIES 

Results of the optimization studies of Task 2 were analyzed to demonstrate application 
and validity of the TPS sizing methodology. The same semi-smooth Rene' 41 panel 
with open corrugations was investigated to show the effects of heating uncertainties . 

The nominal case considered the application of conservative heating amplification fac- 
tors to account for the heating reduction and increase respectively for a separated and 
reattached boundary layer. Two off-nominal cases were considered next: (1) no heating 
amplification, i.e. , flat plate heating for an attached boundary layer, and (2) amplifica- 
tion upon the nominal case of 25 percent for turbulent heating and 10 percent for laminar 
heating. A more detailed description of this study and computed results and sensitivi- 
ties are presented in Section 4. 
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1.4 TASK 4 - PROGRAM DEMONSTRATION AND DOCUMENTATION 

The computer code was developed in. Fortran IV for the CDC 6400 at Convair Aerospace 
in San Diego. Section 5 of this report discusses a few elementary differences noted be- 
tween the Fortran IV and V programs which run on the GDC and MSC machines . The 
conversion and demonstration at MSC offered no difficulties. The program documentation 
is shown in detail in the companion volume , the User's Manual. 

1.5 TASK 5 - RECOMMENDATIONS AND PROGRAM REFINEMENTS 

The following paragraphs summarize recommendations for an extension to the work of 
the present contract. Details are presented on Section 6. The objective of a continued 
program is (1) to refine development of a computer program for use in optimizing the 
thermal-structural design of the thermal protection system of the space shuttle vehicle, 
and (2) to employ both the computer program developed under Contract NAS 9-10956 and 
an improved version to generate parametric weight and cost data for optimum TPS for 
a variety of materials, panel and support configurations, vehicle locations, and vehicle 
trajectories. This investigation is divided into four tasks, the first two of which up- 
date the computer program to include new TPS concepts and to improve computational 
speed, the third is the design of optimum TPS at local areas on the shuttle vehicle, and 
the fourth is documentation. 

1.5.1 COMPUTER PROGRAM IMPROVEMENTS . The primary function of the task is 
to employ refined numerical and mathematical models to decrease computer run time 
and make application of the computer program for design optimization easier for the 
program user. The mathematical and numerical methods improvements include refine- 
ments in varying degrees of difficulty to all of the major subroutines of the analysis: 
thermodynamic, stress, fatigue, weight, and cost. In addition, current TPS concepts 
on the space shuttle will be included in the sizing program. These include reusable 
surface insulation, carbon-carbon systems, mass transfer cooling, and ablators. 

Improved thermodynamic analysis involves incorporation of newly developed aerodyna- 
mic heating prediction at high angles of attack including real gas effects, cross-flow, 
and boundary layer transition. Internally, the adoption of any one of a number of im- 
plicit heat transfer equations presently available at Convair Aerospace is recommended 
after a short evaluation study has been conducted of a few typical configurations and 
problems. The major changes and additions proposed for the TPS stress analysis in- 
clude reorganization and refinement of existing analyses of the metallic panels and their 
support, computation of stress redistributions due to creep and panel deflections, and 
development of additional models for reusable surface insulation, carbon -carbon sys- 
tems, mass transfer cooling , and ablators . Fatigue analysis refinements will include 
an assessment of the effects of damping of the panel due to edge constraints and edge 
members, supporting structure, and the insulation. Also a low frequency analysis 
(panel flutter and stability) will be adopted to include the effects of this structural re- 
quirement on TPS redesign. 
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Weights and cost analyses will increase the depth of their considerations to include 
variations in manufacturing and development techniques and newly developed data on 
refurbishment costs. Mathematical models will be developed to predict weight and cost 
of the new TPS's in sufficient depth to allow valid comparisons of these important pa- 
rameters. Finally, an intensive effort will be made in a search of the literature to ob- 
tain better material property data; this will include the effects of material degradation 
with time. Also, whatever new data (either S-N fatigue data, or weight/cost data) are 
developed and published will be included in a form for quick and ready access to the 
computer input. 

1.5.2 PROGRAM APPLICATIONS AND OPTIMIZATION TECHNIQUES. This task will 
be divided into two subtasks: short term improvements and formal optimization proce- 
dures. In the former, techniques will be evaluated and employed (if practicable) to 
automate the computational procedures. These might include the stacking of succes- 
sive cases, the curve fitting of resulting weights and costs, and minimization of the 
resulting curves to determine the optimum configuration. In the latter task, MSC's 
previous experience in the application of optimization procedures to the thermal sizing 
of TPS will be assessed to determine if there is a rational but economically feasible 
approach to combine the TPS sizing program with a formal optimization procedure . 
Finally, additional TPS designs, possibly including active cooling systems for localized 
areas and probably leading edge concepts, will be incorporated into the size procedure. 

1.5.3 OPTIMIZATION AND SENSITIVITY STUDIES. Detailed optimization and sensi- 
tivity studies of configuration and materials already available to the program will be 
begun at program go-ahead. The parametric studies will consider typical vehicle areas, 
materials, and trajectories and will in fact be an organized concerted effort to size the 
shuttle TPS locally. Particular attention will be paid to one vehicle area and TPS con- 
cept to gather information on the time and costs of running a detailed sizing study with 
the techniques developed under this contract. The remaining study time can then be 
scheduled to develop valid parametric data in a cost effective manner. 

1.5.4 PROGRAM DEMONSTRATION AND DOCUMENTATION . The TPS Optimization 
Computer Program will be developed at Convair Aerospace in Fortran IV for the CDC 
6400. The resulting program will then be checked out on the Univac 1108 in Fortran V. 

Early in the program development, programming instructions will be established to en- 
sure a minimum number of changes both in computer languages and systems. The 
documentation of the User's Manual will include all mathematical descriptions, methods 
of solution, program listings, flow charts, list of symbols, diagnostic messages for 
typical failures, and sample problems. Documentation of the entire study will include, 
in addition, results of the optimization and sensitivity studies and all conclusions and 
recommendations . 
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SECTION 2 


COMPUTER PROGRAM 

This chapter describes the prediction methods and analyses adapted for use in the TPS 
sizing computer program. The complete listing and operating instructions for the 
program are presented in a companion volume entitled User's Manual. 

2. 1 PROGRAM ORGANIZATION 

The thermal protection system sizing program consists of a number of modules or sub- 
routines each of which is designed to perform a specific function of design and/or com- 
munication with other subroutines as well as the main or driver program. The 
major modules perform the tasks of program input, program output, thermody- 
namic, stress, sonic fatigue analyses, thermodynamic and structural redesign, and 
weight and cost evaluations. A number of smaller subroutines have been developed to 
perform specific auxiliary tasks for each of the major analyses. For example, these 
smaller problems include nodal specification of the stress analysis model, determina- 
tion of applied and allowable dynamic stresses for the acoustic fatigue analysis, and 
storage of various data blocks in the weight/cost subroutine to name just a few. 

The following sections describe each of the major analyses alluded to above with the 
discussion of input and output (including an extensive description of a number of thermo- 
structural properties required for program operation) presented in Section 2. 8. 

2.2 . AEROTHERMODYNAMIC ENVIRONMENT 

The prediction of aerothermodynamic environment is based on Convair Aerodynamic/ 
Structural Heating Program 3020 (References 1 and 2). The local environment is 
established by the prediction of the flow field about simple geometric shapes such as 
flat plates, wedges, cones, cylinders, or spheres. 

Prediction of aeroheating is classified into two regimes: high and low local angle of 
attack. For low angle -of-attack applications, the shock waves are assumed attached 
to the body, and flow field properties are computed from tangent wedge/cone techniques. 
Using these local properties, the algorithm then computes local heating rates using 
either the Eckert reference enthalpy method or the Spalding -Chi technique. Transitional 
heating between the laminar and turbulent boundary layers will then be calculated as a 
linear interpolation of turbulent and laminar heating values, the degree of turbulence 
depending on the "turbulent fraction" exhibited by the boundary layer with respect to 
values of Reynolds number for transition onset and end. 

At high angles of attack, the flow field cannot be predicted so conveniently as at low 
angles of attack. Thus, current state-of-the-art techniques recommend aeroheating 
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rate calculation by swept cylinder methods, either laminar or turbulent. At the 
moment, no transition criterion has been established for the switch from laminar to 
turbulent swept cylinder heating prediction techniques. 

Detailed descriptions of the aeroheating and (where applicable) pressure prediction 
equations are given in the following: 

2.2.1 F LOW FIELD CALCULATIONS. The effective angle of attack, t is defined as 

sin r = DNX cos a cos 0 + DNY sin $ + DNZ cos 0 sin a (2-1) 

where 

a = angle of attack 
j3 = yaw angle 

DNX, DNY, DNZ = direction cosines of outer normal from surface, and the 
high angle of attack is defined arbitrarily for t > 35°. 

2.2. 1. 1 Low Angle of Attack. The thermodynamic properties of the air ahead of the 
shock, P ro , To,, and and the free stream Mach number M^, are determined from 
the given flight conditions using the 1963 Patrick AFB atmosphere. Once T^, and the 
surface temperature T w are known, the free stream and wall enthalpies can be deter- 
mined using Figure 2-1. A set of empirical equations (Reference 3), shown in Figures 
2-2 through 2-7, is used in determining the shock layer thermodynamic properties at 
the boundary layer edge as functions of the hypersonic similarity parameters Moo sin a 
and Mj.a (Reference 4). For the flat plate configuration, a = 0 degrees, it is assumed 
that the free stream conditions exist at the boundary layer edge . 

2.2. 1.2 High Angle of Attack. At high angle of attack the pressure is calculated by 
a method developed by Convair Aerodynamics Group (Reference 5). Other properties 
are not available for the present program. 

The method which was developed is a modified flat plate approach in that the pressure 
at a point depends only on Mach number and on the local slope of that point. However, 
the local slope depends not only on the geometry of the point, angle of attack, and yaw 
but also is corrected for a boundary layer displacement term which accounts for wall 
temperature, Mach number, unit Reynolds number, and distance from a leading edge. 

Thus the pressure model contains three main parts: 

a. Boundary layer displacement correction. 

b. Windward pressure model. 

c. Leeward pressure model. 


2-2 



ENTHALPY (Btu/lb) 


i = 2.2864 + 0.2274T + 0.11571x10 




REFERENCE 

TABLES OF THERMAL PROPERTIES OF 
GASES (N. B. S. CIRCULAR 564) 


TEMPERATURE ( R) 


Figure 2-1. Enthalpy Vs. Temperature for Air 


Boundary Layer Displacement Correction . Failure to derive correlations of pressure 
using combinations of impact pressure laws, viscous induced pressure laws, and blast 
induced pressure laws led to the belief that the pressure data could be correlated by a 
common pressure law in which the effect of the viscous terms is accounted for in a 
revised local slope. Thus: 


where 


K (M, r ' ) 


= pressure coefficient 


M = Mach number 


K(M, T 1 ) = unknown pressure law 
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ENTHALPY RATIO 


EQUATIONS : 



Figure 2-3. Wedge Enthalpy Ratio Vs. Sin a 
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M sin a 
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Figure 2-4. Wedge Velocity Parameter Vs. Sina 
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EQUATIONS 


Mo, sin q; < 1*5 : P /P^ = 1.007 + 0.3816 Moo sin a + 1.522 (Moo sin <*) 

8 - 0.1593 (Moo sin a) 3 

1.5 < Moo sin a < 5 *° : p B / p oo = 0*2397 + 1.161 Moo p sin a - 

+ 1.06 (Moo sin ar + 0.0489 (Moo sin a) 

o 

5.0 £ Moo sin a : P /V = -3.182 + 4.177 Moo sin a + 0.8373 (Moo sin ar) 

8 + 0.0216 (Moc sin a)3 


1000 


100 


Pii 

\ 



0 

.5 

1.0 

1.5 

2.0 

2.5 

3.0 

0 

1.0 

2.0 

3.0 

4.0 

5.0 

6.0 

0 

5.0 

10. 

15. 

20. 

25. 

30. 




Moo sin a 





Figure 2-5. Cone Pressure Ratio Vs. Sin a 
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ENTHALPY RATIO 


EQUATIONS : 

2 

Moo sin a - 8 *^ : * /ioo = 1.03 + 0.0827 Moo sin a + 0.2354 (M m sin a) 

S - 0.6956 x 10~ 3 (Moo sin a) 3 

o 

Mqo sina > 8.0 : i /i^ = 1.106 - 0.3685 Moo s * n a + 0.3466 (Moo sin a ) 

8 - 0.007766 (Moo sin “) 3 
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Figure 2-6. Cone Enthalpy Ratio Vs. Sina 
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M = M„ x 10 

OO 


-3 


EQUATIONS: 


0 < sin a < 4.4 : (V 2 - V 2 ) x 10 ® = 0.157 + 0.75 M r sin a 

+ 0.9861 (M sin a ) 2 + 0.06944 (M sin a)^ 
r r 

4.4 < M r sin a : (V 2 - V 2 ) x 10' 6 = 6.187 - 1.038 M r sina 

+ 1.414 (Mr sin a) 2 - 0.0062 (Mr sin a)' 5 



M sin a 
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Figure 2-7. Cone Velocity Parameter Vs. M r Sina 
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and 


r ’ = T + At (approximately) (2-2) 

At = viscous effect correction 

Equation 2-2 is an approximation because the boundary layer correction should be ap- 
plied to the direction cosines of the outer normal at the point in question and the 
correct expression for corrected local slope is 


CST = 

cos (A T ) 


(2-3) 

SST = 

sin (At) 

DNY 

(2-4) 

^/(DNZ 2 + DNY 2 ) 

BST = 

sin (AT) 

DNY 

(2-5) 

J( DNZ 2 + DNY 2 ) 

DNX* 

= DNX * 

CST + DNZ • SST \ 


DNZ' 

= DNZ * 

CST - DNX • SST > 

(2-6) 

DNY' 

= DNY • 

CST - DNZ * BST ) 


sin t ' 

= coso! cos 0 DNX' + sin0 DNY' + cos j3 sin a: DNZ' 

(2-7) 


where Equations 2-3 to 2-5 are used to give the proper magnitude and sign to the Y 
and Z components of the boundary layer corrections. 

The correction was derived from the Chapman-Rube sin laminar displacement thickness 
expression for a flat plate in the following way: 


/Be T 

6 V CX = 1,721 + 0,322 (y ' 


(2-8) 


and 


AT = tan 


- m') 


(2-9) 


A T = tan 


-1 


1 / C W 2 -1/2 

I\/r~ a- 721 — + °- 0332 (r - 1) M m ) x 
V el 03 
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where 


Re i = unit Reynolds number 

C = Chapman-Rube sin constant 

T = wall temperature 

w 

T = ambient temperature 

M = free stream Mach number 

X> 

X = distance from leading edge 

A turbulent correction has not been developed yet so that the laminar correction is 
used at all conditions at present. 

An arbitrary limit is presently applied to prevent this term from generating very large 
corrections at the leading edges of the vehicle (where X goes to zero) so that the 
largest correction which can be made is 30 degrees. This correction also causes some 
points to have positive slopes where geometrically they have negative slopes and thus 
can result in positive pressure on the leeward side of vehicles which agrees with test 
data. 


Windward Pressure Model. An equivalent body shape has been defined and the correc- 
ted local slope at the point in question is positive so that it is a windward point. A 
modified tangent wedge pressure formula was derived (Reference 5). 


, 2 

sin (t 




where 


K = M sin r ' 

2 2 
0 = M - 1 

T 1 = corrected local slope 


Cp = pressure coefficient 

y = ratio of specific heats 

M = free stream Mach number 
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This pressure model predicts pressures very close to Van Dykes tangent wedge at 
small inclinations to the stream and closely follows modified Newtonian at large 
inclinations to the stream. 

Leeward Pressure Model. The leeward pressure model is composed of three parts: 

a. Small disturbance theory law. 

b. Two dimensional base pressures. 

c . Limit turning angle . 

The small disturbance theory pressure formula is used to predict pressures on surfaces 
from zero inclination to the stream to negative inclinations defined by a limit turning 
angle. The small disturbance theory law is 



. 2 

sin r' 





K = j3 sinr' 


( 2 - 11 ) 


The limit turning angle is defined from Love' s limit turning angle of Reference 6 and 
is shown in Figure 2-8. Two-dimensional base pressures are used on surfaces which 
have larger negative inclinations than the limit turning angle and the pressure is com- 
puted using the following empirical fit: 


C p = (-0.3008/ (M-0. 5434)) + 0.01132M -0.05252 
for Mach numbers less than 6 and 


C 

P 


-1 


0.7 M 


2 


(2-12) 


(2-13) 


for Mach numbers greater than 6. 

2.2.2 FLAT PLATE, WEDGE, AND CONE AERODYNAMIC HEATING. At low angle 
of attack, the flow field properties are computed by tangent wedge/cone techniques as 
described in Section 2.2.1; aerodynamic heating rates are computed in two combina- 
tions: 

a. Eckert laminar with Eckert turbulent. 

b. Eckert laminar with Spalding -Chi turbulent. 

Eckert reference enthalpy method (Reference 7) depends upon the assumption that the 
incompressible mass, momentum, and energy equations can be used for compressible 
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CRITICAL TURNING ANGLE, (deg) 



Figure 2-8. Limit Turning Angle 

flow solutions provided the thermodynamic and transport properties of the gas are 
known and are evaluated at an appropriate reference enthalpy. Real gas effects, 
including dissociation, are taken into account in the determination of the properties 
of the gas just outside the boundary layer. The Eckert reference enthalpy is 

i* = i + 0.5 (i - i ) + 0.22 (i - i ) (2- 

s ws rs v 


where i g is the shock layer (just outside the boundary layer) static enthalpy, ^ is the 
enthalpy of air at the wall temperature, and i r is the boundary layer recovery enthalpy, 
given by 


i 

r 


i + r 
s 


V 


2g c J 


(2-15) 


where the flow recovery factor, r, is 0. 84 for laminar flow and 0.89 for turbulent 
flow. 
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2.2.2. 1 Laminar Boundary Layer. The Blasius solution for an incompressible 
laminar boundary layer gives a skin friction coefficient 


0.664 

f yifc* 


Subsitution into the Reynolds analog relation 


Pr* 2//3 St 



(2-16) 


(2-17) 


yields the solution for Stanton number 
_ 0.332 

y^Re* Pr *2/3 

where the Stanton number is defined by 


h 

i 

p* V 

s c 


(2-18) 


(2-19) 


and the Prandtl and Reynolds numbers are 


Pr * = 


u* C * 
_ P_ 

k* 


Re* = 


p* V X 
s 


(2-20) 


( 2 - 21 ) 


A constant value of 0.71 is assumed for Pr*, based on the data of Hansen (Reference 8), 
as a convenient approximation. The heat transfer coefficient is defined by 



(2-22) 


Using the relationships given above, the heat transfer coefficient can be expressed as 
13.43 Vp*p* V 

h = p= (2-23) 

i V X 

The heat transfer coefficient, h^, is multiplied by the factor 1. 73 for conical flow 
(Reference 9). 
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2 . 2 . 2 . 2 Turbulent Boundary Layer 


Eckert's Reference Enthalpy Method. The Schultz -Grunow solution for the local skin 
friction coefficient is (Reference 10). 


7.485 p* V 


c f - 


(log 1Q Re*) 


s 

2.584 


(2-25) 


For conical flow (Reference 11), the heat transfer coefficient is multiplied by the 
factor 1. 176. 


Spalding -Chi Heating and Drag Technique ■ The original intent of the Spalding -Chi 
method was to develop a procedure whereby skin friction for the compressible flow 
over a flat plate could be computed rapidly by hand (Reference 12). The theoretical 
development is based upon the postulate that, for incompressible flow, skin friction 
can be expressed functionally in terms of Reynolds number based on momentum thick- 
ness as Cf = 'I'g (Re 0 ) where the tilde denotes incompressible quantities. The relation- 
ship between incompressible and compressible quantities is such that 

c = c, F and Re = F _ Re rt 
f f c 9 r0 0 


where the transformation variables F c and F f0 are functions of temperature ratio 
(free stream to wall) and Mach number. (The values of these parameters approach 
unity for incompressible flow, i.e., Tw/Te = 1, Me = 0). Another expression can 
be derived theoretically to yield 1/2 Cj F c = (F tx Re x ) with a new coefficient F rx . 
The transformation coefficients were evaluated from equations suggested by previously 
successful correlation techniques. Hence, the expression 
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is evaluated across the isobaric boundary layer using the Crocco relation and the 
perfect gas law. Similarly, it was postulated that 



(2-27) 


and the exponents were evaluated for a wide variety of experimental data to be p = 

-0. 702 and q = 0. 722. The parameters F c , F r0 , and F rx were tabulated as functions 
of Mach number and temperature ratio. The computational technique is comprised of 
evaluating the transformation parameters for certain flight conditions and then, with 
the proper Reynolds number, computing the compressible skin friction coefficient. 
Application of a suitable Reynolds analogy factor leads to calculation of the film heating 
coefficient. 0 



Wallace (Reference 13) suggested that the technique could be improved for high free 
stream enthalpy by employing enthalpy instead of temperature ratios to evaluate F r g. 
With a view toward computer application, Komar (Reference 14) avoided the lengthy 
double interpolations necessary to determine F c , F rx , and F r g by curve fitting the 
product Cf F c as a function of the products F rx Re x and F r g Reg. The incompressible 
skin friction coefficient Cf F c is computed as 


c 


f 


F 

c 




(2-28) 


where 


A = log (F Re ) 
e rx x 


(2-29) 


The equations necessary to compute the parameter F are outlined below (Reference 
14). Thus C 


F = 
c 


i i L- 1 

- (2 c + b) 

-1 

- sin 

-b 

1 ^ ^ sin 

-J b 2 -4ac . 

. J b 2 - 4ac. 


-2 


(2-30) 


where 


a = 


w 


e 


and 


r 2 w 

b = 1+- (y-1) M - — ; — 
2 e i 

e 


(2-31) 


r 2 

c = “ J (y - !) 


Then 


,-0.702 


F n = 
r0 \ l 


w 
e ' 


,0.72 


aw 


w 


(2-32) 


F = F /F 
rx rg c 


A = = ^g < F rx ^ 


(2-33) 

(2-34) 
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n 


°f F 0 = exp (£ g j 

' 3-1 


j-A 


(2-35) 


The coefficients gj are those developed by Komar for the curve fits. 
g l : 
g 4 : 

e„ : 


g 


10 


9.2809 

S 2 

= -4.7340 

g 3 = 

6.6859 x 10 -1 

-2 

-4. 1877 x 10 

g 5 

= -5.5055X10 -4 

g 6 = 

2.8367 xlO -4 

-2. 1250 x 10~ 5 
-10 

1.3236 x 10 

g 8 

-7 

= 8.0162x10 

g 9 = 

-1. 5901 x 10~ 8 


Successive calculation of the parameters of Equations 2-31, 2-30, 2-32, 2-33, 2-29, 
and 2-28 gives Cf. Reynolds analogy leads to the heating coefficient. A recommended 
equation for the Reynolds analogy factor is that of Karman as proposed by Bertram 
(Reference 15). 


S 



5 Pr + 1 \ 
6 ) 


(2-36) 


The only parameter still undertermined in this analysis is the shock layer Mach 
number. For a Prandtl- Meyer expansion; this computation was already made but 
not printed as output. For oblique shock waves, the local speed of sound is com- 
puted from an approximation (Figure 2-9) of the data of Hansen (Reference 8). The 
sound speed is computed from the relation 


2 

t-Re - m T + b 
P e 

for T < 2700°K, m=-8.75 (10)~ 4 /°K, b = 1.432 

for T > 2700°K, m = 0, b = 1.2 

2.2.2. 3 Transitional Boundary Layer. The transitional boundary layer is represen- 
ted as a linear transition from laminar into turbulent flow. Transition is assumed to 
begin at a specified value of the shock layer Reynolds number, denoted Re^ r , and end 
at a specified value of the shock layer Reynolds number, denoted Reg. The heating 
rate is computed as the linear interpolation between laminar and turbulent values, the 
virtual origin of both being the leading edge. Interpolation is performed using the 
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Figure 2-9. Speed of Sound Parameter for Air as a 
Function of Temperature 

value of Reynolds number that occurs between the input values specifying the onset 
and end of transition. Hence, for example, the film heating coefficient is given by 


Re - Re 


h = 


tr 


Re - Re 
E tr 


h + 
T 


Re - Re. 


1 - 


tr 


Re„ - Re^ 
E tr 


* f h '+ f h T 
T T L L 


6 6 

The values usually used are Re, - 10 and Re_ = 2 x 10 such that 

tr E 


(2-38) 


f 

tr 




2 


(2-39) 


These values are reasonable (Figures 2-10 and 2-11) according to Masaki and Yakura 
(Reference 16). 

2.2.3 SPHERE AND CYLINDER AERODYNAMIC HEATING. At high angle of attack, 
swept cylinder theories are used to predict the aerodynamic heating. For laminar 
flow, heating is computed by method of Kemp and Riddell, while the method of Beck- 
with and Gallagher is used to compute the turbulent heating. 
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Figure 2-11. Effect of Mach Number on 
the Transition Zone Reyn- 
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2.2.3. 1 Laminar Swept Cylinder. Laminar heat transfer to a swept cylinder is 
computed through a transformation of spherical heating rates. Heat transfer to a 
sphere is calculated by the Kemp-Riddell (Reference 17) expression 
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Cylindrical heating rates are then obtained by adjusting for sweep by 

q cyl = °- 75 q sph (OOSA)1 ' 2 


(2-41) 


where the sweep angle 


A = 90° - r 

2 . 2 . 3 . 2 Turbulent Swept Cylinder. Turbulent heat transfer to a cylinder is com' 
puted by the method of Beckwick and Gallagher (Reference 18). Equation for heat 
transfer to the stagnation line of a swept cylinder is 
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where a and n are constant in the Blasius skin friction law and were taken as a = 0.0228 
and n = 4. Ur is the free stream velocity vector and u is the component of Ut, 

1X 9 CO CO -tV> OD 

normal to the cylinder stagnation line; thus 


u = U R (cos A) 

CO AV I CO 


(2-44) 


Using this relationship and the definition of Stanton number gives 
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(2-45) 


Sutherland's equation for the viscosity of air is used. Static pressure at the wall, p , 

b 

is that which would be sensed by a pitot tube placed normal to the bow shock when u^ 
is supersonic. This pressure is 


P 


s 


where 


7 7 , 

t ■■ ■ ■ r -i -■ ^ 


y+ 1 M 2 

7-1 

7 + 1 

7-1 

2 N, “ 


27 M 2 - (7 “ 1 ). 
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M„ = M m cos A. 
N, 00 00 


The stagnation line velocity gradient can be expressed as 



1/2 


(2-46) 


(2-47) 
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where 


2 

a N 0 

T = -1-- (2-48) 

s y R 


and a N is the normal component of the free stream stagnation sonic velocity (Ref- 
ence 1^); therefore 


T = T 
s 00 



(2-49) 


Reference conditions T and u were evaluated at T in Equation 2-45. 

r r s 


The final expression for heat transfer to the stagnation line was given as 


q cyl, T 



(2-50) 


where values of recovery factor used in evaluating i r were obtained by curve fit of the 
data in Beckwith and Gallagher (Figure 2-12). 


(a) A = 0" 



Figure 2-12. Variation of Recovery Temperature Around Circumference 
of Cylinder at Station 2 
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2 . 2 . 3 . 3 Correction for Flat Face Velocity Gradients. The methods of Bertram and 
Henderson (Reference 19) are used to correct the flat face velocity gradients. Since 
stagnation heating is proportional to the square root of free stream velocity gradient, 
it was postulated that 


q 

q 



I d Ue 



flat 

A dx J 

1 x = 0 

flat 

cyl ^ 

Jm 

^ x = 0 

cyl 


(2-51) 


For laminar boundary layer of a flat disc 
(Figure 2-13) 

(ir)x = 0 = 0.745 + 3.14(^) (2-52, 

the laminar centerline heating rates are with the form 


Figure 2-13. Geometry of Flat Face 
Velocity Gradient 

Correction ilat cyl 

For turbulent boundary layer a 1/5 power correction was used. 


'0. 745 + 3 


2.315 


■“«) 


(2-53) 
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= q 
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745 + 3. 14 (r/D)1 1/5 


2.315 


■J 


(2-54) 


2 . 2 . 3. 4 Transport Properties . The transport properties given by Hansen (Reference 8) 
are shown in Figures 2-14 and 2-15. From these, curve fit equations were obtained for 
p and fJ. as functions of enthalpy and pressure: 


i < 1300, 

-4 _ .-0.849 
p = 0.576 x 10 P i 

s 

(2-53) 


H = 0.42642 x 10" 7 i°‘ 493 

(2-54) 

i * 1300, 

p = 0.865 xlO- 5 P i -0 - 584 

s 

(2-55) 


-6 0.228 

fi = 0.28428x10 l 

(2-56) 


where base values of and evaluated at T^ = 400R were used as reference values in 
the jj,/pk and p |j/p b ^ ratios shown on the figures. The base value of density pj> is 
adjusted as a function of the local static pressure. 


2-22 



VISCOSITY RATIO (v/fi.) 



Figure 2-14. Viscosity Ratio Vs. Enthalpy 
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DENSITY -VISCOSITY PRODUCT (PfVP M- ) 



Figure 2-15. Density -Viscosity Product Vs. Enthalpy 
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Sutherland’s formula is used to compute the viscosity at free stream conditions. 




2.27 x 10 


t 3/2 

T + 198.6 


(2-57) 


2.3 STRUCTURAL TEMPERATURE DETERMINATION 

The temperature at any point in the structure is a function of external and internal 
convection and radiation rates and the conductivity and thermal inertia properties of 
the structure itself. In this program, structural temperature distributions are 
evaluated through use of the lumped parameter method of finite differences. 

The surface and/or structure is divided into an arbitrary number of small segments. 
The segments are arranged in rows parallel and columns perpendicular to the surface 
as shown in Figure 2-16. 


x 


Y 


AERODYNAMICALLY HEATED 



INNER SURFACE 


Figure 2-16. Surface/Structure Segmentation 

Then, for some small time increment (At), the net heat flux to each surface segment 
is determined. For a segment i, then, the temperature change from time to time, 
t + At, is 

T i (t + At)-T i ,t) - [Vt A BL- (<)A) INSD- S(kA f 1 'n' S (<l rad A >d 


At 

(WC ). 

pi 

where 


A = area perpendicular to direction N 
N 

k = effective thermal conductivity in direction N 
N 


(2-58) 
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X N = length of conduction path in direction N 

AT^ = temperature difference between adjacent elements in direction N 
at time t 

q , = heat transfer by radiation between nodes 
rad 

W = weight of element i 


C = specific heat of element i 
P 

It is assumed that all the mass is concentrated in a point at the centroid of the segment, 
A n is equal to the segment interface area, and is equal to the distance between 
centroids in direction N. 

The net heat transfer is given by 



where q A is the boundary layer convective heat transfer rate, q rad represents the 
energy loss due to surface radiation, and is an optional multiplying factor. 

This factor may be used to approximately allow for the effects of shock wave inter- 
actions, flow divergence, etc. 

The term (qA)jNSD provides an optional capability to include internal convective 
cooling of the structure. The cooling occurs on the backface of the last segment. 

Terms in Equation 2-58 which do not apply to a given element are dropped out for 
that element. Thus, the q net term applies only to elements of the first row, i.e. , 
those representing the surface, and the (qA) INgD term applies only to elements of 
the last row, i.e. those representing the backface. Internal radiation heat transfer 
is also taken into account in this program. The accuracy of Equation 2-58 is depen- 
dent upon the size of the segments and the computation interval and improves as 
these parameters are decreased. The program may be used for either one -dimensional 
or two-dimensional arrays of segments. 

In the following, the two-dimensional case is discussed; the one -dimensional problem 
is treated as a special case of a single column. The structure is set up in a matrix 
shown in Figure 2-17. In addition to the rectangular coordinate system, the cylindrical 
coordinate is also included as an option. The material may not be homogeneous in 
either the x-direction or the y-direction. The maximum number of rows and columns 
is nine each. Configurations are given in Figure 2-18, and nodal breakdowns for the 
thermodynamic and stress analyses are given in Figure 2-19. 
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Figure 2-17. Structure Heat Transfer Matrix 


The temperature at any segment at time t + At is dependent upon the summation of 
direction heat transfer rates at time t. The general equation is given by Equation 
2-58. From Figure 2-17, for atypical node, the temperature, t, changes from 
time t to t + At by 


T (t + At) 
ij 


■ t ij (t) + 


EL e 


IJ 


Q 

k=l 


(2-60) 


Aij is the thermal mass of the node as defined in Table 2-1, together with other 

parameters. The equations for computing Q 's are listed in the following: 

K 

a. Conduction heat transfer 


Q = (T. . - T. ) / E. 

1 i-1, j i, i-1, j 

2 v i+1, j ij i, j 


(2-61) 

(2-62) 


2-27 






2. FLAT CORRUGATION WITH SKIN 



— »-j |-« — T C 

3. RIB-STIFFENED PANEL 



4„ SKIN-STRINGER 



5, OPEN CORRUGATION 6. OPEN CORRUGATION 

(CIRCULAR ARC CORRUGATION) 


Figure 2-18. Panel Geometries 
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Figure 2-19. Configurations for Thermodynamic and Stress Analysis 

2-29 

















Table 2-1. Notations for Two Coordinate Systems 


Symbol 

Coordinate System 

Rectangular Cylindrical 

x' 

l 

X 

i 

" < r i - r i 2 -i> 

y' . 

i. 3 

y j 

2tt r. y. 
1 3 

A 

i, i 

< P Vij x i' 

B 

i, J 

X / 

In [ r /(r. - 0.5 xi)J 

2K..y. 
13 3 

2 TT K.. y. 
13 3 

C. 

i, 3 

y. 

3 

2 K. . x. 

13 1 

D. . 
1,3 


In l(r. -0.5 xi)/r. 1 
1 1 l-lJ 

2 V K.. y. 

13 3 

E. . 

1,3 

B + B 
ij i+1, 3 

B. . + D , . 

1 , 3 i+1, 3 

F 

i, 3 

c. + c 

1 , 3 i, 3+1 


Q = (T. . - -T .) /F. . 

3 M-l ij i, J-l 

Q = (T. . - T. .) / F. . 

4 l, J+l i, 3 i, J 


Radiation heat transfer 

Q = aSLJ U (T 4 -T 4 
K a,b m,n i,j) 

where 


L = x.' (for k = 3,4) 


= (for k = 1,2) 


(2-63) 

(2-64) 


(2-65) 
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c. Nodes on the first column 


Q = 0 

1 

d. Nodes on the last column 



= 0 


e. Nodes on the first row 



x 1 q 

i net. 

l 


f. Nodes on the last row 


Q = h x* (T - T ) 
4 ins l ins ij 


(2-66 


(2-67) 


( 2 - 68 ) 


(2-69) 


2.4 STRUCTURAL ANALYSIS 

This discussion covers the short-time static strength and creep analysis of TPS panels 
and support members. The analysis is restricted to simply -supported panels with ex- 
pansion joints which permit free thermal expansion. The loadings considered are 
bending due to aerodynamic pressure, and the internal forces induced by temperature 
gradients within the panel cross section. 

Given the instantaneous normal pressure and temperature distribution on the TPS panel 
at each of a series of times throughout a flight trajectory, the internal stresses are 
determined by stress analysis and hence the instantaneous creep rate for each element 
comprising the cross section is obtained from Larson-Miller curves. To determine the 
critical trajectory point or time for each mode of failure considered in the static strength 
analysis, a number of structural indices relating the applied and allowable stresses 
throughout the cross section are computed. The critical time for each mode of failure 
is determined upon completion of the trajectory by selecting the maximum value of this 
appropriate index. 

The failure modes considered are yielding, ultimate tension, ultimate compression and/ 
or crippling of the various elements comprising each panel cross section. Margins of 
safety are computed at the critical trajectory points for each of these modes. 

The total creep strain accumulated in each element of the panel cross section during 
the trajectory is obtained by the summation of the product of the average creep rate in 
each time interval and the time increment. 
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The final step in the analysis is a redesign procedure which increments the panel 
thickness in the event that the minimum margin of safety is negative and/or the maxi- 
mum creep strain exceeds the permissible value. Each step in the analysis procedure 
is further discussed in the following paragraphs. 

2.4.1 LOADS ANALYSIS AND NODAL BREAKDOWN OF PANEL. The TPS panels 
considered in this investigation are of sheet stringer or corrugated construction (both 
exposed and covered corrugations) with simply-supported edges and joints which permit 
free thermal expansion. All configurations are illustrated in Figure 2-18. Each stif- 
fening element or corrugation in the panel is assumed to behave as a simply -supported 
beam subjected to normal pressure loading and to a self-equilibrating system of internal 
loads due to temperature gradients across the section. Figure 2-20 illustrates a cor- 
rugated section attached to support configuration number one in which the panel rests 



DISTRIBUTED LOAD 
W = p X A $ 


BEAM SHEAR 

S MAX = 0-SxWxLp 


W lb/in. 


^ 1 v l w i r j r y w 


A 


A 


LOADING 


a MAX 


flTTw 


SHEAR 


AR 



Figure 2-20. Panel Loads Analysis 
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on rails which in turn are supported by heat posts. Shear stresses and bending mo- 
ments resulting from the loading are also shown. Section properties and the nodal 
breakdown are computed first. As an illustration, panel configuration number one is 
shown in Figure 2-21. Since the section is symmetrical about its centerline, only half 
of it need be considered. For n-segments of the area under consideration, the section 
area, controid, and moments of inertia are given respectively by the general expressions 

n 

A = L As. t. 
i=i i i 





Figure 2-21. Geometry of Configuration No. 1 
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zz 
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xz 
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Geometry is computed for this particular configuration more specifically in the follow- 
ing manner. First the angle 0 is computed as 


0 = cos 


■‘(¥1 


Next, the areas of the first five segments of half the cross section are determined as 


A. = — - (A -B_) t i = 1, • • * , 5 
l 10 v s F s 

and their orientation in terms of the lateral and vertical coordinates x and z is given 
respectively by the distances to the centroid 


X. = 


X + 2 X 
i-1 1 


Z. = Z 


i = 2, ' ' 5 


where 


X, = 0.05 (A - B ) 

1 s F 


Z, = 1/2 t 

1 s 
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A g and Bp are the width of the full corrugation and flange, not just half widths. The 
sixth and seventh nodes consider the combined thicknesses of the skin such that 

A. = 1/4 B_ (t + 1 ) 

i F v s c ' 


Note that the portion of the panel being analyzed extends from centerline of the corruga- 
tion to the centerline of the flange areas. The lateral distances are then determined as 


X 6 

= — (A - B„ ) 

+ 8 B F 2 A S “ 9 B F 

Z 6 

= -7 (t + t ) 
2 s c ’ 


and 



X 

7 

= X + - B 
6 4 F 


Z 

7 

= Z 

6 


The next calculations in the nodal breakdown concern areas and orientation of the curved 
portion of corrugation. This area is separated into seven segments of equal length and 

area. Hence 



A 

i 

= 1/7 (R 0 t ) 

C 1 

i = 8, 9, * • *, 14 


where the parameter 

2 R sin 0 

tc l ~ 2 R0 c 


considers the fact that the corrugation is formed by stretching a flat plate of width 2 R 
sin 0 to a circular arc of length 2 R 0, a constant volume process. The angle theta is 
now broken into seven equal segments and centers of the segment are located by the 
angular increment 0/14. Hence 

a = e/14 

a 

and 

X = R sin a 

8 

Z = R cos a - (R - H + 1 ) 

o ' o 


as seen from, the geometry of Figure 2-19. 



The last values are incremented 


a. = a. , - 2 A a 
1 i=l 


X = Esina. i = 9, 10, * * * 14 
i i 


Z = Rcosa. - (R - W - t ) 
i i s 


The other configurations of Figure 2-19 are broken up into finite elements by similar 
computational procedures. 


Stresses on the section are given by two loadings: (1) the applied bending moment due to 
normal pressure, and (2) a self-equilibrating system of internal loads due to temperature 
gradients across the section. The stress of any of the discrete elements of Figure 2-21 
is given by 


f 

i 


M (Z. - Z) 
max l + 

l 

xx 



where f-p. is the thermal stress. The critical buckling stress (an allowable stress) for 
the flat plate section of the corrugation (discrete elements 1 through 5 of Figure 2-21) 
is given by 

f cr - K ” E (f) 


where b is the width of the flat section (i. e. , b = Ag - B F ). The constant K is the 
buckling stress coefficient which depends on the edge conditions. For simply -supported 
edges with no restraints against lateral expansion, a value of K = 3. 62 is used. The 
term rj is the plasticity correction factor (see the accompanying sketch) which is given by 


r\ 



JL 

2 




The tangent and secant moduli, and E g 
respectively, are given by 


E 
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STRAIN, 




and 



where E, n, and F 7 are material properties dependent upon temperature. The expres- 
sions for the critical stress and the plasticity correction factor T) are solved 
iteratively in terms of the shear and tangent moduli. If the critical stress exceeds an 
upper limit of either 100 percent of the critical yield stress or the ultimate tensile 
stress, whichever is greater, the flat skin of the corrugation has buckled. This portion 
of the panel cross section is considered ineffective in bending, and the stress analysis 
is performed as if these elements of the cross section were not involved. Note that the 
buckling of the panel skin is applicable only to panel configurations 1 and 2, the corru- 
gation stiffened panels. 

2.4.2 TEMPERATURE INTERPOLATION. Temperatures through the panel cross- 
section members are calculated throughout the trajectory in the thermodynamic analy- 
sis. External heating rates are either computed or input, and the temperature response 
of finite elemental volumes are predicted by numerical solution of the energy equation. 

In the stress analysis, the temperatures of each of the nodes set up for the discrete ele- 
ment analysis for stress determination are predicted by linear interpolation of the ther- 
modynamic and stress models are shown in Figure 2-19 where the stress nodal indices 
are superposed on the two-dimensional conduction matrix. 

2.4.3 COMPUTATION OF STRESSES. (Thermoelastic Analysis of Statically Determ- 
inant Beams.) The stress analysis procedure is based on the finite sum method of 
Reference 20. This sub-section presents the method for determining the deformation 
and stresses of an unrestrained beam subjected to temperature variations through the 
beam cross section. The section properties for the general case of variable modulus 
(due to temperature and construction) are given in integral form for the purpose of 
determining the response of the cross section to both temperature and load. The general 
solution is then presented in integral form and is evaluated by the methods of finite sum. 

2 . 4. 3 . 1 General Solution. The following paragraphs present the general thermo - 
elastic solution for an unrestrained beam in integral form, as derived in Reference 20. 
Evaluation of these integrals, which causes most of the difficulty in obtaining numerical 
solutions to specific problems, is discussed in a subsequent section. 

The following assumptions and limitations apply to the thermoelastic solution of the 
unrestrained beam: 

a. Plane cross sections before bending remain plane after bending. 

b. Material is linear elastic at any temperature. Thus, a single relationship of stress 
to strain (a = E c) can be utilized to connect the equations of deformation and 
equilibrium and the principle of superposition can be employed. 
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c. The variation of the cross section and temperature along the length of the beam 
is both continuous and smooth and does not produce any significant shear forces. 

The unrestrained beam with temperature variation in the cross section is analyzed by 
subjecting the beam to a set of force systems which satisfy equilibrium and produce 
deformations which are compatible with the requirement of plane cross sections re- 
maining plane after bending. 

Consider a unit length of beam. Initially, each fiber of the beam is liberated from the 
influence of its neighbors. The temperature distribution is then applied to the beam 
which causes each fiber to expand by an amount oT. In general, the thermal expan- 
sion of the fibers will cause the cross-sectional plane formed by the ends of the fibers 
to warp. To satisfy the requirement of plane cross sections remaining plane, a 
pressure loading of - EaT is applied to eliminate the thermal expansion and return 
the cross section to its original position and condition (plane). This pressure loading 
upsets the equilibrium of the cross section. An axial load (F* = ftTda), equal in 
magnitude but opposite in direction to the force on the cross section due to the pressure, 
is applied at the elastic centroid of the cross section. This balancing axial force causes 
pure translation without rotation of the cross section plane so that the requirement of 
plane cross sections remaining plane is not violated. Rotational equilibrium still re- 
mains to be satisfied since the location of the resultant restoring force y, in general, 
will not coincide with the centroid of elastic area y. Equilibrium is achieved by applying 
a balancing moment to the cross section of sufficient magnitude (M' = [y - y ] F') to 
cause pure rotation of the cross-sectional plane. 

The superposed force systems now satisfy equilibrium and produce deformations which 
are compatible with the requirement of plane cross sections remaining plane. Thus, 
the procedure outlined above must result in a stress -deformation distribution for an 
unrestrained beam which is consistent and unique under the assumptions. 


Elastic Section Properties of Cross Section. The structural response to both mechan- 
ical loads and thermal stimulae is governed by the effective bending (El) and axial (EA) 
stiffness of the cross section. The equations for these properties, which are stated 
below, are identical to those for cross sections having constant E except that E is re- 
tained within the integral sign since it is permitted to vary over the cross section. 


EA = J EdA 

y = jEydA/J EdA 

z = j’EzdA/J EdA 


(2-70) 

(2-71) 
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= J E z^dA - z 2 J EdA 
= J E y 2 dA - y 2 J EdA 
= J E yzdA - y z J EdA 


(2-72) 


EI— + EI— 

(EI , EI ) = -J *— ” ± 

uu w 2 


'EI EI— \ 2 

yy zz\ 


— 2 
+ (EI— ) 

yz 


(2-73) 


The distances y and z to the elastic centroid are given by Equation 2-71. Thus, the 
elastic centroid, Equation 2-70, is the centroid of the effective elastic area EA, not 
of the geometric area. Similarly, the geometric moments of inertia are of no signif- 
icance when E varies over the cross section. The effective bending stiffnesses. 
Equation 2-72, must be employed. 


Equation 2-73 expresses the bending stiffnesses about the elastic centroid principal 
axes in terms of the bending stiffnesses about arbitrary centroidal axes. The elastic 
principal axes are defined as those orthogonal axes for which 


EI = P EuvdA = 0 
uv J 


2. 4. 3. 2 Stresses and Deformations of the Cross Section Due to Temperature. No 
thermal stresses are caused in an unrestrained isotropic, homogeneous, linearly 
elastic body by an a T distribution which varies linearly in a rectangular coordinate 
system. Clearly, in the particular case of an unrestrained beam, a linear aT distri- 
bution over the cross section produces free thermal expansions which cause cross 
sections to remain plane after deformation. No stresses are required to maintain the 
plane cross section. 


In general, a non-linear temperature distribution over the cross section would cause 
free thermal expansions of the beam fibers which would warp a plane cross section out 
of plane. Thermal stresses are produced which restore the plane cross section. In 
this respect, a temperature distribution can be considered as a 'thermal load" which, 
when applied in the absence of mechanical load, still produces stresses and deforma- 
tions. 


Following are the equations for the stresses and deformations in an unrestrained beam 
due to thermal load. 


F' 

EA 
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M' — = (z - z ) F' 

yy 

2.4.3. 3 Stresses and Deformations of the Cross Section Due to Combined Mechanical 
and Thermal Loading. For a linearly elastic beam, mechanical loads can be com- 
bined with thermal loads simply by superposition. Thus 
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cr = E [- 01 T + € t + (w z ) t (y - y ) + (w ) t (z - z) ] 

Moments about yy and zz centroidal axes are positive when their sense is such that 
they tend to cause compressive stresses in the positive quadrant (quadrant where 
both y-y and z-z have positive values). "F" is positive when tensile, and "T" 
is positive when above a datum value . 


2-40 



2 . 4. 3 . 4 Evaluation of Integrals ■ The solution of the thermoelastic beam problem 
requires the evaluation of the integrals discussed in the previous sections. The cross 
section is broken up into a finite number of elemental areas, selected so that the var- 
iation of o;T and E in each element is small. The procedure is adaptable to all cross 
sections and the degree of accuracy increases with the number of elements. 

The finite sum solution for the deformations of the cross section is based on an approx- 
imating geometry consisting of a finite number of points of concentrated elastic area 
located at the centroids of elements. Once the deformations (e, Wy, w z ) have been 
calculated from the tabular solution, stresses can be obtained at any points on the 
cross section. 

2.4.4 CALCULATION OF DESIGN FACTORS. Applied stresses vary throughout the 
trajectory depending on the gradients across the panel induced by aerodynamic pressure 
and heating, whereas the allowable stresses and elastic module vary with temperature. 
Each mode of failure of the static strength analysis is considered by determining the 
stress ratios (the ratio of applied to allowable stress) for each element in the cross 
section at each point in the trajectory at which the stress analysis is performed. Pres- 
ently, the stress analysis is undertaken at each print-out interval of the computer 
program. The most critical time for each mode of failure is determined as the maxi- 
mum value of the particular stress ratio. The computed ratios and the corresponding 
failure models are given in the following table. 



Failure Mode 

Stress Ratio 

Fortran Name 

1. 

Ultimate Tension 

f/F m 

ri 

2. 

Ultimate Compression 

f/F 

CY 

R2C 

3. 

Yielding 

f/F CY 

R2 

4. 

Crippling 

f/(F cy D) 1/2 

R3 

5. 

Elastic Stability 

f/E 

R4 

where f is the local stress of the incremental structural model node, and Frjvy> Fq Y » 
and E are ultimate tensile stress, ultimate compressive stress in yield, and modulus 


of elasticity. 

These maximum values, and the times at which they occur, determine the design point 
for the particular configuration. All elements comprising the cross section are consid- 
ered for failure modes 1, 2 and 3. The expressions for the corresponding margins of 
safely are given as follows : 
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1 . 


Ultimate Tension: 


M.S. = < f tu ./<UF) * f.) - 1 

2. Ultimate Compression: 

M.S. = (1.1 F /(UF) • f.) - 1 

U 1 

If 1. 1 F > F then the above expression is replaced by: 

Lx XU 

M.S. = (F TU /(UF)f.) - 1 
i 

3. Yielding (tension and/or compression): 

M.S. = (F TY /f.) - 1 


(2-74) 


(2-75) 


(2-76) 


(2-77) 


Note that 


a. The compression yield stress F ) is assumed equal to the tension yield stress 
( f TY) since for most materials of interest, F is not available. 

VX 

b. (UF) = ultimate factor, a factor of safety applied to the limit loads to determine 
the ultimate loads. 

The remaining failure modes (4 and 5) are applicable only to selected elements and 
configurations. For example, consider the analysis of configuration 1, Figure 2-20. 


2.4.4. 1 Crippling of Flange (Width - BF). The analysis is made for the trajectory 
point which gives a maximum value of the stress ratio 



i = 6, 7 


(2-78) 


The crippling stress is given by 


1/2 808 

F' cc = 1.385 f cy /[(F cy /E) * BF/t]* (2-79) 

When this value of the allowable crippling stress exceeds 1. 1 times the crippling yield 
stress, then the latter value is used as the allowable crippling stress. When, in turn, 
1. 1 Fqy exceeds the ultimate allowable tensile stress, then the ultimate allowable 
tensile stress is used for the allowable crippling stress. Hence, mathematically 
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(2-80) 


F 

F» 

F' s 1.1 F 

cc 

cc 

cc CY 

F 

1.1 F 

1.1 F__ < F' 

cc 

CY 

CY cc 

F 

F 

F mTT s 1.1 F_„ 

cc 

TU 

TU CY 


The panel thickness t in Equation 2-79 is given by 

t = t + t 
s c 

and the margin of safety in crippling is 


M.S. 


F 

cc 

f (UF) 
i 


2.4.4. 2 Buckling of the Semicircular Arc Corrugation (Elastic Stability) . This 
analysis is performed for the trajectory point which gives a maximum value of the 
stress ratio fj/Ej for nodes i = 8, 9 • • • , 14 of configuration No. 1. The buckling 
stress is given by 



kE t 
s 

R 


(2-81) 


where 


E is material secant modulus 
s 

k is buckling stress coefficient 
t is the corrugation thickness. 


This value of the allowable crippling stress is determined by the same procedure as 

in the previous section. When the value of the allowable crippling stress given by 

Equation 2-81 exceeds 110 percent of the crippling yield stress, then the latter value 

is used as the allowable crippling stress. When, in turn, 1. 1 F .y exceeds the ultimate 

allowable tensile stress, then this value is used for F . Hence, 

Civ 


F 

cc 


F' for F’ £ 1. 1 F„„, 
cc cc CY 


F 

cc 


1. 1 F for 1. 1 F < F' 
CY CY cc 
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F 

cc 


= F m for F ^ 1.1 F ^ F f 
TU TU CY cc 


The margin of safety (MS) is given by 


MS 


CR 

f. (UF) 


- 1 


In the computer program, if all the elements in a flange or corrugation only experience 
tensile loading, the compression crippling or buckling analyses are not applicable; in 
this event the corresponding MS is equated to TOO. If the analysis is not applicable to 
a given configuration the MS is equated to 1000. 


2.4.5 CALCULATION OF CREEP STRAIN. The total creep strain accumulated in 
each elemental area of the panel cross section is evaluated throughout the trajectory 
under the following simplifying assumptions: 


a. The creep strain accumulated in each flight is identical. 

b. The stress distributions are identical in every flight (i.e. , stress redistribution 
due to accumulated plastic strains is not considered). 

During the stress analysis, instantaneous creep rates for each element are determined 
from relevant Larson-Miller data for the particular material under consideration. The 
computational procedure is outlined below. Strain data is presented as a function of 
stress and the Larson-Miller parameter (LMP) where 

LMP = (460 + T) (20 + log t) 

and 

T = temperature in °F 
t = time in seconds 


The problem is to find the creep rate, € = de/dt, at the stress and temperature levels 
being experienced by the particular element of interest. The derivative is approximated 
by 


dc _ A€_ 
dt ~ At 

where Ac denotes the strain difference between the two curves of Figure 2-22, co- 
ordinates of which are input to the program. The difference in time, At, is determined 
from the two values of time corresponding to the two values of the Larson-Miller para- 
meter at the specified temperature and stress level. 
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Hence, since 



t 


then 
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LMP 
460 +T 



( 


LMP 
460 +T 



Figure 2-22. 


where the subscripts 1 and 2 correspond 
to the two strain curves, 1 and 2, shown 
in Figure 2-22. Upon completion of the 
trajectory, the total creep is derived by 
integrating the creep rates over the trajectory by the trapezoidal rule. Then 


Larson-Miller Plot of 
Creep Data 


e 


DST DST 



€ = NF r e. 

i 1 


where NF is number of flights. 


2.5 FATIGUE ANALYSIS 


2. 5. 1 PREDICTION OF FATIGUE LIFE. Fatigue, which results from an accumula- 
tion of alternating stress levels greater than the endurance limit of a material or 
structure, can always be induced if the applied loads are sufficiently large. (Endurance 
limit is defined as the highest stress under which repeated application can be endured 
indefinitely.) In acoustic fatigue, however, the applied loads are small. Therefore, 
the high stresses required to produce fatigue must be generated by other means. The 
prime contributor to the problem is the phenomenon of resonance. It is axiomatic 
that nearly all aircraft structural fatigue failures resulting from steady -state inputs 
such as noise or vibration are due to resonance. A structure will resonate when there 
is coincidence between one or more of its natural frequencies and frequencies of applied 
alternating forces. From the fatigue standpoint, the significant attribute of structural 
resonance is that of dynamic stress amplitude magnification, Q. 

For a linear oscillator with small damping being driven sinusoidally at resonance, 

Q = 1/2 -£L, where is the ratio of equivalent viscous damping to critical damping. 

This shows that the dynamic magnification is directly related to the damping. 
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When the oscillator (structure) is driven by a broad band random forcing function, only 
the energy near the resonant frequency of the structure performs significant work. 

This effective bandwidth, which is defined by the ’’half -power" points, is approximately 
equal to twice the damping coefficient, 2 5 ^. For example, if the resonant frequency 
(f n ) is 100 Hz and the damping is 2 percent of critical, the effective bandwidth is 4 Hz, 
i.e. , 98 to 102 Hz. Thus, increasing the damping of a structure decreases its dynamic 
stress amplification, but increases the effective driving force level because of the 
increase in effective bandwidth. If the damping of the above structure were increased 
to 4 percent of critical, the Q would be reduced by a factor of 2 while the effective 
force would increase by the factor ,/2 (assuming a flat input spectral density). Hence, 
increasing the damping of a structure is advantageous, as shown in Figure 2-23. 

Another factor which must not be neglected 
in any fatigue problem is that of stress 
raisers. Experience at Convair Aero- 
space has shown, and has been confirmed 
by independent tests under NASA contract, 
that stress raisers in conventional air- 
craft structure correspond closely to 
Krj. = 4.0 for acoustic fatigue problems. 
Experience at Convair Aerospace has also 
shown that for conventional aircraft 
structures, reduction of a stress raiser 
in itself (by changing methods of fastening, 
for example) may provide the difference 
between a structure with an adequate or 
inadequate fatigue life. This arises from the 
the fact that for a structure adequately 
designed for other service loads, acoustic 
fatigue failures almost always are initiated 
at fastener lines (because of the stress 
raisers) and then propagate into the skin 
and secondary structure (see Figure 2-24). 

In estimating the acoustic fatigue life of a TPS panel, there is a logical, straight- 
forward procedure which is carried out. Up to a point this procedure is independent 
of any cumulative damage estimation theory. 

The first step is to calculate the flexural natural frequencies of the panel. A rapid, 
approximate method has been found with errors in the lower modes of less then 10 per- 
cent. This method assumes that, in flexural vibration with small amplitudes, the 
natural frequencies are primarily a function of the flexural bending wavelengths in the 
panel. In the case of an isotropic, rectangular, simply-supported panel, the general 
expression for the frequencies of the natural modes of vibration is: 



Figure 2-23. Relationship Between Q 
and Effective Bandwidth 
(f /Q) 
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Figure 2-24. Acoustic Fatigue Failure 
Showing Crack Propaga- 
tion Along Rivet Line 



Figure 2-25. Mode Shapes of First 
Five Modes of an 
Integrally Stiffened Panel 



where 

D is the flexural stiffness per unit area 
of the panel 

M is the mass per unit area of the 
panel 

k , k are panel wave numbers in 
orthogonal coordinate directions 

For a clamped panel, wavelength equivalenc- 
ing factors can be used with the above ex- 
pression . For two -dimensionally -stiffened 
panels, the basic flat panel stiffness is used 
and to it is added the averaged stiffness of 
the reinforcing beams or corrugations, as 
may be applicable. In general, this pro- 
cedure can be employed for any structure 
that can be said to have a well defined mode 
shape which can be expressed in terms of 
simple mode shape functions (Figure 2-25). 
The number of frequencies and modes which 
must be calculated is directly related to the 
type of acoustic excitation involved. For 
the analysis of this program, only the fun- 
damental mode of each panel is considered 
as the panel's resonant frequency. 

It is next necessary to estimate the energy 
available to drive the panel at its resonant 
frequency(s). This is a function of the spec- 
tral density of the acoustic pressures and 
the damping of the panel (Figure 2-23). 
Assuming that the effective rms acoustic 
pressure is a static pressure, the equivalent 
rms static stress of the panel is calculated 
by standard procedures. This equivalent 
static stress is then multiplied by the Q of 
the panel and again by the appropriate local 
stress raiser. This, then, yields the 
dynamic rms stress at the resonant 
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frequency. Test data has shown that the statistical distribution of the instantaneous 
acoustic pressures due to rocket noise, for example, is approximately normal. Thus, 
for rocket noise, the peak pressure distribution which is of significance in the fatigue 
problem and which is the same as the peak stress distribution for a linear structural 
system can be represented by an integrated Rayleigh distribution: 



where [p > s p ] is the probability of exceeding s p in percent, and Sp/s Q is the ratio of 
peak stress/rms stress. Knowing the temporal life requirement for the panel and the 
stress distributions at the important resonant frequencies, the fatigue life can now be 
expressed in terms of numbers of significant stress reversals. 


At this point, assuming that the required S-N data is available, it is ready to be 
combined with a cumulative damage estimation procedure, as shown in Figure 2-26. 



Figure 2-26. Acoustic Fatigue Estimation Procedure 

2.5.2 LIFE PREDICTION METHODS. All theories relating to fatigue damage under 
random loading have in them, either explicitly or implicitly, a concept of cumulative 
damage concept or cumulative partial damage. The simplest and hence the most at- 
tractive theory is the linear cumulative damage concept, or Palmgren- Miner theory. 

This states that the sum of the partial damage is equal to unity, £ — = l, where n. is 
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ALTERNATING STRESS (KPSI) 


the actual number of stress reversals at stress level Sj, and Nj is the number of stress 
reversals required to cause fatigue failure at this level. This is shown in Figure 2-27. 

The overall trend of test results indicates 
(particularly for reversed bending) that 
the linear rule generally overestimates 
fatigue lives according to Freudenthal 
(References 21 and 22). The actual values 
in any given circumstances will depend 
on the maximum and minimum stress 
levels involved in the random loading 
distribution. 



Figure 2-27. 


Linear Damage Rule 
Diagram 


According to Freudenthal, the unconserva- 
tism of the linear rule results from neglect 
of consideration of the interaction between 
the infrequent high-stress amplitudes and 
the frequent low-stress amplitudes which 
produces a disproportionately high degree 
of damage at the low-stress amplitudes. 


To salvage the basic idea with its simplistic approach, Freudenthal has proposed a 
quasilinear rule which employs stress interaction factors as functions of the stress 
spectrum. Thus, for the linear rule 


V R ■ 1/£ «>i /V S i> 


and for the quasilinear rule 


K = vs <p. m./v > 


R 


1 1 SI 


where 


V' 

R 

P. 


l 

CO. 



fatigue life estimated on basis of linear damage rule 
fatigue life estimated on basis of quasilinear damage rule 
relative frequency ratio of cycles of stress amplitude S. in spectrum 
stress interaction factors 

"characteristic" values of extreme value distributions of fatigue 
lives N 

si 
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ALTERNATING STRESS (KPSI) 


Other workers have attempted to handle the stress interaction problem in a similar 
manner with varying degrees of success (References 23, 24, 25, 26, 27, 28, and 29). 
At Convair Aerospace, another approach which has proved successful has been taken 
with respect to acoustic fatigue life prediction. 


Since the phenomenon of fatigue deals with applied oscillating loads and structural 
vibratory responses, it may be considered as one of energy absorption and/or dissipa- 
tion. Above its endurance limit (if one exists) a material (or structure) has only a 
finite or limited capacity for absorbing energy or having work done on it before failure 
will occur. Where the applied loads (and structural responses) vary sinusoidally, the 
standard S-N curve defines the amount of work which can be accomplished without 
failure. With this basic concept in mind, one can consider its extension with the idea 
that if there is a limit to the energy which a structure can absorb without failure, then, 
within certain restrictions which will be discussed later, the limiting energy for work 
should be independent of the rate, sequence, or level of application. 

If a continuous, randomly variable schedule of oscillating peak stresses is considered 
(such as an integrated Rayleigh distribution) and plotted on S-N coordinates, then it is 
apparent that the curve envelope at any given time represents the cumulative partial 
damage, as shown in Figure 2-28. However, since total damage information is gener- 
ally available in the form of S-N curves 
for sinusoidally varying stresses, the 



N (STRESS REVERSALS) 

Figure 2-28. Equivalent Fatigue Damage 
Diagram, Random Loading 


partial damage must be interpreted in 
terms of an equivalent sinusoidal stress 
level. 

If the random stress distribution envelope 
is plotted for the required number of 
stress reversals (Figure 2-28) and a curve 
B parallel to the applicable S-N curve A 
is plotted so as to be tangent to it, it is 
seen that the stress level at the point of 
tangency will, at any given time, have 
contributed the greatest partial damage. 
This stress level is referred to as the 
"critical” stress level, S cr , and if the 
point of tangency were with the actual 


S-N curve, failure obviously would al- 
ready have occurred. However, if the critical stress level, S cr , is extended to N cr 
stress reversals so that the area S CI> N cr is equal to that under the random stress 
envelope, the work accomplished at N C r reversals by stress S cr is equal to that accom- 
plished at the required number of stress reversals by the randomly varying stresses. 

If N cr is less than the required service life, no failure would be predicted. Life expec- 
tancy would be prorated on the basis of the ratio between and the required service 
life. 
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Any simplistic view of cumulative fatigue damage ignores many parameters which 
could modify the analytical results. For example, the improvement in endurance limit 
of some materials by exposure to very low stress levels is an obvious point. Or again, 
the sequence in which a structure is exposed to a wide range of stress levels may be 
important. One point that cannot be ignored, however, is that if stress maxima ap- 
proach yield values, no simple theory will hold. Another important consideration is 
that if the range of stress levels is such as to cause redistribution of the structural 
stress pattern, as the stresses vary from minimum to maximum, again no simple 
theory will hold. 

Convair Aerospace has successfully used its equivalent work acoustic fatigue analysis 
procedure on its own 880 and 990 commercial jet aircraft, on the North American A3J 
engine inlets, and on the C-141 and C-5A empennages. The procedure is presently 
being used on NASA Contract NAS1-9793 (LRC), "Coated Columbium Alloy Heat Shields 
for Space Shuttle Application. " It is also being applied to the acoustic fatigue evaluation 
of the Convair Aerospace SSV booster. 

For the Convair Aerospace technique, the digital computer program computes sound 
pressure levels from basic engine and aerodynamic data, computes the dynamic re- 
sponse characteristics of the TPS panels as a function of their material properties, 
geometry, damping and restraint, computes equivalent static and dynamic stresses 
and computes fatigue life. All input data, except S-N data, is programmed. S-N data 
is included by a curve-fitting process for the individual materials under consideration. 
Details of the computatioial procedure are given in the following paragraphs. 

2.5.3 PROGRAM ORGANIZATION. The sonic fatigue analysis subroutines are 
organized into a number of basic functional tasks which include input, determination of 
panel fundamental frequencies, noise source computations, and the calculation of (1) 
dynamic stresses, (2) the number of stress reversals, and (3) the critical stress levels. 
Each of these computations is explained in the following subsections. 

2. 5.3. 1 Input. Input parameters are read into the computer in records which allow 
description of each of the four noise sources as well as panel geometry and the allowable 
S-N data of the material in question. The latter are coefficients to a cubic least-squares 
curve fit of stress (in thousands of pounds per square inch) as a function of the number 
of stress reversals to failure. The moments of inertia which are input to the acoustic 
fatigue analysis are for the panel cross section under investigation and lie in the plane 
of the panel. The parameter is used to compute dynamic bending stresses in the 
classical manner. At the moment, the panel moments of inertia must be computed 
external to the program since the computation itself is rather lengthy and a function of 
complex geometry, the size and number of corrugations, and the like. The moment of 
inertia for the cross section of the panel can be determined by superposition of the 
moments of inertia of the simpler areas of components comprising the total structure . 
The technique can be found in any handbook or fundamental text of stress analysis. The 
list of input parameters is given in Table 2-2 along with the appropriate units of each 
parameter. 2-51 



Table 2-2. Input Parameters for Sonic Fatigue Analysis 


DT(1) 

XL 

REY 

VU 

QL 

AMACH 

DT(2) 

TT 

WE It 

D 

VS 

XI 

DVEH 

YCL 

DREF 

IPAD 

DT(3) 

AE 

UJ 

WEJ 

VV 

TJ 

XJ 

YP 


Period of turbulent boundary layer noise excitation, sec 

Run length of turbulent boundary layer, ft 

Local Reynolds number 

Local velocity, ft/sec 

Local dynamic pressure, psf 

Local Mach number 

Period of rocket engine noise excitation 
Rocket engine thrust, lb 
Rocket engine weight flow, lb/sec 
Rocket nozzle exit diameter, ft 
Local speed of sound, ft/sec 

Distance between point of interest and rocket engine exit 
Vehicle diameter, ft 

Y -distance from vehicle centerline to point of 
interest, ft 

Distance from rocket exhaust plane to reflecting surface 
on pad, ft 

= 0, vehicle not on pad 
/ 0, vehicle on pad 

Period of jet (flyback) engine noise excitation, sec 
Nozzle exit area, ft^ 

Jet velocity, ft/sec 
Jet engine weight flow, lb/sec 
Vehicle velocity at flyback cruise, ft/sec 
Jet engine thrust, lb 

Axial distance from point of interest to jet engine exit 
nozzle, ft 

Radial distance from point of interest to jet engine exit 
nozzle (less than 200 ft but greater than the nozzle exit 
diameter), ft 
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Table 2-2. Input Parameters for Sonic Fatigue Analysis, Contd 


DT(4) 

HP AN 
HC 
AI 
AIY 

AW 

BW 

EP 

KFLEX 

NPAN 


C(l) 

C(2) 

C(3) 

RHOP 


Period of jet (flyback) engine scrubbing noise excita- 
tion, sec 

Panel thickness, in. 

Core thickness, in. (for honeycomb sandwich) 

4 

Panel moment of inertia, in. 

Panel moment of inertia, in.^ (for normal direction 
of corrugated panel) 

Panel length, ft 

Panel width, ft 

Modulus of elasticity for panel, psi 

Flexural rigidity index 

= 0, if structure symmetrical, rigid 

= 1, if structure unsymmetrical and/or flexible 

Panel configuration index 

= 1, flat plate 

= 2, honeycomb sandwich 

= 3, integrally stiffened 

= 4, corrugated 

First coefficient of a least-squares, third-order curve 
fit of allowable S-N data 

Second coefficient of a least squares, third order curve 
fit of allowable S-N data 

Third coefficient of a least-squares, third-order curve 
fit of allowable S-N data 

NOTE : The curve must be fitted to the allowable S-N 

data with the ordinate, S, in thousands of 
pounds per in. 2 and the abscissa as the log- 
arithm of the number of cycles, N. 

Panel density, lb/ft^ 
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2. 5. 3. 2 Calculation of Fundamental Frequencies. For the present analysis, it is 
assumed that the panel resonates (and suffers its fatigue damage) at its fundamental 
frequency. Four different types of rectangular panels are currently considered: 

(1) isotropic, (2) honeycomb, (3) integrally-stiffened, and (4) corrugated. The pro- 
gram itself determines which configuration to utilize. The technique for calculating 
the fundamental mode of each of these configurations is essentially the same. 

Rectangular Isotropic Panel. The fundamental frequency is computed as a function 
of the panel unit mass, its flexural rigidity, and the length of its shorter side. Hence 


f = 


K 


2 IT a 2 y/ M/D 


Hertz 


M, the mass per unit area of the panel, is given by 


M = 


ph_ 

g„ 


where 


3 

p = panel density in lb/ft 

h = panel thickness in inches 

g Q = gravitational constant = 32.2 ft/sec 


(2-82) 


(2-83) 


and the flexural rigidity D is given by 


D 


with 


E h" 


12 (1 - v ) 


E = modulus of elasticity 


(2-84) 


v = Poisson's ratio (for materials of interest, this is approximately 0.3) 

It is assumed, based on actual airplane experience , that the plate vibrates as an 
average between a clamped edge and a simply supported edge along all four edges. 

For the rectangular isotropic panel, this coefficient is given by the following table. 

In addition, if the plate supporting the structure is symmetrical and rigid compared 
to the plate, the frequency is as calculated from Equation 2-82. If, however, the sup- 
porting structure is unsymmetrical and/or flexible compared to the plate, then the fun- 
damental frequency is reduced by a factor of 4, i.e. , fj = f/4. 
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b/a K avg 


1.0 

27.89 

1.5 

20.63 

2.0 

18.45 

2.5 

17.61 

3.0 

17.08 

00 

16. 12 


Rectangular Honeycomb Panel. The fundamental fre- 
quency is computed by Equation 2-82, the same as the 
rectangular isotropic plate, except now the flexural 
rigidity is given by 


D = 


E i \ ‘ f * t 

12 12 


(Vc) 


4 (1 - 


v >( E f *1 +E f \ \ 
\ 1 1 2 2 / 


where the subscripts fj and denote the outer and inner face sheets respectively 
(Reference 31). The geometry is given in Figure 2-29. Under the assumption that the 
face sheets are the same thickness and Poisson's ratio is approximately 0.3 for 
materials of interest, this expression reduces to 



Et, t t 
f c 

1.82 
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Figure 2-29. Rectangular Honeycomb Panel Geometry 

It is assumed that the sandwich panel vibrates as simply -supported on all four edges. 
Thus, the coefficient K of equation (2-82) is given by the following table for the honey- 


comb panel. The panel frequency is then given by 


b/a 

K 

combining Equations 2-82 and 2-85. 

1.0 

19.74 

Integrally Stiffened Panel. This configuration is often 

1.5 

14.26 

referred to as a machined plank. The geometry is 

2.0 

12.34 

given by Figure 2-30. If pure flexure is considered 
and all torsion neglected, the overall flexural stiffness 

2.5 

11.45 

of the panel can be taken as the effects of the panel and 

3.0 

10.97 

stringer stiffness in parallel (Reference 32). The total 
flexural stiffness is the sum of that of a homogeneous 

oo 

9.87 

panel plus that of a number of beams averaged over the 


panel width. The assumption is valid for evenly dis- 
tributed stiffeners and is applicable only to the fundamental mode where all stringers 
can be assumed in flexure. The flexural rigidity for a uniform panel is given by 



12 (1 - v 2 ) 
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Figure 2-30. Integrally-Stiffened Panel 
Geometry 

a homogeneous panel, the flexural stiffness 



which, for metallics of interest where 
v = 0.3, simplifies to 


10.9 

The panel wave numbers in the x- and 
y-coordinate directions respectively are 
given by k m = mn /a and 1^ = nrr/b 
where m and n are the mode numbers 
or the number of half wave lengths. 

For the simply-supported panel in the 
fundamental mode, m = n = 1. For 
is given by 


and if the stiffnesses in the two coordinate directions are uniform but not the same, 
the flexural stiffness is given by 


, 4 2 2 4 

D 7 = D k + 2 D kk+Dk 
x m xy m n y n 

where D x and Dy are flexural rigidities in the coordinate directions, and Dxy is a cross 
flexural rigidity. If, in addition, stiffness is added in one direction by stiffeners or 
stringers, then the rigidity in that direction and the cross-flexural rigidity will be 
altered. Assuming that the stiffness is added in the y-direction, then the stiffness in 
the x-direction, D x , remains constant. Assuming thin stringers, the stiffness in the 
y-direction becomes the sum of flexural stiffnesses of the homogeneous panel, D, and 
of the stiffeners per unit panel width if they are spread evenly across the panel surface. 
(One is assuming that the stringers do not affect the panel mode shape. ) Thus 


D 


x 


= D + 


El 



Eh" 


12 (1 - v) 


Eh 


3 


10.9 


With the cross flexural rigidity given by 



/ 
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the total flexural stiffness of the structure is given by 


K. 


■ [< 


1/2 2 1/2 
(D ) ' k + (D ) 
x m y 




For the mass per unit area of the configuration M-p being the sum of the masses per 
unit area of the unstiffened panel and the uniformly distributed stringers, then the 
circular frequency of the integrally stiffened panel in purely flexural modes is simply 





2 

k + 
m 




( 2 - 86 ) 


and the natural frequency is given by 

f = Hertz 

2 rr 

Corrugated Panel. The computation of the natural frequency of the corrugated panel 
(Figure 2-31) is similar to that of the integrally stiffened panel. It is assumed that 

a. Flexural properties are uniform but not necessarily equal along the x and y axes. 

b. The flexural rigidities along the x and y axes are given respectively by D and D . 

x y 



X 


Figure 2-31. Corrugated Panel Geometry 

c. The panel is vibrating in its fundamental mode such that the wave numbers are 

given respectively by k = 77 /a and k = ir/b. 

m n 

1/2 

d. The cross -rigidity is given by D = (D D ) 

xy x y 
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Then the flexural stiffness is given by 





(2-87) 


where Dy = EI x /a and D x = Ely/b are the flexural stiffnesses per unit panel width. 
For a mass per unit area M-p given approximately by 


the circular frequency of the panel is then 



and the natural frequency is 


f = ~~ Hertz 

27 7 


(2-88) 


2. 5.3. 3 Calculation of Noise . The following paragraphs describe the computation 
of noise produced by four different sources: the turbulent boundary layer, the boost 
rocket engines, and the flyback (jet) engines (including noise on the vehicle and scrub- 
bing of the jet exhaust on the panel or location of interest). 

Turbulent Boundary Layer. The theory developed for these computations relates to 
the attached turbulent boundary layer of a large aircraft structure. It is based on 
wind tunnel and aircraft data (Reference 32) and has been modified to fit analytical 
expressions (Reference 33). 


The overall intensity is given by 


_p_ = 0.012 

q 1+0.14 M 2 

where p and q are the overall and dynamic pressures respectively in psf, and M is 
the local Mach number. 
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The sound pressure level is given by 


FPL 

f 


20 log 

10 


P (f) 


i41. 8 x 10 



db 


where 


P (f) 


q 


0.012 


\l + 0.14 M 




(2-89) 


(2-90) 


The parameter f is the fundamental frequency of the panel, whereas the characteristic 
frequency of the boundary layer f is given by 



with u being the local velocity in feet per second, x the local run length in feet, and the 
ratio of boundary layer thickness to local run length being given in terms of Reynolds 
number as 


6 ^ 

x 


0.37 Rey 


-1/5 


1 + 


Rey 


,2.9 x 10 


0.1 


Rocket Noise on Vehicle at Lift-off. The sound pressure level measured in db in 
octave bands over the surface of a rocket-powered vehicle on the launch pad is given by 


OBSPL = 10 log i() (0.676 mV 2 ) + 


SPL - 20 log, n R — AB' 
o 10 


(2-91) 


where 


9 2 

(1) 0.676 mV^ = 0.676 t g 0 /w is the mechanical stream power of the 
rocket exhaust in watts when t is the total thrust in pounds, w is the 
total weight flow in pounds per second, g 0 is the gravitational con- 
stant in ft/sec^, and V = tg 0 /w is the gas velocity at the nozzle exit 
in ft/sec. 

(2) SPL 0 is a reference octave band sound pressure level at a distance 
R of one foot. 
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(3) R is the distance from the point of interest on the vehicle surface to 
the noise source location in the jet stream. 

(4) AB' is a correction for ambient conditions to obtain levels actually 
existing at the vehicle surface. 


Calculation of the mechanical stream power of the rocket exhaust is straight-forward 
and simple. However, determination of the reference sound pressure level SPLq is 
a more complicated matter. This is dependent on the geometry of the booster -orbiter 
configuration and characteristics of the rocket exhaust as represented in the Strouhal 
number fD/V where f is the fundamental frequency of the panel, V is the rocket ex- 
haust velocity, and D is the nozzle diameter. The reference octave band sound pres- 
sure level is then given by 



(2-92) 


The term involving R, the average distance from the point of interest on the TPS panel 
to the noise source in the rocket exhaust stream, is best explained by referring to 
Figure 2-32. The average distance from the noise source of frequency f to the rocket 
nozzle exit plane is the distance X Q . Again, in terms of the Strouhal number fD/V, 
this is given by the expressions: 



0.222 + 1.315 log 1Q 



and 



0. 175 



= 0.625 - 0.202 log (y) 



^ 0.175 
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ORBITER 


NOTE: 

LATERAL DISTANCE 



In flight, when the exhaust stream extends directly 
back of the vehicle, the characteristic distance R 
is given by 

R = X + X 


However, on the pad, and when X Q > a (i.e. , 
when the noise source is downstream of the de- 
flector), 


R 


= [(X + a) 2 + (X q - a) 2 ] 


1/2 


The correction for ambient conditions is neglected 
since the reference conditions are assumed to be 
520°R and 14.7 psia. However, for panel funda- 
mental frequencies 


IT D 


veh 

Figure 2-32. Booster/Noise Source 

Relationship where C a is the local sound speed and D yeli is 

the vehicle diameter; a correction of 6 db is added to obtain levels actually existing 
at the vehicle surface. Thus 


A B' 


6 db 


for 


f > 


C 

a 


77 D 


veh 


Finally, the sound pressure level at a given center frequency is computed (with the 
octave band sound pressure level as a reference) as 

SPL = OBSPL - 10 log (BW) db 
f 6 10 

where BW is the bandwidth which corresponds to the particular center frequency as 
given in Table 2-3. The spectrum pressure is computed as shown above, where the 
BW corresponds to the particular octave bandwidth. 

Flyback (Jet) Engine Noise on the Vehicle. At the present time there is no general 
method for predicting the near -field noise of either subsonic or supersonic jets. 

The procedure employed here, which is empirical, is based on a valid far -field pre- 
diction (Reference 34); it is then corrected to near -field conditions based on experimental 
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Table 2-3. 

Bandwidth 

Frequencies 

f 

BW 

(Hz) 

(Hz) 

2 

1.37 

4 

2.75 

8 

5.5 

16 

11.0 

31.5 

22.5 

63 

45 

125 

90 

250 

180 

500 

355 

1000 

700 

200 

1400 

4000 

2800 

8000 

5600 


data (References 35 and 36). Values of the near- 
field sound pressure levels obtained by this procedure 
will generally apply at a distance aft of the jet nozzle 
corresponding to termination of the jet core (i.e. , 
X/Ds= 5) (Figure 2-33). Values will be typical (±5 db) 



Figure 2-33. Schematic of Jet Flow Field 


to about twice the distance of the jet core termination (i.e. , X/D as 10). At distances 
greater than X/D as 10, and at side distances of y as 2, the jet boundary will wipe the 
airframe and pressures in this region will be controlled by the actual jet pressure at 
the particular location of interest. These calculations will be outlined in the following 
section entitled "Flyback (Jet) Exhaust Scrubbing. " The procedure outlined herein is 
valid for jet velocities between 1000 and 2500 feet per second. 


The overall sound pressure level at a distance of y = 200 feet and at an angle of maxi- 
mum radiation is given by the expression 

OASPL 200 = 10 lo ^io (f (V R )) + 10log io (p2 A e ) (2_94) 

where 



A 

e 


= relative jet velocity (i.e. , jet velocity minus aircraft forward 
velocity) ft/sec 

= nozzle exit area in ft 2 


w 


p = rr = 

e J 


engine exhaust density where w is the engine weight flow in lb/sec, 
A e is the nozzle exit area in ft^, and Vj is the jet velocity in ft/sec. 
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The function 10 log f(V r> ) is given by 

XV 

10 '°ho f(V e> - 145 + 100 Iog 10 (xe^r) 

which has been normalized for unit density and unit exit area. 

The overall sound pressure level in the near field must now be computed by correcting 
the far -field value (i.e. , y = 200 ft). This is done in terms of the dimensionless para- 
meter y/D where D is the nozzle exit plane diameter. The correction takes the form 

A db = B log 1Q 

where 

B = 20 db/decade for 30.0 5 y/D 

B = 16 db/decade for 2 . 5 < y/D <30.5 

B = 14 db/decade for 1.0 s y/D < 2.5 

The terms y£ and y 2 are dummy parameters to cover the distance from y 2 = 200 feet 
to the point of interest; i.e. , y. 



A sample calculation of this correction is given below for the case of the radial distance 
from the engine centerline to the point of interest of 6 feet (i.e. , y^ = 6 ft) and a jet 
engine nozzle exit plane diameter of 4 feet (i.e. , D = 4 ft). Thus, the total correction 
is given by 


A db 



+ 16 log 

10 


30 

2.5 


+ 14 log 
6 10 



An additional correction of 6 db is made to account for the sound pressure level actually 
experienced by the structure. Hence, the near-field overall sound pressure level is 
given by 

OASPL . = OASPL + Adb + 6 db 
nf 200 


Finally, for the case of the jet flyback engine, the octave band sound pressure level is 
computed by correcting the overall sound pressure level for the bandwidth of the center 
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frequency. The octave band sound pressure level (OBSPL) will have a maximum value 
when the Strouhal number is 0. 8, i.e. , a characteristic frequency is given by the 
expression f / D/V = 0.8 where D and V are nozzle plane diameter and velocity. Next 
the center frequency of the bandwidth in which this characteristic frequency f ' falls is 
computed from Table 2-3, and the level of the maximum octave band sound pressure 
level is calculated as 


OBSPL = OASPL „ - 5 db 

max nf 

Below this maximum value, the OBSPLs decrease at the rate of 4 db/octave, and above 
it at 3. 5 db/octave. This decrease from the maximum is computed by determining how 
many octaves the fundamental panel frequency is above or below the center frequency 
of the octave band in which i‘ falls. 

Next, the spectrum pressure level SPL^ is computed as 


SPL f = OBSPL — 10 log iQ (BW) 

where the bandwidth BW is given in Table 2-3 for the band in which the panel fundamental 
mode, f, falls. Finally, the spectrum pressure is computed as 


p(f> = 10 


(SPLf 


2.09 x 10 



psf (rms) 


Flyback (Jet) Engine Exhaust Scrubbing. For the case in which the jet exhaust stream 
actually impinges on the panel of interest (i.e. , the panel falls within a seven-degree 
half angle), the component of actual jet pressure is computed as follows. First, the 
normal pressure in the exit plane is given by 


T 1.275 T 



e 


where T is engine thrust in pounds and D is the exhaust nozzle exit diameter, 
grazing jet pressure P is given by 

X 


p 

x 


0.155 T 



x 


The 


where D is the jet diameter at a station x feet from the exhaust plane, i.e. , 

X 

D = D (l + 0.244 — \ 
x e \ D e / 
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The overall sound pressure level can now be computed from the grazing pressure as 


OASPL 

x 


20 log 


10 ^ 2. 09 x 10~ 8, 


db 


Next, as was the case for noise due to the jet engine alone, the maximum octave band 
sound pressure level is computed at the characteristic frequency f 'which occurs for 
the value of the Strouhal number, 

f' D 

- 0.4 

V. 

J 

Using the center frequency for the octave band in which f' falls (Table 2-3), the maxi- 
mum octave band sound pressure level is computed 

OBSPL = OASPL - 5 db 
max x 

and the deviation from this maximum (either 4 db/octave below or 3. 5 db/octave above) 
is computed by comparing the panel' s fundamental frequency to the center frequency 
computed above. Finally, the spectrum pressure level and the spectrum pressure are 
computed respectively from 


SPL f = OBSPL - 10 log iQ (BW) db 
and 




x 10' 8 SPL.) 


psf (rms) 


where the bandwidth BW (given by Table 2-3) psf (rms) corresponds to the panel's 
fundamental frequency. 


2 . 5 . 3 . 4 Calculation of Dynamic Stresses in a Resonating Panel. The fundamental 
frequencies for each of the panel configurations of this study have been calculated as 
outlined in the previous sections, and all noise sources have been investigated to 
determine sound pressure levels and rms acoustic pressures for each source. Next, 
the dynamic stresses for each type of panel are to be calculated for each noise source. 
It is assumed that the acoustic pressures will be phase correlated over the entire 
surface of the plate. 

Rectangular Isotropic Panel. For the case in which the panel is symmetrical and the 
supporting structure is rigid, the deflections and stresses are calculated for a clamped 
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panel. The maximum bending stresses occur at the center, and the stress and deflec- 
tion are given respectively by 


S 

max 


0.5 P a 2 


lb/in" 


rms 


Y 

max 


0.0284 P a 4 
Eh 3 


in. 


rms 


where p is the rms acoustic pressure in psi 


p = p(f) y BW 

The acoustic pressures p(f) have already been computed as shown in previous sections, 
and the bandwidth BW is given by 

BW ■ 2 (t) f 

where f is the panel fundamental frequency and the term (c c /c ) is the critical damping 
coefficient. However, the stress S^x is an equivalent static stress, and since the 
plate is resonating, it must be multiplied by the dynamic magnification factor. 



For mechanically fastened plate (i.e., rivets, screws, etc.), atypical value is 
c/c = 0.02 and Q = 25. The dynamic stress is then given by, 

S, = 25 S 

d max 


A local stress raiser of four is next included at the fastener (Reference 37). Thus 
the local dynamic stress is given by the equation 

S = 100 S lb/in 2 rms 

id max 

The stresses are considered to be distributed according to an integrated Rayleigh 
distribution with a maximum peak to rms ratio of four based on experience. Thus, 
the maximum peak stress will be given by 



max 


400 S 

max 
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For the case in which the rectangular isotropic panel is unsymmetrical and/or has a 
flexible supporting structure, it is assumed that alternate frames will twist such that 
the panel wave length is 2a instead of just a (the case for the symmetrical panel and/or 
rigid supporting structure). 

The length of the panel that will now have phase correlated pressures is 1.25a, and 
this will be treated as a clamped panel. Thus, the maximum stress and deflection 
are given by 

2 

0.5 p (1.25a) 

^max ,2 

h 

4 

Y = °- 0284 P (1.25a) 

max , 3 

E h 

The dynamic stress S d> the local dynamic stress Sj^, and the maximum peak stress 

Sr, are all computed from the maximum static stress S just as outlined in the 
Pmax max 

previous paragraphs for the symmetrical panel (Reference 30). 

Stresses in a Resonating Honeycomb Sandwich Panel. As in the case for the rectangu- 
lar isotropic panel, it is assumed that the acoustic pressures will be phase correlated 
over the entire surface of the panel. For the simply-supported panel, the maximum 
bending stresses and deflection occur at the panel center. Hence, 


- 2 
0. 75 p a 


max , 2 

h eff 


psi (rms) 


0.01421 p a 4 . 


max 


in. (rms) 


E h' 


eff 


where the effective thickness h ^ is given by 


h = 1.817 (h h h) 

eff c f 


1/3 


The dynamic stress S d , the local dynamic stress S^, and the maximum peak 

stress S n are computed as in the case of the rectangular isotropic panel. 

Pmax 

Stresses in a Resonating Integrally Stiffened Panel. For a simply -supported panel for 
which the acoustic pressures are phase correlated over the surface of the panel, the 
maximum stress and deflection are 
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psi (rms) 


max 


0 . 75 p a 2 


eff 


0.0142 p a^ 


max 


E h 


in. (rms) 


eff 


where here the effective panel thickness is given by a somewhat more elaborate calcu- 
lation than was necessary in the two previous cases. The task is to find an effective 
thickness of an isotropic panel whose stiffness equals that of an integrally stiffened one. 
For the isotropic panel, stiffness is given by 

2 2 2 

k. = D (k + k ) 
l x m n 

whereas, for the stiffened panel 


f 1/2 2 1/2 212 

k = (D ) ' k + (D ) ' k 
s V x m ' y nJ 


where 


El 

D = D + — 
y x b 


The stiffness ratio is defined as 


r s ' k 


D / , the flexural rigidity of the integrally stiffened panel, can be expressed as 

D' = D r 
x s 

where D is the flexural rigidity of the panel without stiffeners, then 
x 


D = 
x 


E h" 


E h 


12 (1 - v ) 


D' = 


eff 


2 s 
12 (1 - v ) 


which combine to give 


h = \Jv h 
eff v s 
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Stresses in a Resonating Corrugated Panel. Once again, for the simply-supported 
panel, the maximum stress and deflection are given by 


S 

max 


Y 

max 


„ „ - 2 
0.75 p a 


eff 

0.0142 p a 4 


E h 


eff 


where, as in the case of the integrally stiffened panel, the effective thickness is 
determined by equating stiffnesses of a rectangular isotropic and a corrugated panel. 
The stiffnesses for the corrugated panel and the equivalent isotropic one are given 
respectively by 


- [<D , 1/2 X 2 ♦ (D > 1/2 k 2 f 
Lx m y nj 


and 


,2 2\ 2 
k = D (k + k 
m 


if 


Equating these expressions (assuming equal stiffness) and solving for D such that 


D = 


Eh .. 
eff 

12 (1 - v ) 


gives, after suitable substitution for D , D , k , and k , 

x y m n 


eff 


10.9 


x 1/2 , T 1/2 

( j) (?)»& 6 ) 

(7-7) 


1/3 


2. 5. 3. 5 Computation of Critical Stress Levels and Number of Stress Reversals. 
In applying the Convair Aerospace particle damage theory outlined in the opening 
paragraphs of this section, the next task is to determine the critical stress levels 
and number of stress reversals corresponding to each local dynamic stress on the 
panel of interest. It is assumed that the applied stresses obey a random Rayleigh 
distribution; that is, the distribution has the form 
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2 



where [p > S p ] denotes the probability of exceeding the stress S p and (Sp/S^d) is 
ratio of peak stress to rms stress. Thus, the distribution of the number of times a 
stress reversal corresponding to a stress of a given magnitude will occur is given by 



The ratio of probable stress to the local dynamic stress as a function of the number of 
times this stress will occur over the lifetime number of stress reversals is given by the 
inverse of Equation 2-95. Such a curve is shown in dimensionless form in Figure 2-34. 


0.0001 0.001 



Figure 2-34. Rayleigh Stress Distribution 
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As a part of the present effort to predict the effects of acoustic fatigue, both the Ray- 
leigh distribution and its inverse have been curve fitted for automatic computations 
which will be explained in later paragraphs. 

The Rayleigh distribution for the applied random stress is given by the curve fits 



where 

C = 2.631102 xio” 2 C = -2.635 x 10" 2 C = -4.939331x10 

1 2 3 


Since the stresses are distributed as described, the purpose of the next calculations 
is to determine the stress which is contributing the greatest partial damage to the TPS 
panel (Figure 2-35). This clearly occurs at that critical stress (S cr ) on the Rayleigh 
distribution which is closest to the allowable S-N curve (i.e. , the locus of actual 
material failures). This is found mathematically to be the point at which the applied 
S-N curve (the Rayleigh stress distribution) is tangent to the allowable S-N curve 
(established by experiment). A critical stress is found in this manner for each of 
the local dynamic stresses (corresponding to each noise source), and a composite 
critical stress is next determined as 
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1/2 



S 


n 


£ S er. 
j=! J 


where n is the total number of local 
dynamic stresses being considered for 
the particular case at hand. Finally, 
the equivalent number of stress rever- 
sals for the sinesoidally applied effective 
stress (N cr ) is determined by equating 
the work applied by the critical stress 
S to the work done by the random stres- 
ses (i.e. , the area under the stress 
envelope) . Hence , 


Figure 2-35. Equivalent Fatigue Damage 
Diagram for Random Loading 


N = 


1.261 

§ 


n 


£ 




where the constant 1.261 has been evaluated by numerically integrating the Rayleigh 
stress curve. 


The point (S, N) is compared to the allowable S-N curve. If the point falls above the 
allowable S-N curve, then failure has occurred - the number of stress reversals at 
the composite critical stress level S has exceeded those at which the material failed 
during testing. If the point falls below the allowable S-N curve, failure due to sonic 
fatigue is predicted not to occur. 

2 . 6 TPS SECTION REDESIGN 


The TPS redesign is performed on two separate components of the system — the 
metallic re-radiative panel and the underlying insulation. The thickness of the metallic 
cover panel is varied to provide just enough strength to transfer the mechanical loads 
of air pressure to the supporting structure and balance the thermal stresses of the 
panel itself. The thickness of the underlying insulation is varied to allow enough re- 
sistance to heat transfer to maintain a specified allowable temperature in the under- 
lying structure. Only the thicknesses of the cover panel and the insulation are variable; 
all other dimensions and properties are constant. 

2.6.1 STRESS REDESIGN OF THE PANEL. A particular configuration is input to 
the computer program by specifying the configuration number (e.g., circular corruga- 
tion with skin, flat corrugation with skin, rib-stiffened panel, skin stringer, open 
corrugation, or circular corrugation) and other geometric quantities such as panel 
length, width, number and depth corrugations, radius of corrugations, and the like. 
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Throughout the trajectory, indices measuring the margins of safety for various modes 
of panel failure, as described in Section 2.3, are stored for evaluation at the end of 
the trajectory. At that time, the redesign procedure of the stress analysis is activated 
by the program. The margins are scanned and the one with the largest negative value 
identified. This, then, becomes the design point. The thickness of the panel is in- 
creased by ten thousandths (0.010) inch, and the stress analysis using the loads and 
temperatures stored for each node of the stress model for that time in the trajectory 
is performed again. It has been ascertained during program development that tem- 
perature gradients (the physical parameters which govern structural design of the 
panel) do not change significantly (many times not at all) with panel thickness. The 
margin or design factor from the new stress analysis is then investigated to see if 
the increased thickness is sufficient to handle the load. If so, the design is complete 
and program control passes to the weights/cost analysis. If not, a new panel thickness 
is chosen by extrapolation of a Newton-Raphson-type iteration. After each stress 
analysis, the margins are rechecked until the design point shows a positive margin of 
safety. A design is considered optimum when the minimum design factor is between 
0 and 10 percent. Results should be checked carefully since the panel thickness is not 
decreased during the redesign procedure; an over-design panel may result which can 
be eliminated by a subsequent computer run. 

2.6.2 THERMODYNAMIC REDESIGN OF PANEL. Heat transfer to the underlying 
structure can be reduced in two ways for a given panel configuration: the metallic 
panel thickness may be increased or the insulation thickness may be increased, both 
cases thereby increasing the heat capacity of the respective material. Of the two 
techniques, the latter is the more efficient in terms of system weight. For this 
reason, only the insulation thickness can be varied to affect the temperature of the 
structure . 

The computational procedure is as follows: For each material of the TPS, an associa- 
ted allowable maximum temperature is input. For example, an aluminum node will 
not be allowed to get above 300° F. One particular material is identified as being the 
insulation. At each new temperature calculation, material temperatures are compared 
to allowables. If one is exceeded, the insulation thickness is increased, the program 
is returned to the starting point of the trajectory, and calculations are resumed. When 
the temperature allowables no longer exceed the actuals, the design is considered 
finished as far as thermodynamic considerations are concerned. At the present time, 
the insulation variation is performed by doubling the insulation thickness, but this 
technique is currently under refinement. A typical sketch of the thermodynamic re- 
sizing procedure is shown in Figure 2-36. 

When an actual temperature exceeds the maximum allowable, the insulation thickness 
is doubled, the program returns to the starting point of the calculations, and the com- 
putations are begun anew. The insulation thickness is continually doubled until either 
the actual temperature of the critical material is less than the maximum allowable, or 
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TIME 


Figure 2-36. Thermodynamic Resizing 
Procedure 

2 . 7 WEIGHTS/COST ANALYSIS 


until the insulation thickness reaches one 
foot. For the latter case, the run is 
terminated. Here, the program user 
should again be critical of the printed 
results. This first generation computer 
program does not decrease the insulation 
thickness once a value is found which 
satisfies the temperature constraints. 
Hence, the insulation thickness may be 
over design, and only perusal of the re- 
sulting temperature distribution will 
determine this fact. A subsequent com- 
puter run with a different initial insulation 
thickness may be necessary. Messages 
describing all changes in insulation thick- 
nesses and/or termination of the case 
(along with appropriate temperatures) are 
printed out. 


2.7.1 PARTS LISTING. The parts listing procedure requires that a library be main- 
tained to call out an associated list of detailed parts when a given panel or structure 
configuration is specified. Each part has associated with it a list of manufacturing 
processes and amounts of material required to produce the part. In this way, a method 
of deriving the manufacturing costs and material costs of the complete TPS system has 
been developed with its basis at the detail parts level. The parts listing process occurs 
as a function of the panel concept (specified by KINDP = 1, 2, 3 as shown in Figure 2-37). 
A detail parts listing for the specified configuration is included in the manufacturing 
cost summary. Figure 2-38. 



Figure 2-37. Panel Concepts 
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Figure 2-38. Manufacturing Cost Summary 


2.7.2 WEIGHTS. Actual and purchase weights are computed for each detail specified 
in the parts listing process; the results are listed in the output under the headings 
ACTUAL WEIGHT and MATL WEIGHT. Actual weight is just what the name implies. 

It is computed based on the actual geometry of the finished detail part, taking into con- 
sideration all the necessary manufacturing and design requirements that normally go 
into producing a real part. Purchase weight is the weight of the raw material that 
must be purchased to produce each detail part. It is always larger than the actual 
weight. Calculation of the actual and purchase weights for a TPS detail part of the 
thermal protection system uses an equation of the form demonstrated below. 

Actual weight = density * length * width * thickness 

Purchase weight = density * (length + C ) * (width + C ) * (thickness + C ) 

1 a O 

where C^, C 2 , and Cg are incremental distances added to the part geometry to account 
for the material removed during the cutting to size manufacturing operations. The 
actual equations can be found in the subroutine WTTPS in the program listing. 

2.7.3 MANUFACTURING PROCESSES. To predict the manufacturing processes re- 
quired for each detail part, a library of shop orders and shop planning records was 
established. These processes, along with the San Diego operation standard hour data. 
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were studied and used to identify basic standard shop operations and to correlate each 
detail part in the parts list with a set of these standard shop operations. At present, 
this information is stored in the subroutine BLOCK DATA in tables. The data are used 
to compute the standard hours necessary for each shop process required in the produc- 
tion sequence of each detail part. Standard horns and their actual method of computa- 
tion are explained in later sections of this report. Seven standard shop operations 
were identified for this study: 

Manufacturing Operation Index Manufacturing Operation 

1 Cutoff 

2 Milling 

3 Forming 

4 Drilling, Routing, Deburring 

5 Surface Treatments 

6 Heat Treatments 

7 Painting, Identification 

Cutoff includes issue of the stock material to the shop and all cutting or sawing opera- 
tions required during production. Milling includes milling, boring, and turning opera- 
tions; forming includes all forming, stretch forming, and bending operations. Drilling, 
routing, and deburring operations are considered as a related group. Surface treatments 
include chemical milling, etching, anodizing, and peening. Heat treatments, painting, 
and part identification are also considered. An example of a callout for the manufacturing 
operations required to produce a corner post is: 

Corner Post (Titanium) Manufacturing Operations 
KRandKS = 1.2. 4. 6. 7 

Thus, the operations include cutoff, milling, drilling, heat treating, and identification. 

The set of operation indexes KR and KS are then used by the program along with a 
material form index KK to call out values from the arrays KSETUP (KK, KS) and KRUN 
(KK, KR) located in the K-TABLES with BLOCK DATA. These values (in hours per 
pound actual weight) are then used to estimate a setup time and a run time required to 
manufacture each detail part. It is assumed that one basic machine setup is required 
per detail part, and that the setup for each additional like part is included with the run 
time. The equations take the following form: 

Setup Time (SETUP): ACWT * KSETUP (KK, KS) 

Run Time (RUNTM) : ACWT * KRUNTIM (KK, KR) 
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where ACWT is the actual weight of the detail part. The setup and run times calculated 
are in standard hours as discussed in the following section. 

To predict the operations required to assemble the detail parts into the basic subcom- 
ponents and then to develop a complete final assembly, a library of shop assembly 
planning records was established. Basic processes were identified and correlated 
with each subassembly and with the final assembly. This information is stored BLOCK 
DATA as constants which are called out for use with equations to compute the standard 
hours for subassembly and final assembly. 


Seven standard assembly operations were identified for this study. They are given in 
the table below along with their appropriate assembly operation. 


Assembly Operation Index 
1 
2 

3 

4 

5 

6 
7 


Assembly Operation 
Setup 
Clamping 
Drilling 
Securing 
Inspection 
Disassembly 
Cleaning 


Setup includes the mounting of parts to assembly fixtures, and clamping involves the 
clamping together of the parts to be assembled. Drilling considers the location, 
center punching, and the drilling of all required holes. Securing consists of the 
actual attachment of required fasteners, and inspection is self explanatory. Cleaning 
involves the final cleaning of the finished assembly and any cleanup required during 
the assembly process, such as deburring of the holes after drilling. In some cases 
after drilling, it is necessary to disassemble the assembled parts for inspection and 
cleanup purposes. The disassembly process is essentially the reverse of the setup 
and clamping processes, and includes reassembly before the final securing is done. 
The program calculates a value in hours for each required subassembly and for the 
final assembly of the finished part. Currently, these calculations are based on con- 
stants stored with the BLOCK DATA. The actual equations take the following form. 


Setup Time: 
Clamping Time : 
Drilling Time: 


HOLES * TIME6 
CLAMP * TIME1 
HOLES * VALUE 
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Securing Time: 

HOLES 

* VALUE 

Inspection Time: 

HOLES 

* TIME2 

Disassembly Time : 

CLAMP 

* TIME 5 

Cleaning Time : 

HOLES 

* TIME 3 


Where: HOLES equals the number of fasteners (fastener holes) 

CLAMP equals the number of clamps 

VALUE equals the volume of material removed from all the holes 
x a material complexity factor based on the material type 
x TIME 7 


The material complexity factor is assumed to be 3.0 for titanium. The constants 
TIMEX where X = 1 through 7 are stored in the BLOCK DATA as: 


TIME1 
TIME2 
TIME3 
TIME4 
TIME5 
TIME6 
TIME 7 


0.035 hour/clamp 
0.0008 hour /hole 
0.01 hour /hole 
0.0015 hour /hole 
0.012 hour /clamp 
0.0012 hour /hole 
0.52 hour /cubic inch 


Provision has been built into the program so that the constants TIMEX can be replaced 
with equations. Values replacing the constants will be computed using these equations, 
which will be based on data collected and curve fit during the course of further study. 

An example of the type of data available is shown in Figure 2-39 where the hours re- 
quired for drilling in titanium are shown as a function of the volume of material removed. 


2.7.4 STANDARD HOURS. For each detail called out in the parts listing process, 
there is a corresponding list of required manufacturing operations in the cost sub- 
routine. For each manufacturing operation a calculation of required standard hours 
is made. Standard hours are defined as a standard time, measured in hours, repre- 
senting an optimum required to perform a task. Standard hours for each production 
process are established by the industrial engineering department by analysis of time 
and motion studies of typical tasks. They are used as a means of measuring perfor- 
mance by determination of realization factors (or efficiencies) when compared with 
actual labor hours acquired through labor accounting processes. 
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Figure 2-39. Drilling Time in Titanium 


Within the program, standard hours are 
calculated for two separate cost areas : 
factory production and final assembly. 
Factory production includes the manu- 
facturing and assembly of detailed parts 
into the basic subcomponents such as 
posts and beams. Final assembly involves 
the final assembly of the subcomponents 
into a finished TPS system. A listing of 
standard hours is output as part of the 
cost data under the heading STD HOURS 
(see Figure 2-38). 

Standard -hour data have been collected 
for each of the seven manufacturing pro- 
cesses and assembly operations discussed; 
they have been adapted for use with this 
program. Standard-hour values are used 
to estimate the actual labor hours re- 


quired for each of the production processes, and then to estimate the actual labor costs. 
Figure 2-40 presents an example of a (San Diego operation) shop planning order. Listed 
are the various manufacturing processes required to produce a brace and the correspon- 
ding number of standard hours required for setup and running each shop process. The 
object of the standard-hour calculation technique in the program is, in effect, to be 
able to predict the planning order . 


Calculation of the standard hours for each of the manufacturing processes is performed 
in two parts: setup time and actual run time. Setup time is derived based on the size 
of the part, the complexity of the required setup, and the type of machine to be used. 
Run time is dependent on the machine operation rate, amount of material to be removed, 
depth of cut, and the surface area to be covered. The manufacturing standard hour 
equation is: 

STDHR = SETUP + RUNTM * KT * SHIPSET 

where 


SETUP is the setup time in standard hours 
RUNTM is the run time in standard hours 

KT is the number of parts required of a given detail to produce a 
complete component 

SHIPSET is the number of shipsets of that component produced 
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Figure 2-40. Example of a Shop Planning Order for a Brace 
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The setup and run time standard hours are calculated from stored data. These values 
are constants in terms of hours per pound of material, which were determined for 
each manufacturing process required to produce each detail part. Data in the tables 
at present are valid only for titanium, but provision has been made for any number of 
new materials. 

Calculation of the standard hours for each subassembly operation and for the final 
assembly is performed using the equations discussed in the previous section. The 
value of standard hours determined for subassembly processes are summed and added 
to the manufacturing process standard hours and listed for each detail part. The stan- 
dard hours determined for the final assembly of the completed component are listed 
separately after each assembly process. 

2.7. 5 REALIZATION. Realization factors are the ratio of standard hours assigned 
for a given process to the actual hours required. Realization for a process is used as 
a means of measuring actual performance against a standard time that is typically an 
ideal number of hours. It is a measure of shop efficiency and, as such, varies from 
department to department and from day to day within a department. Realization data 
for the various departments involved in production tasks at the San Diego operation 
have been collected, studied, and adapted for use with the program. Since these data 
take into account the effects of the learning curve for a specific task, realizations 
can be specified either as an average value or as a time -dependent variable. Some of 
the factors affecting realization are : 

a. Inaccurate planning of the required work, setup times, or run times. 

b. Machine breakdown. 

c. Change in machine, tools, or procedure which are not reflected by corresponding 
change in standard hour estimates. 

d. Tool breakage and part spoilage . 

e . Availability of previous setups . 

f. Use of special supervision. 

g. Ability and level of effort of individuals assigned to the task. 

Some of these factors are subject to control by managers and foremen, but they can 
also vary based on the current shop work load. Realization factors are useful in 
determining the overall effect of deviations from standard hour estimates. Realization 
for a particular task is defined as standard hours divided by actual labor hours. 

In the program, values for realization are stored for each manufacturing and assembly 
process, and since standard hours can be estimated, the actual labor hours can be 
calculated: 

Actual labor hours = standard hours/realization. 
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Realization factors, along with labor rates and overhead ratios, are stored within the 
program. Provision has been made in the program for a different realization factor 
for each of the manufacturing processes, for subassembly, and for final assembly 
operations. However, at present, a constant realization factor of 0.40 is used for 
all operations, but data is readily available to establish an individual realization for 
each operation (although these factors would be strictly valid only within the San Diego 
operation shops). 

Average realization factors, along with the computed labor hours, are listed for each 
part in the cost output data under the headings REAL FACT, LABOR HOURS, and 
TOTAL HOURS. Figure 2-41 illustrates some typical realization factors and manu- 
facturing standard hours plotted as a function of the number of ship sets. 






Figure 2-41. Typical Realization Factors and Standard Hours 

2.7.6 LABOR AND OVERHEAD RATES. Labor and overhead rates are used within 
the program to calculate appropriate costs, based on the number of actual labor hours 
required for each manufacturing and assembly process. 


Labor rates reflect the wages paid directly to the individual employees for each hour 
of clock time. The rates do not include fringe benefits or company contributions to 
retirement. Social Security, and state unemployment; these are considered part of the 
overhead cost. Also included as part of overhead are indirect labor costs, maintenance, 
supplies, taxes, insurance, and depreciation. 


2-82 



LABOR RATE ($/hr) 


Labor rates are largely uncontrollable by management, being a function instead of 
union/management agreements and reflecting current labor supply and demand, general 
economic conditions, and inflation. Labor rates are a function of time and are readily 
predictable over the short term, although the incorporation of time -dependent equations 
has been left for future work. At present, the program can store a labor rate corres- 
ponding to each manufacturing process and the assembly operations. However, a con- 
stant manufacturing labor rate of $4.75 per hour is used, and for the assembly operations 
a constant rate of $5.20 per hour is used. 

The overhead ratio is the ratio of overhead costs to labor costs and can, therefore , be 
used to determine an effective overhead rate : 

Overhead rate = Overhead Ratio * Labor Rate 

The overhead ratio is a useful tool for estimating purposes, and is readily available 
based on past labor cost and overhead cost data. The program can store an overhead 
ratio corresponding to each manufacturing process and the assembly operations. A 
constant ratio of 1.75 is used at present. Representative values of labor rates and 
overhead ratios as a function of time are presented in Figure 2-42. Values for the 
realization factor, labor rate, and overhead ratio are stored in the program. 

The average labor rate for each part 
is output with the cost data, along with 
the corresponding average overhead 
rate (overhead ratio for final assembly 
costs). Labor cost and overhead 
costs are calculated for each part as 
the product of the labor hours and 
labor rate , and the labor hours and 
overhead rate, respectively. Factory 
costs and assembly costs, then, are 
the sum of the corresponding labor 
and overhead costs. 

2.7.7 MATERIAL. Material costs are 
computed based on the raw material 
purchase weight of the given material 
type and form considered. Calculation 
of material costs required the deriva- 
tion of a unit material cost (in dollars 
per pound purchased), which is then multiplied by the purchase weight and by the man- 
ufacturing usage variance factor : 

Material Cost = COSWT * MAWT * MUV 



Figure 2-42. Typical Factory Direct Labor 
Rates and Overhead Ratios 
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Total material cost is the sum of the material costs computed for material type and 
form required to produce each detail part. Material costs are listed in the output for 
each detail part and are summed for the complete structure. The sum of material 
cost and factory cost is the fabrication cost, which is also listed in the output for each 
detail part and each component of the complete structure . The total manufacturing 
cost is the sum of the total fabrication cost and the total assembly cost. 

Material purchase weight is discussed in the section labeled WEIGHTS. The unit 
material cost is a function of several factors including: 

Type: Titanium, steel, etc. 

Form: Bar, sheet, plate, casting, extrusion, 

forging, fastener, etc. 

Size, shape, complexity: Standard stock, special lengths or widths, 

special gages or tolerances, complex 
cross sections, etc. 

Alloy, temper: Standard stock, special alloying or heat 

treating, etc. 

Availability : Standard stock or special mill run with 

corresponding setup changes and die or 
mold costs. 

Quantity, delivery schedule : Number of shipsets required minimum 

buy requirements, special large volume 
rates, need time, changes in unit cost 
with time, etc. 

Packing, shipping: Standard or special handling. 

Material cost as a function of these items has been collected and reduced to a cost per 
unit weight basis depending on material type and form, and stored in the program 
along with the corresponding manufacturing usage variance factor. Material cost data 
are stored in an array; they represent typical values for material unit cost for each of 
the 16 material forms considered to date by this study (Table 2-4). Material cost (in 
dollars per pound purchased) is called out by a material form index and a material 
type index, respectively. At present, values for material cost are stored for titanium 
only. The capability is present to add any number of additional materials. Current 
plans call for the addition of a logic sequence that will compute a material unit cost 
using variables instead of constants for the previously listed items. This capability 
will also include time. Hence, assuming the labor and overhead rates are also time 
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Table 2-4. Summary of the Available Material Forms and the 
Corresponding Material Form Index 


Material 
Form Index KK 

Material Form 

Part Reference 

1 

Sheet Tl 

Edge pieces 

2 

Sheet T2 

Skins, seal strips 

3 

Sheet T3 

Beams 

4 

Sheet T4 

Corners, doublers 

5 

Sheet T5 

Ribs, posts 

11 

Tubing D1 

Posts 

12 

Tubing D2 

Post support rings 

13 

Tubing D3 

Post support tubes 

21 

Rod R1 

Plugs 

51 

Honeycomb 

Honeycomb core 

61 

Insulators 

Insulators 

71 

Honeycomb inserts 

Honeycomb inserts 

81 

Fasteners 

Fasteners 

91 

Insulation 1 

Insulation 1 

92 

Insulation 2 

Insulation 2 

93 

Insulation 3 

Insulation 3 


dependent, the complete cost portion of the program will be able to account for changing 
costs with time . Also included will be a summation of each general type of material 
form; the corresponding cost will be computed based on a lot buy for the total material 
requirements of that form, reflecting minimum-buy penalties, quantity -buy discounts, 
special mill charges, die costs, form complexity factors, and tolerance requirements. 
This type of data has been accumulated, and it is only necessary to curve-fit them and 
write in the corresponding equations to adapt this logic to the program. Table 2-5 is 
an example of available data reflecting mill pricing policies with respect to the quantity 
of material bought for aluminum. 

The manufacturing usage variance is the ratio of the actual material purchased to the 
original estimated material required. The variance factor is established by accounting 
practice, and is the result of material and part overbuying, losses, spoilage, duplication, 
substitution, and changes. Table 2-6 shows some typical values for the manufacturing 
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Table 2-5. Quantity Buy Price Differentials for Typical 
Aluminum Extruded Items 


Quantity Ordered Mill Standard Items 

(Pounds) (Dollars added per pound) 


30,000 and over (basic price)* 

$0.0 

20,000 thru 29,999 

0.004 

10,000 thru 19,999 

0.008 

5,000 thru 

9,999 

0.030 

2,000 thru 

4, 999 

0.055 

1,000 thru 

1,999 

0.080 

500 thru 

999 

0.110 

300 thru 

499 

0.175 

100 thru 

99 

0.375 

25 thru 

49 

0.575 

less than 25 


0.727 


*The base price applied to 30,000 pounds of any single item, or to 
30,000 pounds of grouped items with a minimum of 200 pounds of 
any one. 


substitution and changes. Table 2-6 shows some typical values for the manufactur- 
ing usage variance factor for a past commercial transport aircraft program. The 
subroutine uses a constant value of 1. 10 throughout, corresponding to a 10 percent 
variance between actual and estimated material costs. Values for the manufacturing 
usage variance factor are stored in an array KMUV (KK, JJ), also located with the 
BLOCK DATA and called out by the indexes KK and JJ. The K-TABLES KCOSWT 
and KMUV are discussed under K-TABLES. Table 2-7 summarizes the values stored 
in the KCOSWT table. 

2.7.8 PROGRAM COST SUMMARY. An example output of a program cost summary 
is presented in Figure 2-43, The total program costs are derived on the basis of the 
TPS gross weight, usually in terms of a Theoretical First Unit Cost or TFU. This is 
the predicted production cost of the prototype article including manufacturing planning, 
fabrication, subassembly, sustaining engineering, sustaining tooling, quality control, 
materials, and subcontracted parts. It is calculated using an equation of the form 

TFU = C 1 = (WTTPS) a 
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Table 2-6. Typical Manufacturing Usage Variance Factors for a Past 
Commercial Transport Aircraft Program 


Contract 
Lot No. 

Actual 

Material Costs 
A 

(millions of $) 

Original Estimated 
Material Costs 
E 

(millions of $) 

Percent Variance 
(A - E/A) * 100 
(percent) 

Manufacturing 
Usage Variance 
Factor 

1 

40.65 

34.74 

17.0 

1.170 

2 

4.61 

4.25 

8.5 

1.085 

3 

16.67 

14.39 

13.8 

1.138 

4 

22.69 

21.40 

6.1 

1.061 

5 

16.28 

15.84 

2.8 

1.028 

6 

66.50 

62.15 

7.0 

1.070 

7 

10.22 

9.84 

3.9 

1.039 

8 

68.71 

61.94 

10.9 

1.109 


Table 2-7. Summary of the Values Stored in the KCOSWT 
Table for Titanium 

Material 
Form Index 
KK 

Material 

Material Type Index JJ = 1 Form Index 

($/lb purchase wt) KK 

Material Type Index JJ = 1 
($/lb purchase wt) 

1 

15.00 

21 

10.50 

2 

14.50 

51 

65.00 

3 

14.00 

61 

25.00 

4 

13.50 

71 

95.00 

5 

13.00 

81 

132.71 

11 

12.00 

91 

50.00 

12 

12.50 

92 

50.00 

13 

14.00 

93 

50.00 
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where is a complexity factor for the 
material and the panel configuration, a, isj 
a constant, and WTTPS represents the 
complete TPS system gross weight. At 
present an area of 22,000 square feet is as- 
sumed for the complete TPS system, and 
this area is, in turn, used to compute the 
complete TPS system gross weight based 
on the unitized panel weight. 

Included under nonrecurring costs are engi- 
neering design and development, tooling, 
ground test hardware, flight test article, 
and flight test spares costs. Within these 
categories, engineering design and develop- 
ment costs include the engineering, design, 
development, laboratory test, support activ- 
Figure 2-43. Program Cost Summary ities for subsystem development, and the 

cost of hardware required to support 

laboratory development testing and component qualification testing. Tooling costs in- 
clude initial tooling for subsystem fabrication and general structural tooling. Ground 
test hardware, flight test articles, and flight test spares reflect the production costs of 
all the hardware used in subsystem and combined subsystem development testing. En- 
gineering design and development (ED&D) costs are derived using an equation similar 
in form to that used to calculate the TFU: 

ED&D = C * (WTTPS) C 1 

where C ^ and C-^ are constants. Tooling costs are derived based on an equation of the form; 
Tooling = C 3 * WTTPS 

where Cg is a complexity factor representing tooling dollars per pound. Ground test 
hardware, flight test articles, and flight test spares are derived based on an equation 
of the form: 

Cost = C 4 * TFU 

where C4 is a complexity factor representing the equivalent ground test hardware, 
number of flight test articles, or equivalent flight test spares, respectively. 

Included under recurring production costs are sustaining engineering and tooling costs 
(included as part of the TFU cost), production article costs, and test article conversion 
cost. Production article costs are for general production articles, and test article con- 
version costs reflect the cost of conversion to operational configurations. Both of these 
costs are derived based on an equation of the form: 

Cost = C * TFU 
5 

where Cg is a complexity factor representing the number of production articles and the 
equivalent test articles, respectively. 
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Included under recurring operations cost are refurbishment costs which reflect the 
cost of vehicle servicing and spares, and making it ready for flight validation over a 
ten year life of operation. Refurbishment costs are derived using an equation of the 
form: 


Refurbishment = C * TFU 

6 

where C g is a complexity factor representing the equivalent flight articles. The various 
complexity factors required for the program cost summary were derived based on work 
done at the San Diego operation by the Economic Analysis Department during the course 
of study on several advanced technology vehicles including the space shuttle vehicles. 

2.7.9 K-TABLES. These tables provide storage for values relating to the manufac- 
turing usage variance factor, manufacturing process setup time, manufacturing process 
runtime, material unit cost, and realization, labor rate, and overhead ratios, respec- 
tively. Although these tables currently contain constant values, they will later be 
replaced with equations, located in the COST subroutine, that will calculate the neces- 
sary data. 

KMUV (KK, JJ) refers to a table of manufacturing usage variance factors, which are 
called out by the material form index KK and the material type index JJ. 

Table 2-5 summarizes the presently available material forms referenced within the 
K-TABLES and the corresponding material form index KK. The material type index 
JJ is derived within the program as a function of the material type. At present, the 
KMUV table stores a constant value of 1. 10 for all values of KK and JJ, as discussed 
previously. 

KSETUP (KK, KS) refers to a table of manufacturing process setup time values in hours 
per pound actual weight (ACWT). These are called out by the material form index KK 
and the manufacturing operation index KS, which is set equal to JJ where JJ is a func- 
tion of material type. Values of the index KS for the various manufacturing operations 
required to produce a given part in titanium have been summarized previously (Section 
2.7.3). Once the initial value of KB is established, the corresponding setup time is 
picked from the KSETUP table and the standard hour calculation for the cutoff operation 
is made, as discussed previously. Then the value of KS is incremented and a new setup 
time is taken from the table and the standard hour calculation is made for milling. This 
procedure is continued until each manufacturing operation is considered. Table 2-8 
summarizes the values stored in the KSETUP table for titanium: zero is stored when 

a given operation is not required for a particular part. 

KRUN (KK, KR) refers to a table of manufacturing process run time values in hours 
per pound actual weight (ACWT) , which are called out try the material form index KK 
and the manufacturing operation index KR, which is analogous to the value of KS. 
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Table 2-8. Summary of the Values Stored in the KSETUP Table for Titanium 


Material 
Form Index 
KK 

Manufacturing Process Index KS* 

1 

2 3 4 5 6 

Standard Hours per Pound Actual Weight 

7 

1 

0.0050 

0 

0.0120 

0.0200 

0.0100 

0.0120 

0.0080 

2 

0.0060 

0 

0 

0.0250 

0.0100 

0.0140 

0.0080 

3 

0.0070 

0 

0.0160 

0.0300 

0.0100 

0.0160 

0.0090 

4 

0.0080 

0 

0.0180 

0.0350 

0.0100 

0.0180 

0.0090 

5 

0.0100 

0 

0.0200 

0.0400 

0.0100 

0.0200 

0.0100 

11 

0.0160 

0.1000 

0 

0.0500 

0 

0.0400 

0.0140 

12 

0.0200 

0.1100 

0 

0.0700 

0 

0.0600 

0.0160 

13 

0.0240 

0. 1200 

0 

0.0900 

0 

0.0800 

0.0180 

21 

0.0160 

0.1000 

0 

0.0500 

0 

0.0400 

0.0140 

51 

0.0300 

0.2000 

0 

0 

0 

0 

0.0200 

61 

0.0010 

0 

0 

0 

0 

0 

0 

71 

0.0010 

0 

0 

0 

0 

0 

0 

81 

0 

0 

0 

0 

0 

0 

0 

91 

0.0080 

0 

0 

0 

0 

0 

0 

92 

0.0010 

0 

0 

0 

0 

0 

0 

93 

0.0010 

0 

0 

0 

0 

0 

0 


* Refer to Table 2-4. 


Section 2.7.3 gives values of KR for the various manufacturing operations, and Table 
2-9 summarizes the values stored in the KRUN table for titanium. The index KR is 
initialized and incremented in the same manner as KS. For each cycle, a rim time in 
standard hours is calculated for a given operation using a value from the KRUN table 
and the procedure discussed previously. 

KCOSWT (KK, JJ) refers to a table of material costs in dollars per pound purchase 
weight (MAWT), which are called out by the material form index KK and the material 
type index JJ. Table 2-7 summarizes the values stored in the KCOSWT table for 
titanium. The calculation of material costs using values of KCOSWT is discussed in 
the section labeled MATERIAL. 
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Table 2-9. Summary of the Values Stored in the KRUN Table for Titanium 


Material 
Form Index 


1 


KK 


Manufacturing Process Index KR* 

2 3 4 5 6 7 

Standard Hours per Pound Actual Weight 


1 

0.0500 

0 

0.1200 

2 

0.0600 

0 

0 

3 

0.0700 

0 

0. 1600 

4 

0.0800 

0 

0.1800 

5 

0.1000 

0 

0.2000 

11 

0.1600 

1.000 

0 

12 

0.2000 

1.100 

0 

13 

0.2400 

1.200 

0 

21 

0.1600 

1.000 

0 

51 

0.3000 

2.000 

0 

61 

0.0100 

0 

0 

71 

0.0100 

0 

0 

81 

0 

0 

0 

91 

0.0800 

0 

0 

92 

0.0800 

0 

0 

93 

0.0800 

0 

0 


0.2000 

0.1000 

0.1200 

0.0800 

0.2500 

0. 1000 

0.1400 

0.0800 

0.3000 

0.1000 

0.1600 

0.0900 

0.3500 

0.1000 

0.1800 

0.0900 

0.4000 

0.1000 

0.2000 

0.1000 

0.5000 

0 

0.4000 

0. 1400 

0.7000 

0 

0.6000 

0.1600 

0.9000 

0 

0.8000 

0.1800 

0.5000 

0 

0.4000 

0. 1400 

0 

0 

0 

0.2000 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


* Refer to Table 2-4. 


KCC (KC, X) refers to a table of realization factors (X = KF), labor rates (X = KL), 
and overhead ratios (X = KV), which are called out by the cost center or operation 
index KC, and the factor index KF, KL, or KV for realization, labor rate, or over- 
head ratio, respectively. The cost center index KC is effectively a manufacturing 
process index and has values from one to seven, each corresponding to one of the 
production operations. Table 2-10 summarizes the values stored in the KCC table. 

2.8 INPUT, OUTPUT, AND MATERIAL PROPERTIES 

The input and output subroutines provide communication between the program user and 
the computer. Both aspects are discussed in detail in the User’s Manual. Input para- 
meters can be described in a number of broad categories: 


2-91 







Table 2-10. Summary of the Values Stored in the KCC Table 


Operation 
Index KC 

Factor Index 

Realization 

KF 

Labor Rate 
KL ($) 

Overhead Ratio 
KV 

1 

0.40 

4.75 

1.75 

2 

0.40 

4.75 

1.75 

3 

0.40 

4.75 

1.75 

4 

0.40 

4.75 

1.75 

5 

0.40 

4.75 

1.75 

6 

0.40 

4.75 

1.75 

7 

0.40 

4.75 

1.75 


a. Analysis indices which indicate the user’s preference for particular options in the 
analysis such as heating prediction method, panel configuration, the choice of one- 
or two-dimensional conduction, etc. 

b. Panel configuration, supporting structure, and vehicle geometry to define the para- 
meters necessary to perform thermodynamic, stress, and sonic fatigue analyses 
as well as to determine the weight and cost of the resulting TPS. 

c. Trajectory data (altitude, velocity, and angles of attack and sideslip as a function 
of time) or heating rates and panel pressure. 

d. Overall interchange factors for radiation between nodes of the thermodynamic 
configuration. 

e. Material property data. 

Input parameters for the first three categories of data are well explained in the User's 
Manual; the computation of overall interchange factors can be performed externally by 
any number of suitable programs available in the industry which take into account direct 
and indirect radiation exchange between all combinations of components. The remainder 
of this chapter will be devoted to presenting a summary of material property data which 
are applicable to the investigation of metallic heat shield panels. The thermodynamic 
and mechanical properties needed for the TPS design are shown in Table 2-11. 

Candidate metallic heat shield materials applicable to this study are given in Table 
2-12 along with estimated temperature and heat transfer rate limits for their use. 

The material properties for the first six of Table 2-11 are presented in plots 2-44 
through 2-49 as functions of temperature. A less extensive collection of material 








Table 2-11. Thermomechanical Properties 


Item 

Number 

Property 

1 

Thermal conductivity 

2 

Heat capacity 

3 

Modules of elasticity 

4 

Coefficient of thermal expansion 

5 

Yield strength 

6 

Ultimate tensile strength 

7 

Larson-Miller data for creep deformation 

8 

S-N curves for fatigue evaluation 


Table 2-12. Candidate Heat Shield Materials 


Maximum Use 

Temperature Heat Flux 


Metallic 

Titanium Alloys 

6 A1 - 4 V; 8 A1 - 1 Mo - 1 V 

To 1000 °F, 

1.6 

Btu/ft 2 sec 

Nickel Alloys 

Rene 41; Inconel 718; Inconel 625 

To 1700 °F 

9.0 

Btu/ft^ sec 

Cobalt Alloys 

L-605; Haynes - 188 

To 2000°F, 

15.0 

2 

Btu/ft sec 

Dispersioned Strengthened 
TD Ni; TD Ni Cr 

To 2200°F, 

20.0 

Btu/ft 2 sec 

Columbium Alloy 
Cb - 752 

To 2500 °F, 

35.0 

/ 2 

Btu/ft sec 

Tantalum Alloys 
T-222; Ta-lOW 

To 2800°F, 

60.0 

, 2 

Btu/ft sec 
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SPECIFIC HEAT, BTU/lb 


50 



TEMPERATURE. °F 

Figure 2-44. Thermal Conductivity for Candidate Metals for Heat Shield Applications 



TEMPERATURE, °F 


Figure 2-45. Specific Heat for Candidate Metals for Heat Shield Applications 
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Figure 2-46. Modulus of Elasticity for Candidate Metals for Heat Shield Applications 
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Figure 2-47. Coefficient of Thermal Expansion for Candidate Metals for 
Heat Shield Applications 
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ULTIMATE TENSILE STRENGTH, ksi 
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Figure 2-48. Yield Strength for Candidate Metals for Heat Shield Applications 



0 500 1000 1500 2000 2500 3000 3500 


TEMPERATURE, °F 

Figure 2-49. Ultimate Tensile Strength for Candidate Metals for Heat 
Shield Applications 
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properties for Larson-Miller creep deformation data and S-N fatigue data is currently 
available in the literature. A partial summary of creep data (Table 2-13) is given in 
Figures 2-50 through 2-66. 


Table 2-13. Larson-Miller Creep Rupture Data 


2219-T6 Aluminum 

(Annealed) Titanium, Ti-SAl-2.5 Sn 

Titanium, Ti-8Al-lMo-lV (Mill Annealed) 

718 Nickel, Aged 
718 Nickel, 20% CW + Aged 
Rene 41, ST + Aged 
L-605 Cobalt, Annealed 

Columbium, Cb-5V-5Mo-lZr (B-66), Stress Relieved, RT to 3000°F, 
Strength at Temperature 

Columbium, Cb-5V-5Mo-lZr (B-66) Recrystallized, RT to 3000° F 
Strength at Temperature 

Molybdenum, Mo-0. 5 Ti-0. 8Zr-0.03C (TZM) Stress Relieved Sheet, RT 
to 3000 °F, Strength at Temperature 

Molybdenum, Mo-0. 5Ti-0. 8Zr-0.03C (TZM) Recrystallized Sheet, RT to 
3000° F, Strength at Temperature 

Molybdenum, Mo-0. 5Ti-0. 8Zr-0.03C (TZM) Coated Sheet (Si, Cr, B 
Pack Cementation), RT to 2700 °F, Strength at Temperature 

Tantalum, Ta-8W-2Hf (T-lll), Stress Relieved, RTto 3000°F, Strength 
at Temperature 

Tantalum, Ta-8W-2Hf (T-lll) Recrystallized, RT to 3000°F, Strength at 
Temperature 


S-N data for fatigue analyses are even more difficult to locate than are creep results. 
For this report, two standard reports are cited (Reference 38 and 39). Some of these 
data which are applicable to the present study are presented herein in Figures 2-67 
through 2-74. The materials included are Rene 41 for a number of different heat 
treatments, and a number of different titanium alloys. Under Contract NAS8-27017 
with MSFC in Huntsville , Convair Aerospace is currently conducting a literature 
survey under a contract to perform fatigue evaluation of thermal protection systems. 
The result of this survey is a compilation and placement in a compatible format of 
available S-N information at various temperatures for the following materials : 
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Columbium (up to 2500 °F) 

Titanium (up to 650 °F) 

Rene 41 (up to 1600°F) 

Haynes 25 (up to 1900 °F) 

Boron composites j Da ^ a on whatever matrices and temperatures are available. 
Carbon composites ) 


The S-N data are to be presented: 

a. Tabular form, showing fatigue strengths at three values (minimum) of stress 
cycles (N). The 3NPs selected shall be a minimum value, a maximum value, 
and an average value. 

b. The tabulated data shall also be shown, as may be applicable, in graph form as 
S-N curves (showing effects of stress raisers), Goodman diagrams, plots of 
alternating vs. mean stress, etc. 

As a further task under the fatigue evaluation contract, the tables and graphs will show 
explicitly gaps in the data required for the space shuttle vehicle. As these new experi- 
mental data are developed and tabulated, they will be made available for the present 
study of TPS optimization. It is further anticipated that data sources for other TPS 
materials (both metallic and non-metallic) will be recognized so that available data can 
be presented although no new data will be developed. 

Thermodynamic properties for a number of insulation materials applicable to the TPS 
sizing program are included in Table 2-14 as a function of temperature. No mechanical 
properties are necessary since the insulation does not transfer a mechanical load to 
the structure. 
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Table 2-14. Insulation Properties 


Insulation 

Mean 

Temperature 

(°F) 

Thermal 
Conductivity, 
(Btu - in/hr -ft 2 °F) 

Specific 

Heat, 

(Btu/lb) 

pk, 

(lb/ft^) 
(Btu - in/ft 1 

Fibrous 





Microquarts 

0 

0.18 

0.170 

0.635 

3.5 pcf 

400 

0.336 

0.229 

1 • 18 


800 

0.515 

0.261 

1.80 


1200 

0.765 

0.278 

2.68 


1600 

1.22 

0.288 

4.28 

Dynaflex 

200 

0.31 

0.195 

1.24 

4 pcf 

400 

0.45 

0.215 

1.80 


800 

0.82 

0.243 

3.28 


1200 

1.34 

0.258 

5.36 


1600 

2.08 

0.268 

8.30 


2000 

3.09 

0.272 

12.40 

Dynaflex 

200 

0.29 

0.195 

1.74 

6 pcf 

400 

0.41 

0.215 

2.46 


800 

0.71 

0.243 

4.25 


1200 

1.09 

0.258 

6.55 


1600 

1.61 

0.268 

9.65 


2000 

2.31 

0.272 

13.82 

Dynaflex 

200 

0.27 

0.195 

3.24 

12 pcf 

400 

0.38 

0.215 

4.55 


800 

0.60 

0.243 

7.20 


1200 

0.85 

0.258 

10.20 


1600 

1.16 

0.268 

13.90 


2000 

1.53 

0.272 

18.33 

Zircar 

500 

0.60 

0.132 

7.2 

12 pcf 

1000 

0.75 

0.143 

9.0 


1500 

0.90 

0.149 

10.8 


2000 

1.25 

0.154 

15.0 


2500 

1.65 

0.144 

19.8 

Packaged 

75 

0.234 

0.203 

2.80 

Powder 

200 

0.238 

0.226 

2.85 


400 

0.254 

0.242 

3.06 

ADL-17 

800 

0.320 

0.258 

3.84 

12 pcf 

1200 

0.435 

0.265 

5.21 


1600 

0.588 

0.269 

7.04 
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Table 2-14. Insulation Properties, Contd 


Mean 

Thermal 

Specific 

pk. 

Insulation Temperature 

Conductivity, 

Heat, 

(lb/ft3) 

<°f) 

(Btu - in/hr -ft 2 °F) 

(Btu/lb) 

(Btu - in/ft 2 hr°F) 


Rigidized 


Silica 

200 

0.41 

0.200 

6.35 

(LI-1500) 

400 

0.46 

0.228 

7.10 

15.5 pcf 

800 

0.57 

0.264 

8.85 


1200 

0.87 

0.277 

13.40 


1600 

1.18 

0.289 

18.20 

Zirconia 

500 

0.55 

0.132 

13.80 

25 pcf 

1200 

0.67 

0.145 

16.65 


1600 

0.88 

0.150 

21.90 


2000 

1.18 

0.154 

29.50 

Foam 

Silica 

200 

0.57 

0.176 

5.70 

10 pcf 

500 

0.99 

0.210 

9.90 


800 

1.61 

0.241 

16.10 

Alumina 

200 

3.68 

0.220 

136.0 

37 pcf 

400 

3.00 

0.245 

111.0 


800 

2.36 

0.275 

87.3 


1200 

2.20 

0.295 

81.4 


1600 

2.24 

0.310 

82.9 


2000 

2.46 

0.320 

91.0 

Zironia 

200 

1.05 

0.122 

42.0 

40 pcf 

400 

1.1 

0.132 

44.0 


800 

1.32 

0.147 

52.8 


1200 

1.7 

0.153 

68.0 


1600 

2.28 

0.158 

91.2 


2000 

2.88 

0.160 

115.2 
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STRESS (ksi) TIME (hr) 



LARSON-MILLER PARAMETER 
P. = (460 + T) (20 + Log t) x lO -3 


Figure 2-50. Creep Data for 2219-T6 Aluminum 
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LARSON-MILLER PARAMETER 
P = (460 + T) (20 + Log t) x 10" 3 


Figure 2-51. Creep Data for (Annealed) Titanium, Ti-5Al-2.5 Sn 
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Figure 2-52. Creep Data for Titanium, Ti-8Al-lMo-lV (Mill Annealed) 
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Figure 2-53. Creep Data for 718 Nickel, Aged 
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LARSON-MILLER PARAMETER 
P = (460 + T) (20 + Log t) x 10" 3 


Figure 2-54. Creep Data for 718 Nickel, 20 Percent CW + Aged 
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Figure 2-55. Creep Data for Rene 41, ST + Aged 
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Figure 2-56. Creep Data for L-605 Cobalt, Annealed 
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P = (460 + T) (20 + Log t) x 10 _3 


Figure 2-57. Creep Data for Hastelloy X, Annealed 
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Figure 2-58. Creep Data for TD Nickel, As Rolled 
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Figure 2-59. Creep Data for B-66 Columbium, Recrystallized, Uncoated 
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LARSON-MILLER PARAMETER 
P = (460 + T) (20 + Log t) x IQ -3 


Figure 2-60. Creep Data for B-66 Columbium, Chromizing Corp 
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Figure 2-61. Creep Data for TZM Molybdenum, Stream Relieved, Uncoated 
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LARSON-MILLER PARAMETER 
P = (460 + T) (20 + Log t) x 10 -3 


Figure 2-62. Creep Data for TZM Molybdenum, Recrystallized, Uncoated 
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Figure 2-63. Creep Data for TZM Molybdenum, As Rolled, Uncoated 

2-114 



STRESS (psi) TIME (hr) 



LARSON-MILLER PARAMETER 

-3 

P = (460 + T) (20 + Log t) x 10 


Figure 2-64. Creep Data for TZM Molybdenum, Si, Cr, B Coated 
(Pack Cementation Process) 
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Figure 2-65. Creep Data for TZM Molybdenum, Disilicide Coated 
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Figure 2-66. Creep Data for T-lil Tantalum, Recrystallized, Uncoated 
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ALTERNATING STRESS (NO STEADY STRESS) (ksi) 




CYCLES TO FAILURE 


Figure 2-67. S-N Curves for Smooth Rene 41 Alloy Two Heat Treatments at Room 

Temperature, 1200°F, 1400°F, and 1600°Fwith Zero Steady Loads (A=oo) 
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Figure 2-70. S-N Curves for Smooth and Notched Specimens of Titanium 
Alloy, RC 55 Type 
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CYCLES TO FAILURE 

Figure 2-72. S-N Curves for Titanium Alloy RC-130B 



CYCLES TO FAILURE 


Figure 2-73. S-N Diagram for Ti-75 A Titanium Alloy Tested at Different 
Speeds With and Without Coolant 
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SECTION 3 


OPTIMIZATION STUDIES 


As part of Task 2 of this contract, a preliminary computer study was undertaken to 
demonstrate the application of the TPS sizing technique to optimization of the space 
shuttle TPS. In general, optimization procedures require evaluation of a performance 
index, which is to be extremized subject to constraints on the system. The constraints 
are reflected as penalties to the index. For the TPS, two performance indices not 
necessarily mutually compatible are considered: weight and cost. Under normal cir- 
cumstances, as system weight decreases, its sophistication and complexity in terms 
of materials, design, manufacturing, installation, inspection, and refurbishment 
often drive system cost up. Using either weight or cost as the performance index, 
a number of control parameters must be considered to evaluate how well the TPS 
functions. These are the constraints on the system; they determine how well the TPS 
protects the vehicle from a hostile aerothermodynamic environment. The insulated 
cover panels provide the aerodynamic surfaces that allow the vehicle to fly by trans- 
ferring aerodynamic loads from the surface to the load-bearing structure. Aerody- 
namic heating elevates the temperature of both the cover panels and the underlying 
structure, a behavior that decreases a materials' capability to handle these loads. 

In addition to the static loads caused by both mechanical and thermal stresses, the 
TPS must withstand dynamic loads produced primarily by engine noise and vibration, 
the turbulent boundary layer, and other sources. Here again performance is meas- 
ured by how well the system can withstand the adverse effects of this environment. 

The variables that determine or "drive" the performance of the TPS are the environ- 
ments themselves, heating rates, static mechanical loads due to aerodynamic pres- 
sure and thermal stresses, and the vibrational environment of acoustic excitation. 

The environments are affected primarily by the trajectory; i. e. , rates and integrated 
heat loads, pressure, and strength and duration of dynamic loads as influenced by 
engine thrust, duration of firing or running, and period of oscillation in the turbulent 
boundary. Other forcing functions such as panel flutter exist (an analysis has already 
been implemented in the next generation TPS sizing procedure), but they are not con- 
sidered in the present methodology. 

Once a given mission or trajectory has been selected and a certain area of the vehicle 
considered, the environments are essentially specified. To evaluate how well the 
system performs requires prediction of the response of the TPS to the environments. 
For the thermodynamic considerations alone, such an analysis is formidable since it 
requires solution of a finite difference statement of the energy equation over a period 
that may last 2500 seconds in real time for the orbiter. One of the simplest constraints 
on this analysis is to specify that some portion of the system, such as the primary 
structure, must not experience a temperature greater than some specified value. 
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If the system has only one component, the trajectory analysis is repeated for different 
thicknesses of this material until the variation of the primary structure temperature 
with thickness can be established. Then the TPS thickness that maintains the assigned 
temperature constraint can be established. For a one -dimensional heat transfer model, 
the problem is solved. However, problems of the shuttle TPS are more difficult and 
the solution is more complex. Flat plate panels , although good subjects for a one- 
dimensional heat transfer analysis, are inefficient in transferring aerodynamic loads. 
Therefore more complex configurations such as both open and closed corrugated panels 
are required. These in turn necessitate at least a two-dimensional heat transfer model 
to adequately predict the conduction and radiation between various elements of the panel. 
For the extremely high heating rates of the orbiter trajectory, the heat capacity of the 
metallic panel itself is an inefficient method of reducing heat transfer to the underlying 
primary structure. In the process of reducing heat transfer through the system, panel 
temperature is increased ’at a sacrifice in mechanical strength. Another extremely 
important factor in establishing metallic panel thickness is panel size, or more to the 
point, the distance between supporting elements of the structure. For the configura- 
tions considered in this contract, in which isolated elements are analyzed by state-of- 
the-art technique as simply support beams under uniform loading, panel thickness is 
determined by the distance between supports. The distance between support dictates 
panel width; hence the problem grows more involved. In an attempt to minimize TPS 
weight locally per unit area, panel size, thickness, and insulation size must be varied 
throughout the trajectory and the performance index, weight, must be evaluated. 

The next consideration to be made concerning a true TPS optimization procedure is the 
configuration of the panel. Should it be an open corrugation similar to a sine wave, or 
should the waves have flat areas ? Should one consider a sheet stringer configuration 
with a flat surface but stiffened by some sort of corrugated back plate ? Or should the 
panel be an integrally stiffened one machined from a single piece of material ? The 
variations are endless. Next, consider the insulation material itself. What is its 
configuration? What is the value of the insulation emissivity? Finally comes 
the problem of material selection. Which of the many metallics should be con- 
sidered? It quickly becomes clear that the number of variations in panel geometry, 
configuration, and material as well as those of the underlying insulation becomes over- 
whelmingly large. The number of variables must be drastically reduced to realistically 
consider a computer program that must be performed in the short run times required. 
The amount of material property data that can be reasonably stored in a computer of 
large but finite core size quickly established the limitation of using only one panel 
material instead of letting the computer select from a number of possibilities. Next, 
the complexity of the discrete element stress analysis and the necessity of transform- 
ing general panel characteristics such as panel length, width, number, and shape of 
corrugations and the like into the finite elemental volumes of the model dictate the 
choice of one configuration per run. The resulting variations in parameters were 
obvious: panel and insulation thicknesses. Since varying the metallic panel thickness 
to change its heat capacity and thus either increase or decrease the temperature of 
the underlying primary structure (the usual thermodynamic constraint) is a rather 
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inefficient technique as opposed to changing insulation thickness, it was decided that 
the latter would be done to satisfy the thermodynamic analysis. The only variations 
that reasonably affect strength characteristics and the ability to transfer static loads 
are either panel thickness or size (length and width). It was determined that of these 
two, thickness is more desirable. Thus, for a particular trajectory and its associated 
heat and aerodynamic loadings, the panel size (length and width), its configuration, 
and its materials are fixed; only the thicknesses are varied. Changing the panel thick- 
ness will admittedly change temperature distribution in the panel cross-section. How- 
ever, for small changes in thickness, the resulting small variations in temperature 
should have little effect on material property degradation. Also, thermal stresses 
are a function of temperature gradients, not temperatures, and experience has shown 
again that for reasonable changes in panel thickness, such as would result from the 
intelligent choice of initial values made by an experienced TPS designer, gradients 
and the resulting thermal stresses do not vary drastically. 

The optimization procedure, then, is performed by the program user who varies all 
parameters of interest. For a given panel material and configuration, the usual pro- 
cedure is to first vary panel size. The user must realize that of the two dimensions, 
only panel length is significant in the stress analysis. This is the dimension normal 
to the panel cross-sections defined in Figure 2-18. Panel length determines the forces, 
hence the internal loads, on the isolated beam used in the stress analysis. For the 
actual case of sizing the TPS panel, its width will be determined by sizing the support- 
ing structure that handles the load on the panel edges parallel to the cross-section. 
However, independent analyses conducted at Convair Aerospace indicate that for the 
high damping ratios generated by the slip joints between panels , minimum gauges are 
the design point for both support rails and posts. Thus, to determine minimum weight 
per unit area for a specified configuration and material type, panel length will be 
varied for a number of panel widths , and the optimum panel dimensions will be those 
that give minimum weight. The best of a number of different materials for the same 
panel configuration is determined first by establishing optimum panel dimensions for 
each material, then comparing weights for each optimum size. For example, the 
minimum weight of 2 pounds per square foot of a 12 X 24-inch panel of material A will 
be compared to the minimum weight of 1. 9 pounds per square foot of an 18 x 36-inch 
panel of material B. In turn, the dimensions of the panels of materials A and B have 
been determined as minimum weights per unit area for variations in length and width 
of both panels. Each panel configuration can then be compared for a variety of ma- 
terials. The optimization procedure may also be performed using cost as the per- 
formance index. 

In this first generation computer program, the determination of panel and insulation 
thicknesses is performed by an elementary iterative technique which will be improved 
and whose speed will be increased as experience with the techniques are gained. It is 
possible that a mathematically rigorous optimization procedure can be adapted to the 
calculation of thicknesses. However, such an inclusion will require some rather 
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sophisticated logic at the expense of computer storage space. Improved iterations 
leading to this investigation will be confined to Newton-Raphson techniques. 

The location of the panel on the bottom centerline of the Convair Aerospace booster 
is shown in Figure 3-1 with an indication of the typical dimensions. Although the 
booster TPS is really a shell construction as indicated in Figure 3-2, the area was 
analyzed as a panel under the same assumptions as the orbiter TPS; i. e. , each panel 
is allowed free thermal expansion so that an isolated section can be properly considered 
as a simply supported beam. The configuration of the panel itself is the so-called 
’’semi-smooth" corrugation of Figure 3-3. The wave length of each corrugation is 6 
inches and the depth is 0.4 inch. For the purposes of this analysis, the corrugations 
were assumed cross -wise to the flow. Two panel widths were considered, one and 
two feet, and panel lengths were varied from one inch to two feet (Figure 3-4). Since 
the booster TPS is uninsulated, the mathematical model showing the nodal geometries 
of the panel and the underlying aluminum cryogenic tank is given in Figure 3-5. For 
simplicity, the radiation view factors between nodes of the panel and the underlying 
tank wall were taken as unity; hence, each node of the corrugation sees only the node 
of the aluminum tank lying immediately below it. (As part of a future study, the 
radiation interchange factors of this configuration will be more rigorously defined, 
and the differences effected by this assumption will be assessed. ) The panel is Rene 
41, and the thermophysical properties used are shown in Figures 3-6 and 3-7. Param- 
eters were input as tabulated values as required by the computer program. The tra- 
jectory flown is given in Figure 3-8. The aerodynamic heating was computed for 
each run, using the plate -cylinder option for a non-uniform wall temperature distri- 
bution (i. e. , the flow field is computed as that of a wedge until the shock wave de- 
taches at high angle of attack; thereafter, heating is computed to a swept cylinder). 

In addition, heating rates to the TPS surface were modified to account for boundary 
layer separation and reattachment (in a conservative fashion) based on the data of 
Bertram (Reference 40). This distribution, used for both laminar and turbulent flow, 
is given in Figure 3-9. Input data required to perform the sonic fatigue analysis are 
summarized in Table 3-1, and input data for the weights /cost analysis are given in 
Table 3-2. 

The results of this investigation in terms of weight and cost per unit area are shown 
in Figures 3-10 and 3-11. For this particular minimum gauge supporting structure, 
with a panel width of two feet, the optimum panel length with respect to weight is 
six inches; theoretical first unit cost and both recurring and non-recurring production 
and operations costs (which are primarily weight driven) give the optimum length 
also at six inches (Figures 3-12 and 3-13). Manufacturing costs show optimum length 
at six to ten inches. The optimum spacing at such a small panel width (i. e. , six inches) 
for this configuration is not surprising; however, Figure 3-14 shows the results of an 
earlier study performed by Lockheed for the wing structure of a hypersonic cruise 
vehicle. For this case also, for a two-footplate, optimum spacing occurred at twelve 
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inches (Reference 41). The parameter presented in the Lockheed data is the effective 
thickness t defined by 

weight = (material density) x (panel area) x (effective thickness) 

Hence, weight per unit area is proportional to t. 




OF INTEREST 

Figure 3-1. Booster Panel Location 
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Figure 3-2. Body TPS Shell Structure 
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Figure 3-5. Mathematical Model 
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Figure 3-9. Heating Multiplication Factor Distribution 
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Table 3-1. Fatigue Analysis Input 


Period of turbulent boundary layer noise excitation 

100 sec 

Run length of turbulent boundary layer 

100 ft 

Local Reynolds number 

10 7 

Local velocity 

8000 ft/sec 

Local dynamic pressure 

1400 lb/ft 2 

Local Mach number 

6 

Period of rocket engine noise excitation 

40 sec 

Rocket engine thrust 

10 6 lb 

Rocket engine weight flow 

1000 lb/sec 

Rocket nozzle exit diameter 

12 ft 

Local sound speed 

1000 ft/sec 

Distance between point of interest and rocket engine exit 

20 ft 

Vehicle diameter 

40 ft 

Period of jet (flyback) engine noise 

20 sec 

Nozzle exit area 

3 ft 

Jet velocity 

1200 ft/sec 

J et engine weight flow 

200 lb/sec 

Vehicle velocity at flyback cruise 

650 ft/sec 

Engine thrust 

1000 lb 

Axial distance from point of interest to jet exit nozzle 

10 ft 

Radial distance from point of interest to jet exit nozzle 

10 ft 

Period of jet (flyback) engine scribbing noise 

10 sec 

Panel moment of inertia 

10 -3 in 4 

Panel length 

0.5 ft 

Panel width 

2 ft 

Panel modulus of elasticity 

10 7 lb/in 2 
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Table 3-2. Weight /Cost Data Input 


Panel type 

Panel length 

Panel width 

Panel length overlap 

Panel width overlap 

Distance between adjacent panels 

Support structure type 

Insulation thicknesses 

Insulation densities 

Number of corrugations across width 

Corrugation chord length 

Corrugation radius 

Skin density 

Weight of bolts 

Weight of nutplates 

Weight of washers 

Thickness of corner piece 

Thickness of post 

Density of corner post material 

Outside diameter of corner post 

Thickness of corner plates 

Thickness of support tube flange 

Thickness of support tube 

Density of support tube material 

Length of support tube 

Thickness of seal strip 


Corrugated 
6 in. 

24 in. 

1 in. 

1 in. 

1/2 in. 

1 

0 

0 

of panel 5 

2 in. 

6 in. 

0.300 ft-lb/in 3 
0. 400 lb/100 
0. 500 lb/100 
0. 100 lb/100 
0. 010 in. 

0. 020 in. 

0. 300 lb/in 3 
1/2 in. 

0. 010 in. 

0. 050 in. 

0. 020 in. 

0.300 lb/in 3 
0. 600 in. 

0. 010 in. 
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Table 3-2. Weight/ Cost Data Input, Contd 


Density of seal strip 

0.300 lb/in 3 

Thickness of long beam doubler channel 

6. 010 in. 

Thickness of corner doubler plate 

0. 010 in. 

Thickness of long beam 

0. 010 in. 

Thickness of short beam 

0. 010 in. 

Height of long and short beams 

0. 500 in. 

Widths of long and short beam 

1. 000 in. 

Density of beam material 

0.300 in. 



Figure 3-10. TPS Panel Unit Weight 
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THEORETICAL FIRST UNIT ($ X 10“ 2 /ft 2 ) MANUFACTURING COST ($/ft 2 ) 



PANEL LENGTH (inches) 


Figure 3-11. TPS Panel Unit Cost 



PANEL LENGTH (inches) 


Figure 3-12. TPS Theoretical First Unit Cost 
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Figure 3-14. Panel Size vs F for Corrugated Heat Shields 
with Hat Sections and Clip Support 
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SECTION 4 


SENSITIVITY STUDIES 


The prediction techniques used to design the TPS always have a degree of uncertainty 
caused by anything from the scatter of data used to correlate heat transfer prediction 
methods to the experimental error introduced while determining material properties. 
The system designer must account for these effects by either providing a margin of 
safety to allow for these contingencies or by limiting or refining mission constraints 
to avoid catastrophic failure due to the increased loads that may occur. The exper- 
ience of previous TPS design work indicates that the compounded conservatism intro- 
duced by simply adding the margins required to provide system safety can lead to 
significant weight and performance penalties of the vehicle. To combat these penalties 
in both the weight and cost of the TPS, two tasks must be performed: (1) the penalties 
due to the performance prediction uncertainties must be established, and (2) a rational 
method of combining these penalties must be established since many of the uncertain- 
ties and associated penalties are not independent. It is the purpose of a sensitivity 
study to establish the penalties in system weight and cost due the uncertainties in 
system design. 

Design uncertainties in the TPS itself can be separated into three categories: The 
environment, physical properties of the system itself, and operational and system 
characteristics. A partial list of these effects is given in Table 4-1. Parametric 
data gathered by establishing the weight and cost penalties due to variations in these 
uncertainties will establish the system sensitivities to these changes. 

Consider, for example, the weight of the TPS per unit area as the parameter whose 
sensitivity is to be established as a function of some uncertainty. First, the weight 
of the system can be determined for the nominal case. This analysis would include 
use of nominal aerodynamic heating and pressure prediction methods with a nominal 
boundary layer transition criterion over a smooth surface. Material and mechanical 
properties of the TPS would be assigned nominal values as would operations and 
system characteristics. For this investigation, the nominal case will be based on 
the optimum panel size for given TPS materials and configuration and a representa- 
tive trajectory. Next, an off nominal case will be established by a variation in some 
driving parameter. For example, for the case of turbulent flat plate heating shown 
by the data of Figure 4-1, a 3a variation in the data amounts to plus or minus 25 
percent on the nominal value predicted by the method of Spalding and Chi. The turbu- 
lent heating value predicted by this method is then varied by these factors of 1. 25 
and 0. 75 and new unit area weights are. deter mined. The change in weight from the 
nominal case can then be expressed as 
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where w denotes system weight and q-p denote turbulent heating. The change in the 
heating rate can be established most easily as a percentage increase or decrease in 
heating rate. The weight penalty and heating increase thus establish the parameter 


/ \ , the partial of weight with respect to heating method. 

\dg T / 


For variations in the 


heating method perturbations, the weight partial can be plotted as a function of these 
drivers. Parametric analysis of the various uncertainites outlined in Table 4-1 will 
then establish system sensitivities. 


Table 4-1. TPS Uncertainty Factors 

Environment 

Laminar heating 

Onset of transition to turbulent flow 

Onset of fully developed turbulent flow 

Turbulent heating 

Separated flow heating 

Local flow conditions 

Trajectory dispersions 

Roughness effects 

Vehicle aerodynamic attitude 

Local aerodynamic pressure loading 

Venting pressure 

Physical Properties 

Material properties (including conductivity, density, specific heat, 
surface emissivity, etc.) 

Material temperature limitations 
Mechanical properties 

Operations and System Characteristics 
Reentry trajectory corridor 
Attitude control system interactions 
Related systems and structure thermal limits 
Manufacturing, assembly and fabrication limitations 
Reuse constraints 
Inspection requirements 
Initial structural temperatures 
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Figure 4-1. Turbulent Flat Plate Heating Data 

The establishment of system weight penalties to a number of uncertainties which are 
mutually compatible is the subject of an extensive study to be undertaken in future 
work. 

As part of Task 3 , some preliminary computer studies were made to illustrate the 
application of the TPS sizing routines to determine the sensitivities of TPS cost and 
weight to various perturbations of input parameters. The TPS considered was the 
same semi-smooth configuration described in the optimization studies. The param- 
eter that was varied was the aerodynamic heating rate computed for the corrugated 
surface. The nominal case consisted of computing aerodynamic heating for the plate- 
wedge-cylinder option modified by the heating magnification factors of Figure 3-8. 

Two perturbations from this case were considered: the first was to assume that 
heating is unaffected by the corrugated surface. This means that the flow remains 
attached to the panel, and there is no heating amplification factor due to boundary 
layer reattachment. The second perturbation considered was a twenty-five percent 
increase in turbulent heating and a ten percent increase in laminar heating. This 
distribution is also shown in Figure 3-9. These heating distributions were applied 
to the case of the 24 X 24-inch panel. In addition, the nominal case and the case for 
uniform heating (with no magnification) were considered for the optimum panel size 
of 24 X 6 inches. 

For the non-optimum case of the 24 x 24-inch panel, the resulting weights and costs j 
are shown to be relatively insensitive to heating rates (Figures 4-2 and 4-3) whereas I 
the optimum panel design (24 x 6 inches) shows a rather weak sensitivity to varia- 
tions in heating rates. Temperature distributions are shown in Figures 4-4 and 4-5. 
Although no temperature constraint was applied to this particular TPS model, a limit i 
of 300 °F on the aluminum tank was met for both the nominal and uniform heating cases. 


4-3 




MAXIMUM HEATING AMPLIFICATION FACTOR 
Figure 4-2. Weight Sensitivity to Heating Amplification 



MAXIMUM HEATING AMPLIFICATION FACTOR 
Figure 4-3. Cost Sensitivity to Heating Amplification 
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PEAK TEMPERATURE (°R) 



Figure 4-5. Peak Tank Temperature 
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SECTION 5 

PROGRAM DEMONSTRATION AND DOCUMENTATION 

Running the TPS sizing computer program (a part of Task 4 of this contract) on the 
Manned Spacecraft Center's Univac 1108 proved to be a straightforward task. The 
major differences between the programs developed for the CDC 6400 and the Univac 
1108 at Convair Aerospace and MSC lay in control cards necessary to set up the runs 
and a few Fortran statements. Essentially, any Fortran IV program can be run in 
Fortran V. A list of exceptions is given in Table 5-1. The only other important ex- 
ceptions concern reading an end of file in the data and the maximum size of a number 
stored in the memory. 

For the case of reading an end of file, the Fortran IV statements for the CDC 6400 
operation are 

READ (Unit, f) 

IF (EOF, Unit) 1 1 , l 

The equivalent Fortran V statement is 

READ (Unit f , ERR = 1 , END - 1 ) 

u JL 


where 

Unit — represents the logical unit involved in the input/output transmissions 
f — represents format references 

ll, V l 3 — are statement labels 

In the Fortran IV statement, program control passes to statement 1^ if an end of file is 
experienced; otherwise , I 2 . In the Fortran V statement, control passes to I 3 if an error 
is detected in execution of the input or to 1 ^ if an end of file card is encountered. 

For the case of the size of a number which can be stored in the MSC 1108, the largest 
value is E30 (10 30 ). 
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Table 5-1. Fortran V Programming Reminders 


1. Do not use parameter names of more than six characters 

2. In computed GO TO statements, the comma after the parentheses is necessary, 
i.e. , GO TO (1, 2, 3, 4), K 

3. Do not use too many parentheses in READ and WRITE statements because the 
MSC 1108 gets confused. Use 

DO 100 J = 1, NSEG 

100 READ (5, 8) ((MAT (I, J), TAMP (I, J)), 1 = 1, K) 
rather than 

READ (5, 8) (((MAT (I,J), TAMP (I, J)), I = 1, K), J= 1, NSEG) 

4. The maximum word size on the MSC 1108 is six as opposed to the CDC's ten. 
Thus, use 10A6 instead of 6A10 

5. Use only Hollerith fields and not *’s in format statements 

6. Do not use multiple statements or multiple assignments on one card, i.e. , 

A=B$C=D$F=G or A = B = C = D = 0.2 

7. Avoid variable names that are the same as MSC library functions, i.e. , be a 
little selective in names. Do not use the word INPUT, STOP, WRITE. 

8. Do not use Convair Aerospace special functions or Fortran n library routines, 
e.g., ATANF(A) 

9. Do not use the two branch logical IF, e.g. , (X.GE.Y) 13, 14 

10. Do not use EXIT but rather call STOP 

11. In READ and WRITE statements, use 

READ (5, 10) X, Y, Z 
and do not use 


READ 10, X, Y, Z 





SECTION 6 


RECOMMENDATIONS AND PROGRAM REFINEMENTS 

The objective of these recommendations, which are a part of Task 5 of this contract, is 
to improve and update the computational methods used in the optimization of thermal 
protection systems (TPS) for the space shuttle vehicle and to generate weight and cost 
data for optimum TPS shuttle designs. 

The end product will be an improved TPS optimization computer program for application 
to the space shuttle . This program will employ refined numerical and mathematical 
models to decrease computer run time and make application of the computer program 
for design optimization easier for the program user. User's manuals, in addition to 
the operational Fortran program, will be delivered to MSC for dissemination to interest- 
ed users. 

The computer program will contain the latest available cost and weight data for shuttle 
TPS. This data, along with current TPS designs, will be used to perform detailed opti- 
mization studies and generate sensitivity coefficients in terms of cost, weight, and 
performance . 

This investigation is divided into four tasks. The first two update the computer pro- 
gram to include new TPS concepts and to improve computational speed. The third, and 
most important, is the design of optimum TPS at local areas on the shuttle vehicle. 

The fourth task is documentation of the investigation and the improved computer program. 

6 . 1 COMPUTER PROGRAM IMPROVEMENTS 

The primary function of this task is to incorporate refinements in the mathematical and 
numerical models of the computer code to either increase computational speed and ac- 
curacy or allow easier and more valid comparisons of the weight and cost of candidate 
TPS materials and configurations. In addition, new current TPS concepts on the space 
shuttle will be included in the sizing program. These include surface insulation, carbon- 
carbon systems, mass transfer cooling, and ablators. Figure 6-1 illustrates the current 
arrangement of the TPS sizing routine. 

6.1.1 MATHEMATICAL AND NUMERICAL METHODS . To speed and simplify compu- 
ter running times and output results , these methods will be incorporated into the appro- 
priate analysis subroutines which describe thermodynamic (aerodynamic heating, inter- 
nal heat transfer, and structural temperature response), stress, panel and support 
fatigue, and weight/cost analyses. 
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Figure 6-1. TPS Design Computer Program 


6. 1.1.1 Thermodynamic Analysis . 

Aeroheating . Typical improvements will include development of real gas oblique and 
conical shock theories at high angle of attack which are suitable for quick, easy, and 
economical computer usage. Also, real air properties will be added to the boundary 
layer calculations . Suitable curve fits (a Mollier chart for air) are available in com- 
puter programs existing at Convair Aerospace . The boundary layer edge properties 
necessary for any transition rationale are also available . Pressure methods them- 
selves will be the subject of a short, concerted investigation and evaluation effort. 

Implicit Heat Transfer . The computation of TPS heat transfer and temperature re- 
sponse in one or two dimensions with internal radiation and in both steady-state and 
time -dependent analysis will be studied to speed computer time. Several alternative 
heat transfer routines are available at Convair Aerospace in a number of heat transfer 
program and thermal analyzers. These consist of Program 2162 using the Crank - 
Nicholson forward-backward differencing scheme; Program 4560 using a similar tech- 
nique; and Program 2772 using the Runge-Kutta-Gill method to solve the differential 
equations. 

It is proposed that two or three typical TPS panels (for example a one -dimensional 
plate, a two-dimensional honeycomb, or a two-dimensional beaded panel) be input to 
each existing computer program and the computational times over a realistic portion 
of a trajectory be compared. A representative sampling of typical problems will select 
the best overall technique. This technique will be adapted to TPSOPT, either as an op- 
tion to the explicit technique or as a replacement. 
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6. 1.1. 2 Structural Analysis . The major changes and additions proposed for the TPS 
stress analysis routines are: 

a. Improvements and extension of existing analyses of the metallic reradiative TPS. 

b. The inclusion of additional analyses for stress redistribution due to creep and 
panel deflection. 

c. The addition of mathematical models and routines for additional TPS for external 
insulation, carbon-carbon, ablators, and mass transfer cooling. 

Details of the work to be performed under each of the other headings are discussed in 
the following paragraphs. 

Panel Analysis and Support Structure . To increase the flexibility of the existing stress 
analysis procedure and to make it readily adaptable to additional configurations, sepa- 
rate subroutines to accomplish each major step in the analysis (e.g. , finite element 
breakdown, thermal stress, creep, deflection, buckling, static strength, and redesign 
procedure) will be organized. The program will be extended to provide additional rou- 
tines for the analysis of panel edges or support points and the support structure (i.e. , 
beams, posts, clips). This will consider bending, crippling, column instability, creep, 
and deflection. 

A simplification of the redesign procedure is possible for panels in which creep is 
found to be noncritical, as determined by an initial trajectory analysis. Subsequent 
iterations of the redesign procedure need be made only for static strength at the most 
critical trajectory point. 

Stress Redistribution Due to Accumulated Creep Strains . In a beam member subjected 
to bending and/or thermal stresses , the accumulated creep strains due to exposure at 
temperature tend to relax the peak stresses and, hence, modify the subsequent creep 
rates and total creep strain. This effect can readily be included in the program by 
combining the cumulative creep strains throughout the cross section with the thermal 
strains (a T) in the existing thermal stress analysis. This stress analysis with creep 
relief will be included as an optional alternative to the existing technique since the 
additional computer run time required for the method may not be desired for all cases. 

Deflection Analysis . A panel creep deflection analysis, combined with the addition of 
support deflections and the existing analyses for deflections due to pressure loading and 
temperature gradient, will allow prediction of deflection throughout the life of the TPS. 

A design constraint based on deflection will be included in the redesign procedure. 

The maximum allowable deflection is given by the following expression: 

ALLOW DEFLECTION = ± (c + KL) 

where c and K are input parameters and L is the panel width. 
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6. 1.1. 3 Dynamics Analysis 


Panel Flutter and Stability. Panel flutter is a self -excited instability of elastic panels 
in supersonic flow, during which the panel oscillates in a direction normal to its plane 
at a constant amplitude; flutter usually results in fatigue damage to the panel. Con- 
siderable theoretical and experimental research has been conducted to understand the 
phenomenon and to develop design criteria to prevent it (Reference 42). The design 
approach recommended by Lemley will be incorporated into the TPS sizing routine , 
TPSOPT. The minimum panel thickness required to preclude panel flutter within the 
flight envelope will thereby be established as a design constraint. The non-dimensional 


flutter parameter 0 = 



is the governing function for the flutter criterion 


and reflects the influence of Mach number (M), dynamic pressure (q), Young's modulus 
(E), and panel thickness (t ) and length 1 . The Mach number correlation parameter, 
f(M) , shown in Figure 6-2, has replaced the familiar compressibility parameter. 



Figure 6-2, Mach Number Correction Factor Versus Mach Number 
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j3 =^M 2 ~ 1 since theory requires a prohibitively large panel thickness to prevent 
flutter at low transonic Mach numbers when jS approaches zero. Figure 6-3 presents 
the quotient of the Mach number correlation factor and the local dynamic pressure as 
a function of trajectory altitude and local Mach number . (A typical shuttle boost tra- 
jectory is given for reference.) The critical flight condition generally occurs at 
[q/f(M)] , although other trajectory points may require investigation, particularly 

if the panels are hotter at other flight points . 

Finally, Figure 6-4 gives the minimum panel thickness, t B , for flutter-free flight in 
terms of dynamic pressure, Mach number correlation parameter, material property 
E, and panel geometry. These curves will be programmed for the TPS sizing routine 
to allow determination of a panel thickness t B in terms of vehicle trajectory and panel 
geometry. 



Figure 6-3 . Plot of q/f(M) Versus Mach Number 
With Parametric Variation in Altitude 
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Figure 6-4. Aerodynamic Parameter q/f(M) Versus i/w 
With Variation in Structural Parameter (t / X) 

r> 

Acoustic Fatigue Analysis . The coustic fatigue analysis, now part of the computer 
program TPSOPT, is considered by Convair Aerospace to be an excellent approach, 
within the defined scope of effort, to the estimation of the acoustically induced loads 
and fatigue life of representative TPS panel configurations. Panel configurations in- 
clude (1) isotropic, (2) beaded, (3) corrugated, (4) honeycomb, and (5) integrally 
stiffened plates . 

Convair Aerospace was recently awarded Contract NAS8-27017, "Fatigue Evaluation 
of Thermal Protection Systems," by the Marshall Space Flight Center. Among the 
objectives of the contract are three pertinent to the effort of a continuing TPS optimi 
zation study. 

a. Obtain and compile basic fatigue (S-N) data for six candidate materials for space 
shuttle thermal protection systems . 
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b. Evaluate fatigue life prediction methods applicable to space shuttle thermal pro- 
tection systems. 

c. Investigate the damping provided by thermal insulation. 

The new acoustic fatigue refinements for the TPS optimization study are described in 
the following paragraphs. 

Higher Panel Modes and Associated Stresses . As in the case for the fundamental mode , 
the maximum stress in any higher mode will occur at a clamped edge of the panel. 
Therefore , the maximum stress can be calculated for each higher mode by considering 
only the clamped edge and the adjacent section of the panel extending to the first mode 
line. This approach provides an effective panel segment whose elements all move in 
phase in response to a correlated pressure. By equating bending moments, it can be 
shown, for example, that a clamped, simply -supported segment can be replaced by an 


equivalent clamped-clamped segment whose 
Figure 6-5. 



PANELS THAT WILL PROVIDE IDENTICAL 
STRESS AT CLAMPED EDGE. 

Figure 6-5. Equivalent Panel Segments 
for Computation of Stress in 
Higher Panel Modes 


length is increased by 25%, as shown in 

Acoustic Pressure Correlation. A re- 
view will be made of pressure correla- 
tion lengths associated with the various 
noise sources and at representative lo- 
cations on the space shuttle vehicle . 

This study will determine the maximum 
number of modes and the upper frequency 
limits that will be required in the inves- 
tigation of higher panel modes . 

Fatigue Characteristics of Panel Support- 
ing Structure . For normalized loads, 
comparative stress levels in representa- 
tive panels and supports will be evaluated 
to determine whether, from the fatigue 
standpoint, one is significant with respect 
to the other. If there is significant cor- 
relation, an attempt will be made to arrive 
at general weighting factors that could be 
applied to fatigue life computations to 
ensure that prediction techniques will be 
nominal. 

Cumulative Damage Theory . From Con- 
tract NAS8-27017, if a significantly im- 
proved damage theory is derived it will 
be incorporated into the TPSOPT program. 
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Damping Evaluation . A review of representative TPS configurations will be made and 
their damping properties estimated. If the spread in damping ratios is small, an aver- 
age value will be assigned to the computer program. If there is a wide variation, it 
may be necessary to assign a nominal value for each type of configuration. 

Fatigue Data Compilation. All applicable S-N data obtained on Contract NAS8-27017 
will be included in the program data bank as it becomes available. 

6 . 1 . 1 .4 Weight/Cost Analysis . The weight and cost analysis portion of the program 
will be expanded in scope, and the existing techniques will be updated and refined to 
simplify program organization and to speed program flow. The overall weight/cost 
subroutine will be simplified to facilitate future additions and changes. The input re- 
quirements will be generalized to increase the flexibility of the subroutine by making 
it possible to combine the structural elements of various existing TPS configurations 
and, thus, be able to synthesize new TPS configurations within the program. (At 
present, only specific configurations may be called at the user's option.) This addi- 
tional capability will allow combining the better elements of several different concepts 
during optimization and sensitivity studies. 

Data storage core requirements will be greatly decreased through the elimination of 
tables. Standard hour data will be curve fit so that standard hours can then be calcu- 
lated directly and used immediately (along with labor rates, overhead ratios, and 
realization factors incorporated into the subroutine) to derive labor and indirect over- 
head costs without any interim storage or retrieval of data. Other tables will be com- 
bined to calculate the cost of the actual purchased material based on the total quantity. 
Consideration will be made based on a lot buy , reflecting minimum-buy penalties , 
quantity buy discounts, special mill charges, die costs, form complexity, and toler- 
ance requirements. The manufacturing usage variance factor, which is established by 
accounting practice, takes into consideration material and part overbuying, losses, 
surplus, spoilage, duplication, substitution, and charges. Material cost data will be 
updated and incorporated into the subroutine . 

One major advantage resulting from the existing weight/cost analysis is derived from 
the fundamental level of approach, using a unique detail part listing process within the 
program. A significant improvement can be made by updating the program in terms 
of detail design and detailed manufacturing and fabrication techniques. Data relative 
to detail part identification, including the latest clip and fastener concepts, will be 
developed and incorporated into the program. A study will be made of state-of-the- 
art manufacturing operations and procedures and assembly operations including diffu- 
sion bonding, electron beam welding, and brazing processes as applied to the various 
TPS concepts and materials . New TPS configurations under consideration will be 
analyzed and included in the subroutine . 

The program cost summary (Figure 6-6) will be expanded to include recently identified 
cost elements and cost data factors . 
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THERMAL PROTECTION SYSTEM 


COST SUMMARY 


COST ($M) 

THEORETICAL FIRST UNIT COSTS — TFU 9.628 

NON RECURRING COST 

ENGINEERING, DESIGN, & DEVELOPMENT 
TOOLING 

GROUND TEST HARDWARE 
FLIGHT TEST ARTICLES 
FLIGHT TEST SPACES 

TOTAL NONRECURRING TPS COST 

RECURRING PRODUCTION COST 

SUSTAINING ENGINEERING - INCLUDED IN TFU 

SUSTAINING TOOLING - INCLUDED IN TFU 

PRODUCTION ARTICLES (1) 9.628 

TEST ARTICLE CONVERSION 2.888 

TOTAL RECURRING PRODUCTION COST 12.516 

RECURRING OPERATIONS COST 

REFURBISHMENT 12.824 

TOTAL RECURRING OPERATIONS COST 12.824 

TOTAL THERMAL PROTECTION SYSTEM PROGRAM COSTS 191.149 

NUMBER OF PRODUCTION UNITS 1 

Figure 6-6. Example Output of a Program Cost Summary 

The theoretical first unit cost which is, in turn, used to estimate the other cost items, 
is the predicted production cost of the prototype article , including manufacturing plan- 
ning, fabrication, subassembly, sustaining engineering, sustaining tooling, quality 
control, materials, and subcontracted parts. The existing program uses a weight- 
driven functional relationship to calculate the TFU. This procedure lacks the depth of 
analysis required to study the sensitivities of design tradeoff choices in terms of 
structural materials and methodology, other than those that change TPS gross weight. 
It also is deficient in cases where cost is an inverse function of weight. Indeed, ad- 
vances in technology are producing components with increased specific strengths and 
greater structural efficiencies; however, such advances require increasingly exotic 
materials and fabrication complexities — hence , greater cost. The existing manufac- 
turing cost calculation procedure contains many cost elements making up the TFU cost, 
including planning, fabrication, subassembly, materials, and subcontracted parts. 

The existing manufacturing cost procedure will be expanded to cover engineering, 
tooling, and quality control considerations. The TFU cost will then be derived based 
on the manufacturing cost portion of the subroutine . In this way the remaining TFU- 
driven cost items will be made sensitive to design tradeoff studies made at the detail 
level , 

Addition of engineering costs requires the development of a method for predicting the 
necessary engineering manhours necessary and the related cost elements. Considera- 
tion will be made for such items as shop liaison, reliability, engineering changes and 

6-9 


53.009 

63.545 

23.588 

19.256 

6.451 

1 65, 849 



fixes, and technical and preliminary design required during the initial concept phases. 
The total engineering costs are obtained for the complete TPS by applying current labor 
rates and overhead ratios to the total labor hours required. The estimating procedure 
will be used based on the engineering tasks to be performed, and the data will be de- 
rived from Convair Aerospace labor accounting records. 

Refurbishment cost data is presently becoming available as a result of studies made 
within Convair Aerospace and at Lockheed and McDonnell-Douglas (References 43 and 
44). These studies estimate the tasks required and the relative reliability of the vari- 
ous TPS concepts and materials. This data will be used to refine the analysis of re- 
furbishment costs within the program, and better define complete system nonrecurring 
operations cost (both refurbishment and maintenance). The existing manufacturing cost 
methodology, based at the detail part level, will be extended to include an analysis of 
the refurbishment tasks necessary for each TPS concept. These tasks include inspec- 
tion, removal, re installation, and certification of TPS panels . Associated costs in- 
clude manhours required to perform the tasks, equipment, indirect overhead, and 
spares hardware. Since detail parts data are now available within the program, it is 
possible to relate each task requirement to specific operations necessary at the detail 
part level, and in this manner to derive the manhour and equipment requirements for 
the refurbishment phase based on the actual work to be done. Overhead and spares 
requirements can also be derived on this basis. Consideration will be incorporated 
into the analysis for the latest available data concerning the various material and sys- 
tem reliabilities. 

In this task, additional weights and costs data will be gathered from the literature 
(Reference 45) and prediction techniques will be developed and refined to include a 
number of both active and passive cooling systems . Prediction techniques will be 
developed in each system in sufficient depth to allow a meaningful cost and weight com- 
parison among the concepts and with the metallic reradiative system already described 
by the existing program. 

Material Property Data . A great deal of mechanical and thermodynamic data are cur- 
rently available in the literature concerning nominal values of density, thermal con- 
ductivity, thermal coefficient of expansion, modulus of elasticity, ultimate strength, 
and yield strength primarily as functions of temperature. Some pressure -dependent 
data are available for insulating materials. It is suggested that a computer routine be 
developed to interpolate values of physical properties which are functions of two vari- 
ables (temperature and pressure). 

Not only are thermodynamic and mechanical properties necessary, but also needed 
are the S-N curves of the fatigue analysis, and, no matter how rudimentary, the cost 
data for raw materials, machining, and fabrication processes. 
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Weight data will also be developed in a similar manner, especially for either metallics 
that oxidize (thus changing panel weights) or metallics that must be coated. Next, the 
typical nonmetallic panel materials should be investigated in the same depth and detail. 
Again, material, manufacturing, development, refurbishment, and similar cost data 
must also be cataloged for input to the sizing program. 

6.1.2 ADDITIONAL TPS DESIGNS . To evaluate TPS concepts that are truly repre- 
sentative of the state of the art in both performance prediction and hardware produc- 
tion, a number of additional configurations of interest and importance to space shuttle 
applications will be incorpbrated into the TPS sizing program. 

6.1.2. 1 Passive Cooling Systems 

Reusable Surface Insulation (RSI). The typical concept shown in Figure 6-7 is RSI 
material bonded to a metal substrate . A gap is shown between it and an internal struc- 
ture such as an aluminum tank wall. A 
supporting structure may also be included 
between substrate and primary structure . 
Other concepts (such as wing surface) 
can have the RSI directly bonded to an 
internal structural panel of titanium or 
other material. The existing thermal 
heating analysis will be used to predict 
the temperature histories and maximum 
temperature of the outer surface, the 
bond line, and the internal structure. 
However, wave band dependence of in- 
depth radiation effects should be con- 
Figure 6-7. RSI Concepts sidered. 

The primary difficulty in analyzing the RSI system is the stress analysis — in particu- 
lar , the deformation of the bond line . The approach to the stress analysis of the RSI 
will be two -pronged: first, finite element methods (for large deflections) will be in- 
vestigated to see if they can be simplified enough for economic analysis of the RSI; 
and, second, the contract monitor and program manager will confer early in the pro- 
gram to review the latest analytical techniques and their application to RSI. 

Within the scope of this contract, acoustic fatigue behavior of the RSI can be predicted 
by existing methods of the TPS sizing routine. Weight and cost data will be developed 
as discussed previously. 

Carbon-Carbon Composites . Oxidation inhibited carbon-carbon composites are being 
considered for use as reusable radiative heat shields for space shuttle vehicles. The 
selection of materials for the leading edge and lower surface areas (Figure 6-8) on 
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these vehicles presents a severe design 
problem based upon the long life and cy- 
clic exposure design criteria established 
for the multimission space shuttle. 

Carbon and graphite composites are 
unique in that they possess excellent 
high-temperature stability and retain 
structural integrity at temperatures above 
4000° F. This property alone makes the 
carbon-graphite class of materials unique 
to other types being considered as candi- 
date heat shield materials . In addition , 
they have low thermal expansion and can 
be tailored for mechanical behavior. Laminated carbon materials and unique 3-D re- 
inforced carbon-carbon integrated structures are being developed for heat shield 
applications. 



Figure 6-8. 


Space Shuttle High- 
Crossrange Orbiter 


One concept uses laminated carbon-carbon composite material as an integrally- 
stiffened shell (Figure 6-9). This idea is based on existing carbon-carbon technology 
developed over the last six to eight years and is severely limited in design flexibility 
by problems associated with the complex tooling required, local reinforcements for 
attachments, and overall reliability. Because of these problems, a 3-D reinforced 
carbon-carbon fluted core was selected for development (Figure 6-10). This is a 
lightweight, double-faced, truss core structure that integrates an external heat shield 
cover panel with a support structure which, in turn, can be attached to a primary load- 
carrying structure. The concept consists of a 3-D woven truss core fabric rigidized 
with high -carbon-yield resins, pyrolyzed, and then carbonized up to 4000° F in an inert 
atmosphere. A typical leading edge design using the carbon-carbon fluted core is also 
shown in Figure 6-10. For the lower surface area identified in the schematic of the 
shuttle vehicle (Figure 6-8) , candidate heat shields may include several types of fluted 
core constructions. Carbon-carbon leading edges and panels will be analyzed on the 
modified version of computer program TPSOPT. 


The structural analysis of the carbon-carbon system will be performed using the exist- 
ing discrete element analysis for the flat panel. The program will be investigated to 
determine applicability of existing mathematical models to the leading edge . 

Heat Sink Concepts. The potential benefits of using heat sink thermal control for the 
space shuttle structure have been shown in a recent study, "Space Shuttle Aluminum 
Booster Study" conducted by Convair Aerospace. 

Heat sink thermal protection will be subjected to design conditions and vehicle operat- j 
ing problems that prior studies have shown are the critical criteria in sizing and eval- j 
uating different design concepts. The thermal loads will size the amount of heat sink 
required once materials have been selected and operating temperature limits set. 
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Figure 6-9. Three-Dimensional Reinforced Carbon-Carbon Core Configurations 
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Thermal stresses in the heat sink tank 
structures using aluminum have been 
found to be a secondary problem be- 
cause the maximum allowable tempera- 
ture used has been 300°F or less and 
heat sinking has kept the temperature 
distribution sufficiently uniform that the 
E a AT product remains less than 4000 
psi. Similar results can be expected 
with other concepts where material stress 
allowable limits require low operating 
temperatures. 


Potential heat sink designs are included 
in two basic arrangements: heat sink 
integral with the airframe structure and 
heat sink separate from the structure . 
For the integral concept, the heat sink 
may be fabricated in the same part as 
the structures or it may be a separate panel or shingle bonded or brazed to the struc- 
ture. The inner surface may be open to ambient conditions, exposed to liquid oxygen, 
or insulated from exposure to liquid hydrogen with cryogenic insulation such as poly- 
phene lyne oxide foam. In case of the separate panel concept, the panel can also con- 
sist of jacketed nonstructural or phase-change materials. 

Table 6-1 lists candidate heat sink materials and material properties. These materi- 
als were selected for their good heat absorption capacities and represent the maximum 
capacity available in materials relevant to this study. 

Along with the geometry of the panel concept and the materials used, operating tem- 
perature is an important parameter to consider. Figure 6-11 indicates the influence 
of operating temperature variation on panel unit weight. The distance between the two 
curves accounts for variations in the heat absorption and emissivities of the different 
materials considered. The reduction in the panel weight with increase in outer wall 
temperature is caused by increased radiation at the higher temperatures. 

The thermo -structural, weight, and cost analyses for the various heat sink TPS con- 
cepts can be carried out using the existing techniques in the TPSOPT program. Some 
modifications must be made to accommodate the heat absorption due to phase change . 

In the transient conduction program, the nodes will be specified to absorb the heat of 
fusion at constant temperature of melting. Additional changes may be made for stress 
and some fatigue analysis. Typical results already developed at Convair Aerospace 
for the aluminum heat sink TPS are shown in Figure 6-12. 



Figure 6-10. Carbon-Carbon Leading 
Edge Designs 
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Table 6-1. Candidate Heat Sink Materials 



HEAT OF 

t melt/ 

DENSITY 
(lb. /In. 3 ) 

SPECIFIC 

THERMAL 

TOTAL HEAT ABSORPTION (Btu/Ib.) 

% VOL. 


MATERIAL 

FUSION 
(Btu/lb. ) 

t boil 

<*F) 

HEAT 
(Btu/lb. - 
•F) 

CONDUCT- 
IVITY (Btu/ 
lb. ft. -*F) ' 

-150* 
ro 250’ F 

70* 

TO 250* F 

70* 

TO 1, 000* F 

70* 

TO 2, 000* F 

CHANGE 

ON 

FUSION 

CORROSIVE- 

NESS 

Phase -Change: 












Lithium 

286 

354/ 

2,400 

0.019 

D 

22.5 

370 

180 

1,216 

2,216 

1.5 


Sodium 

49.5 ; 

208/ 

1,638 

0. 035 


40 

159 

103 

324 


2.6 


Water 

79.7 

32/ 

212 

0. 036 

n 

0.25 

320* 

130* 

— 

— 

-8.3 


Linear Polymer 
(PPO) 

60 

500/- 

0. 033 

0. 45 

0. 19 

166 

81 

478 

— 



Non-Phase 

Change: 












Aluminum 

170 

1,220/ 

3,740 

0. 099 

0.23 

120 

85 

41 

214 

| 

614" 

6.6 


Beryllium 

470 

2,340/ 

5,020 

0. 066 

0.45 

87 

166 

81 

418 

868 

— 


Graphite 

Composite 

— 

— 

0.061 

0.18 

100 

67 

32 

167 

347 

— 


JP-5 

— 

— 

0. 027 

0.60 


185 

90 

465 

965 

— 



• T » 200* F ** Includes heat of fusion 



Figure 6-11. Influence of Operating Temperature 
Variation Upon Panel Unit Weight 
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Figure 6-12 . LHg Tank Wall Temperature 
Vs . Tank Wall Thickness 


6. 1.2.2 Active Cooling Systems 


Regenerative Systems. A technique has been developed under IRAD funds to predict 
performance of a space radiator (Reference 46). This method will be investigated for 
adaptation to the specific applications of convective aeroheating to the shuttle vehicle . 
It may provide capability to evaluate the performance, weight, and cost of one type of 
active TPS to be compared with other systems already described by the computer 
program . 

Mass Transfer Cooling . When mass transfer is used as a means of cooling, there are 
at least two principal heat reducing mechanisms at work. First, heat is absorbed by 
the cooling fluid as it travels against the direction of heat flow from some reservoir 
to a surface of higher temperature and, consequently, lowers the wall temperature; 
second, as the fluid emerges from the surface of the wall, it forms an insulating layer 
between the surface of the wall and the hot gas . 

In one such method of cooling, the wall is manufactured from a porous material and 
the coolant is blown through the pores. The coolant film on the hot-gas side is, there- 
fore, continuously renewed and the cooling effectiveness can be made to stay constant 
along the surface. When a liquid is used as a coolant, a liquid film is created on the 
hot -gas side which is evaporated on its surface , and the heat absorbed by the evapora- 
tion process increases the cooling effectiveness. This process is frequently referred 
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COOLANT FLOW REQUIRED (LB/SEC-IN. 


to in the literature as evaporative -transpiration cooling. The ultimate criterion for 
evaluating coolants is minimization of the system weight required to maintain the wall 
at the desired temperature. Figure 6-13 gives coolant weights required to maintain 
the wall at specified temperatures for some coolants. 



Single- Phase Mass Transfer . Most of the 
analytical investigations on transpiration 
cooling have been concerned with the lami 
nar boundary layer - primarily a result 
of the laminar flow problem being more 
amenable to analytical treatment than the 
turbulent flow. Unfortunately, transpira- 
tion cooling is more likely to be used 
where turbulent flow and, consequently, 
higher heating rates are involved. Even 
if the flow is laminar initially, the dis- 
turbance resulting from the injection of 
the coolant tends to cause transition to 
turbulent flow. 

Transpiration cooling effectiveness will 
be predicted for this study by consider- 
ing the single-phase flow of either a 
liquid or gaseous coolant. The heat 
transfer analysis is for the one-dimen- 
sional case with constant material pro- 
perties . The temperature of the porous 
matrix T through which the coolant flows 
is described by the equation 


Figure 6-13. Coolant Weight Requirements 
Versus Wall Temperature 
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where 0 is a heat source or dissipation 
function. Under the same assumption, the temperature of the coolant, t, is given by 


„ at „ at 
pC P 57 + f"°p sc 


0 


where the conduction of the coolant has been neglected and the material properties are 
those of the coolant. For the case in which either the film heat transfer coefficient 
between the coolant and matrix is very large or the flow rate of the coolant is moder- 
ate , the matrix and coolant are in thermal equilibrium and T = t. Hence, the energy 
equations for the two materials may be combined to give 
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where the subscript m refers to the matrix of the TPS and c refers to the coolant. 

This equation is identical to the heat conduction equation already solved in the conduc- 
tion subroutine of program TPSOPT, with the exception of the heat sink term involv- 
ing the temperature -dependent properties of the single-phase coolant. The modifica- 
tion to the computer program is a minor one . Sizing the TPS will involve computing 
the weight and cost per unit area of the TPS as a function of total coolant volume re- 
flected by coolant velocity, u. 

For this preliminary investigation of the transpiration cooling phenomenon, the be- 
havior of the external coolant effectiveness is assumed predicted by the empirical 
relationship developed by Bartle and Leadon for turbulent flow (Reference 47). 

The coolant effectiveness for laminar flow will be handled in a similar manner; e.g. , 
the data correlation of Figure 6-14 . 
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Figure 6-14 . Laminar Coolant Effectivess 


Ablation Analysis . The computer pro- 
gram will be refined to incorporate at 
least one type of ablation analysis . The 
technique selected will be decided after 
investigation of available methods and 
by considering the most promising type 
of ablative materials for space shuttle 
application . 

Methods of analysis of ablative heat 
shield materials, can be categorized, 
depending on the material melting point 
and oxidation chemistry, as: 

a. Oxidation Controlled. When the 
melting temperature is greater 
than the radiation equilibrium 
temperature. 

b. Simple Sublimers. When the melting 
temperature is lower than the radia- 
tion equilibrium temperature. 

c . Pyrolytic Ablators . When the ma- 
terial decomposition into pyrolysis 
gas and char occurs in depth . 
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Theoreretical methods are available for detailed treatment of these cases. They 
are based on multicomponent mass injection and mass transfer correlations of the 
boundary layer coupled with the wall chemical kinetics. For the present study, the 
treatment for Case 1 is considered too sophisticated to incorporate in the computer 
program. 

It is felt that a less sophisticated method of analysis can be effectively used in the pro- 
posed study; the basic problem to be solved is locating the pyrolysis zone as a function 
of time. Two simplified methods that are candidates for use are discussed in the follow- 
ing paragraphs. 

Heat of Pyrolysis Technique . The degradation of an ablative system, and the attendant 
temperature distribution through the system, can be analyzed in a simplified manner 
by using a property termed heat of pyrolysis. Heat of pyrolysis is defined as the 
amount of heat required to degrade one pound of material and has the units of Btu/lb . 

The thermal model used with this technique will be composed of three zones. Starting 
from the heated surface, the first zone is comprised of char (degraded ablator). The 
second zone is the pyrolysis zone where the degradation of the ablator takes place . 

For purposes of analysis this zone can be considered to be of zero thickness. The 
third zone is the virgin ablator material. The thermal model requires a moving co- 
ordinate system to account for the increase in thickness of the char zone and the de- 
crease in thickness of the virgin material zone . 

At the end of the calculation interval, the amount of heat flowing from the heated bound- 
ary to the pyrolysis zone is used in conjunction with the heat of pyrolysis value to com- 
pute the amount of material degraded during the calculation interval. This, then, de- 
fines the location of the pyrolysis zone for the following calculation interval, and the 
conduction networks in the char and virgin material zones are redefined accordingly. 

Heat of Ablation Technique . A quantity known as heat of ablation can be used to per- 
form a simplified thermal analysis of an ablative TPS. Heat of ablation is a general- 
ized property that accounts for the mass loss from an ablative system as a function of 
heating rate; the units are BtuAb, 

The thermal model for this technique uses a boundary condition that moves in the abla- 
tor coordinate system. This boundary is the pyrolysis zone, which can be treated as 
a constant-temperature boundary. The constant temperature is the pyrolysis tempera- 
ture and is material dependent. 

For a calculation interval, the surface heating rate is used to determine the heat of 
ablation. This is then used to calculate the amount of material degraded, and from 
this is derived the location of the pyrolysis zone for the next calculation interval. A 
one -dimensional conduction solution is performed from the pyrolysis zone boundary 
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through the ablator/structural network. At the start of the next calculation interval, 
the conduction network in the ablator is redefined to account for the movement of the 
pyrolysis zone. , 

6.2 PROGRAM APPLICATIONS AND OPTIMIZATION TECHNIQUES 

This section of the program plan concerns mainly the inclusion of short-term computer 
coding improvements that will speed and facilitate the actual technique of optimizing 
the TPS panels and supporting structure at various places or on various areas of the 
space shuttle vehicle. For the sake of discussion in this program plan, such techniques 
are divided into two categories: short-term improvements and true optimization pro- 
cedures . 

6.2.1 SHORT-TERM IMPROVEMENTS . During the course of checking our computa- 
tional procedures and performing optimization and sensitivity studies , close attention 
will be paid to developing automatic methods to speed the overall assimilation of TPS 
design data and ease the program user's task in generating information that is mean- 
ingful and useful to the space effort. The program resulting from Contract NAS9-10956 
will size the panel and insulation thicknesses for given values of panel length, width, 
and materials, for given thermal and structural constraints, and for a given vehicle 
trajectory (or, alternatively, input heating rates and pressure distribution). For a 
given panel configuration (e.g. , Rene' 41 panel backed by Dynaflex, which is separated 
by an airgap from a titanium primary structure) and trajectory, the program user will 
run a series of cases in which support structure spacing and panel width and length are 
varied. The optimum design for this configuration will be the one that either weighs or 
costs the least. This can be established quite readily by plotting a curve of cost or 
weight as a function of the varying parameter (spacing or panel width or length). The 
design procedure is then repeated for different panel and support configurations, ma- 
terials, trajectories, vehicle locations, and whatever other variables are of interest. 

It will help the procedure considerably if subsequent cases of a particular configura- 
tion can be stacked and the corresponding weights and costs either plotted out automa- 
tically or stored to be curve fitted and the minimum determined analytically later at 
the end of the runs . The possibility of graphically presenting any and all data will be 
thoroughly investigated and adopted where it proves feasible. 

6.2.2 FORMAL OPTIMIZATION PROCEDURES . A great deal of information is cur- 
rently available at MSC concerning the application of mathematical optimization pro- 
cedures to the case of a one -dimensional TPS configuration consisting of up to three 
different materials. The payoff function is usually weight per unit area subject to 
temperature constraints at the material boundaries or interfaces. Early in this pro- 
posed investigation an independent survey will be made of the optimization procedures 
currently available at Houston (e.g. , adoptive creep, pattern search, method of 
steepest descent, Davidon method, etc.), and the possibility of their application to 
the complexities of the present problem will be assessed. It is altogether possible 


6-20 


that the application of one optimization procedure that works well for a wide variety 
of configurations could be implemented to the case in which the TPS weight can be 
varied by only two parameters: panel thickness and insulation thickness. 

6 . 3 OPTIMIZATION AND SENSITIVITY STUDIES 

Detailed optimization and sensitivity studies of those TPS concepts already developed 
in the program TPSOPT will be begun on the first day of the study. For representa- 
tive trajectories and areas of interest on the vehicle, parametric studies will be under- 
taken to establish weights and costs of optimum panel sizes and support member spac- 
ing for the variety of panel configurations (beaded, corrugated, honeycomb, integrally 
stiffened) and supporting structures already available in the program. 

The most promising candidates will be established initially over their expected range 
of application. Typical candidates include titanium alloys to 1000°F, nickel alloys to 
1700°F, cobalt alloys to 2000°F, dispersioned-strengthened such as TD Ni or TD NiCr 
to 2200°F, columbium to 2500°F, and tantalum alloys to 2800°F. One particular panel 
configuration and material will be studied in great detail to determine the effects of a 
variety of insulation materials and orientations. 

The results of these studies will be illustrated as carpet plots of weight and cost as a 
function of heating load, duration, configuration, materials, and vehicle location. 

Such data will then be in a form convenient for inclusion into a number of available 
synthesis programs being developed both at Convair Aerospace and NASA. 
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